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BROADBANDING  TECHNIQUES  FOR  MICROSTRIP  PATC'.  I  ANTENNAS  -  A  REVIEW 


K.C.  Gupta 

Department  of  Electrical  and  Computer  Engineering 
University  of  Colorado 
Boulder,  CO  80309-0425 


Abstract 

The  need  for  development  of  wideband  microstrip  patch  antennas  has  been  very 
well  recognized.  This  article  is  a  review  of  various  techniques  attempted  for  increasing  the 
bandwidth  of  microstrip  patches.  These  include  multiple  resonator  configurations  (with 
patches  located  on  the  same  plane  as  well  as  multiple-layered  substrate  geometries)  and 
wideband  impedance  matching  network  approach.  Various  results  reported  in  the 
literature  are  discussed  and  areas  for  continued  research  efforts  are  identified. 

1 .  INTRODUCTION 

Narrow  bandwidth  available  from  printed  microstrip  patches  has  been  recognized  as 
one  of  the  most  significant  factors  limiting  the  widespread  applications  of  this  class  01 
antennas.  The  operating  bandwidth  of  a  single  linearly  polarized  patch  antenna  is  limited 
by  its  input  VSWR  (or  reflection  loss)  and  is  inversely  proportional  to  the  Q-factor  ot  the 
patch  resonator.  Typical  values  for  the  fractional  bandwidth  are  2  to  4%.  Over  the  last 
decade,  several  attempts  have  been  made  for  improving  the  bandwidth  of  microstrip 
patches.  This  article  is  a  review  of  various  techniques  employed  for  broadband  design  of 
microstrip  antenna  elements  and  arrays  and  attempts  to  present  a  comparison  of  the 
results  reported. 
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Various  techniques  used  to  increase  the  bandwidth  of  microstrip  patch  antennas 
may  be  classified  as  follows:  (i)  Decreasing  the  Q-factor  of  the  patch  by  increasing  the 
substrate  height  and  lowering  the  dielectric  constant,  (ii)  Use  of  multiple  resonators 
located  in  one  plane,  (iii)  Use  of  multilayer  configurations  with  multiple  resonators  stacked 
vertically,  (iv)  Use  of  impedance  matching  networks. 

These  tectiniques  are  discussed  in  Sections  2  through  5.  Special  considerations 
for  the  bandwidth  of  circularly  polarized  patches  are  discussed  in  Section  6  Finally, 
bandwidth  considerations  for  arrays  of  microstrip  patches  are  contained  in  Section  7. 
Comparison  of  vanous  approaches  and  concluding  remarks  are  contained  in  the  last 
section. 

Bandwidth  Definitions 

Since  the  bandwidth  (BW)  of  linearly  polarized  single  microstrip  patches  is  usually 
limited  by  the  reflection  loss  at  the  input  terminal,  the  commonly  used  definition  tor  BW  of 
the  microstrip  patches  is  the  frequency  range  over  which  the  input  VSWR  is  less  than 
certain  specified  value  s.  Usually  s  is  taken  equal  lo  2  but  certain  specifications  choose 
s  =  1.5.  Thus  one  has  to  be  careful  while  comparing  bandwidths  quoted  in  different 
reports. 

A  concept  of  radiation  oandwidth  for  microstrip  auennas  has  been  introduced  in 
[1]  The  radiation  BW  is  defined  as  the  frequency  range  over  which  the  radiated  power  is 
within  3  dB  of  the  incident  power  and  the  radiation  pattern  is  essentially  the  same.  The 
second  part  of  this  definition  can  be  made  more  specitic  by  quantifying  the  allowed 
variations  in  the  radiation  pattern  (for  example  in  beam  direction,  beam  widtn,  or  the 


sidelobe  level).  This  definition  is  more  comprehensive  as  the  radiated  power  is  reduced 


because  of  the  input  reflection  as  well  as  by  other  factors  such  as  changes  in  voltage 
distribution  along  the  edges  oi  the  patch. 

2.  BANDWIDTH  OF  SINGLE  PATCH  ANTENNAS 

Input  impedance  and  the  impedance  bandwidth  of  a  microstrip  patch  antenna  may 
be  computed  using  its  transmission  line  model  or  the  cavity  model.  If  we  represent  the 
resonant  patch  in  terms  of  a  parallel  L-C-G  network,  its  input  admittance  Yjn  can  be  written 

as 

Yjn=G(1  +2jQAf/f0)  (1) 

where  Q  =  cnC/G  and  Af  is  the  deviation  trom  me  center  frequency  fg  Bandwidth  BW 
defined  for  input  VSWR  to  be  better  than  s  is  given  [2]  by 


BW  =  _^L 

QVs  (2) 

For  s=2,  and  substituting  for  Q  in  terms  of  energy  stored  and  power  radiated  we  can  write 


BW- 


V2"hGe 
rtc  Vere  egbe 


(3) 


where  Ge  is  the  edge  conductance,  h  is  the  substrate  thickness,  be  is  the  effective  width 
and  ere  the  effective  dielectric  constant  of  the  equivalent  microstrip  line.  From  (3),  we 
note  that  the  BW  increases  linearly  with  increase  of  h.  Also,  BW  increases  when  t  r  (and 
hence  r.re)  is  reduced.  Recognizing  that  the  edge  conductance  for  a  patch  width  b  is 


given  by  [3], 


byy 


Gp  =  — - — —  for  b/A.n  »l 
6  120  X0  U 


l4) 


and 


"  9o  (b/^°)2  101  b/X°  v<  1 

wo  find  the  fractional  bandwidth  UW  increases  with  the  operating  frequency  Iq 


Typical  variation  of  UW  with  h/X  is  depicted  in  [4)  whereas  variation  with  frequency  is 
shown  in  [2],  Computations  agree  with  the  measured  values.  As  pointed  out  in  [3], 
bandwidths  up  to  15%  can  bo  obtained  at  the  expense  of  an  increase  In  the  antenna 
height  to  about  0.1  wavelengths.  However,  such  thick  substrates  are  not  desirable  or 
practical  becauso  of  tho  following  reasons:  (i)  For  probe-fed  patches,  an  increase  in  the 
substrata  thickness  causes  an  increaso  in  tho  probe  inductance  which,  in  turn,  creates 
input  matching  problems,  (n)  Tor  microstrip  fod  patches,  increased  substrate  thickness 


causes  an  increase  junction  reactance,  which  creates  spurious  radiation  as  we!!  as  the 


input  match  problems,  (iii)  Thick  substrates  aro  not  well  suited  for  microstrip  circuitry  used 
for  signal  distribution  and  phase  shifting  networks  fabricated  on  the  same  substrate,  (iv) 
Thick  substrates  mako  it  mechanically  difficult  to  have  antenna  arrays  contormnl  to  curved 
sudaces  (of  aircraft,  space  crait,  missiles,  etc  ),  (v)  Many  ot  tho  analysis  and  design 
techniques  used  (cavity  model,  etc.)  become  inaccurate  for  thick  substrates. 


Decrease  in  the  OW  of  microstrip  antenna  by  decreasing  tho  substrnto  thickness  h, 
or  decreasing  the  width  b,  or  by  decreasing  the  leng.h  a  (by  increasing  vr)  is  consistent 

with  the  general  relationship  between  tho  Q  and  volume  of  small  antennas  [5) 

A  useful  lormuia  for  bandwidth  of  microstnp  patches  is  given  by  Munson  [0]  as 


nw  =  4f  (  i  ) 

1/32 


(5) 


where  0W  is  the  bandwidth  in  MHz  (for  VSWR  <  2).  t  is  the  operating  Irequency  in  GHz. 
and  t  is  the  thickness  in  inches. 


Non-rectangular  Shaped  Patches 

Bandwidths  available  from  circular  microstrip  patches  [7]  are  of  the  same  order  as 
those  available  from  rectangular  or  square  shaped  antennas.  The  two  other  geometries 
which  have  been  shown  to  yield  wider  bandwidths  are:  circular  ring  [8]  and  square  ring  [9] 
microstrip  antennas  As  pointed  out  in  {8],  TM-j^  mode  of  the  annular  ring  is  a  poor 
radiator  but  TM-|2  mode  of  the  ring  yields  a  bandwidth  about  two  times  larger  than  the 

bandwidth  ol  a  circular  disk.  The  bandwidths  for  the  rectangular  patch  and  for  the 
rectangular  ring  patch  aie  compared  in  {9],  and  in  a  typical  case  (er  =  2.50,  h  =  0.159  cm, 

and  f  =  1 08G  MHz),  the  ring  has  a  BW  value  1.63  times  that  of  the  rectangular  paten. 
Qualitative  reasons  for  improved  bandwidth  are  the  same  in  two  cases.  For  ring 
geometries,  energy  storage  underneath  the  patch  is  smaller  but  the  power  radiated  from 
the  radiating  edges  is  about  the  same. 

3  SINGLE-SUBSTRATE  MULTIPLE-RESONATOR  CONFIGURATIONS 

Use  o!  multiple  resonant  patches  located  in  the  same  plane  (on  the  top  surface  of  a 
single  substrate)  has  been  explored  by  several  investigators  [ 1 0  16).  The  key  idea  in 
these  configurations  is  that  the  stagger  tuned  coupled  resonators  can  yield  wider 
frequency  response  in  the  same  way  as  obtained  in  the  case  of  multi-stage  tuned 
amplifier  circuits. 

In  [10],  two  narrow  conducting  strips  are  placed  close  (-25  to  3  0  times  the 

substrate  thickness)  to  the  the  non-radiating  edges  of  a  rectangular  patch  (fig.  1). 
Bandwidth  of  5%  (VSWR  <  2,  fg  =  350  MHz,  on  1.27  cm  thick  silicone  fiberglass 


substrate)  was  obtained  expenmentally.  Coupling  ol  parasitic  patches  to  the  radiating 
edges  ot  a  rectangular  patch  has  been  described  in  (11).  In  this  case,  the  parasitic 
patches  are  X/4  long  (shorted  at  the  far  ends,  as  shown  in  Fig.  2)  and  therefore  the 
locations  of  the  radiating  apertures  do  not  change.  A  bandwidth  2.12  times  that  of  a 
single  open  patch  (f0  -1.275  MHz,  er  =  2.5,  h  =  3.18  mm,  VSWR  <  2)  has  been  reported. 

The  characteristics  of  this  configuration  have  been  explained  in  terms  of  an  antiphase 
mode  of  a  pair  of  coupled  resonators  and  it  is  shown  that  the  bandwidth  improvement  is 
independent  of  the  coupling  capacitance, 

Multiport  network  analyses  of  multiresonator  configurations  (Fig.  3  a,b  and  c)  with 
additional  patches  gap  coupled  to  the  radiating  edges,  to  the  non-radiating  edges,  or  to 
all  the  four  edges  of  a  nrobe-fed  central  patch  are  reported  in  [1 2]  and  [13].  Coupling 
gaps  are  modeled  by  multiport  capacitive  networks  and  two-dimensional  analysis 
(segmentation  method)  is  used  to  evaluate  input  impedance  as  well  as  the  voltage 
distributions  around  the  edges.  Bandwidth  (VSWft  <  2)  of  the  five-resonator 
configuration  (Fig  3c)  has  been  evaluated  theoretically  and  experimentally  to  be  6.7 
times  that  ot  a  single  patch  on  the  same  substrate  (cr  =  2.55,  h  =  0.318  cm,  and  tg  =  3.16 

GHz,  BW  =  25.8%).  Alternative  multiport  configurations  wherein  the  additional  resonators 
are  not  gap-coupled  but  connected  to  the  central  patch  through  small  sections  ot 
microstrip  lines  (as  shown  in  Fig.  4)  have  been  discussed  in  [I4j.  Again  multiport  network 
modeling  and  planar  analysis  is  employed  and  experimental  verifications  are  reported.  For 
t'ne  5-patch  configuration  (Fig.  4(c)),  a  bandwidth  improvement  of  7.36  times  over  that  for 
a  single  patch  (VSWR  <  2,  rr  =  2.55,  h  =  0.318  cm,  f0  =  3.18  GHz,  BW  =  24%)  have  been 


reported. 


Another  interesting  multi-patch  antenna  configuration  is  described  in  p5].  This 
configuration  consists  of  4  triangular  patches  as  shown  in  Fig.  5.  The  central  patch  A  is 
probe  ted,  the  lower  patch  B  is  gap  coupled,  and  the  other  two  patches  C  and  D  are 
coupled  by  short  sections  of  microstrip  lines.  The  design  was  optimized  by  experimental 
iterations.  A  bandwidth  (VSWR  <  2)  of  11.3%  (tr  =  2.55,  h  =  1.6  mm,  fg  =  3.19  GHz)  which 

is  5.4  times  the  bandwidth  of  a  single  rectangular  patch  antenna  has  been  reported.  Over 
this  bandwidth,  the  3  dB  beamwidth  varies  from  55°  to  72°  in  the  E-plane  and  40°  to  70° 
in  the  H-plane.  The  cross-pol  lc  -el  is  better  than  14  dB  in  the  H-plane  and  12  dB  in  the  E- 
plane.  Analysis  and  design  of  this  type  of  antenna  may  also  be  carried  out  by  the  multiport 
network  modelling  approach  discussed  for  multiresonator  rectangular  patches  in  (12-14], 
Such  antennas  may  find  use  in  applications  where  a  wide  bandwidth  is  required  but 
some  variations  in  radiation  characteristics  over  the  bandwidth  range  may  be  tolerated. 

There  are  two  problems  associated  with  the  use  of  the  multiresonator 
configurations  described  above:  (i)  larger  area  requirement  and  consequent  difficulty  in 
using  these  configurations  as  array  elements,  and  (ii)  variations  cf  the  radiation  pattern 
over  the  impedance  bandwidth  of  the  configuration. 

A  modification  of  the  ruuuiresonaior  patches  (tor  avoiding  the  above-mentioned  two 
problems)  is  described  in  [1b].  This  configuration  (shown  in  Fig.  5)  is  a  modification  cf  the 
non-radiating-edges  coupled  resonator  configuration  (shown  in  Fig.  3b).  Six  narrow 
patches  are  coupled  through  capacitive  gaps.  Their  lengths  are  equal  but  widths  are 
dirferent.  The  configuration  was  designed  by  experimental  iterations.  Measured 
bandwidth  (VSWR  <  2)  equal  to  10  times  that  of  a  single  patch  of  the  same  overall  area  lias 
been  reported  (cr  =  2.2,  h  -  0.8  mm,  fg  -  875  MHz,  BW  -  6%).  Since  the  size  of  this 

multiresonutor  configuration  is  compact,  the  radiation  pattern  is  consistently  good  over 


the  bandwidth.  Measured  radiation  patterns  at  854  MHz  and  897  MHz  support  this 
contention,  information  regarding  variation  of  gam  over  the  bandwidth  is  not  available  in 
(16].  Analysis  and  design  procedures  for  this  type  of  antenna  have  not  been  reported  so 
far.  Of  course,  the  multiport  network  modeling  approach  used  for  3-patch  configuration  of 
Pig.  3b  can  be  extended  to  the  six-patch  configuration  also.  An  alternative  approach  is  to 
model  the  multiple-resonator  patch  as  a  section  of  a  multiple  coupled  microstrip  line,  and 
to  construct  a  coupled-line  model  (similar  to  the  transmission  line  model  for  a  single 
patch;.  This  concept  of  coupled-line  model  as  applicable  to  a  two-resonator  configuration 
is  discussed  in  [1],  Certainly,  multiple-resonator  configuration  with  narrow  patches 
coupled  along  non-radiating  edges,  is  an  antenna  configuration  for  which  further  detailed 
investigations  are  recommended. 

4.  MULTIPLE-RESONATOR  CONFIGURATIONS  WITH  PATCHES  STACKED 

VERTICALLY 

In  this  approach,  two  or  more  layers  of  dielectric  substrates  are  used.  Resonant 
patches  are  located  on  the  top  of  each  of  the  substrate  layers  and  are  stacked  vertically. 
Two-layer  con'  gurations  are  most  common,  although  three  dielectric  layers!  and  three 
patches)  have  also  been  used.  The  two  patches  may  be  identical  (same  size)  or  of  slightly 
ditferent  sizes  (and  hence  resonant  at  slightly  different  frequencies).  When  the  patches 
have  unequal  dimensions,  two  different  arrangements  are  possible.  Smaller  patches  can 
be  located  on  the  top  layer  (Fig.  7a)  or  the  larger  patch  can  be  placed  on  the  top  (Fig.  7b). 
Both  of  these  arrangements  have  been  used.  When  the  smaller  patch  is  on  the  top, 
edges  of  both  the  smaller  and  the  larger  patches  become  the  sources  of  radiation  with  the 
effective  aperture  shifting  from  the  bigger  patch  to  the  smaller  patch  as  the  trequency  of 
operation  is  increased.  On  the  other  hand,  when  the  larger  patch  is  on  the  top  (and  the 
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smaller  patch  below),  the  boundary  of  the  upper  patch  constitutes  the  radiating  aperture. 
The  lower  patch  helps  in  the  broadband  excitation  ot  the  upper  patch,  and  is  termed  as 
the  feeder  patch.  Of  course,  two  patches  can  have  identical  dimensions  in  which  case 
the  distinction  between  the  feeder  and  the  radiator  patches  disappears  and  the  two 
functions  are  merged  together.  In  most  of  these  two-patch  configurations,  the  lower 
patch  is  fed  via  a  probe  or  a  microstrip  line.  Two  different  ways  of  excitations  ot  the  upper 
patch  are  available.  We  can  use  a  vertical  probe  passing  through  a  hole  in  the  lower  patch 
and  making  contact  to  the  upper  patch.  Alternatively  capacitive  coupling  between  the 
patches  is  used  to  excite  the  upper  patch. 

Double  layer  two-patch  configurations  were  initially  used  for  dual  frequency 


microstrip  antennas  [10,17].  Wide  bandwidth  application  was  reported  by  Hal!  et  a!  in  [IS] 


wherein  two-layer  and  three-layer  configurations  are  described.  In  these  designs,  the 
lower  substrate  layer  is  alumina  (sr  =  9.8,  h  =  0.625  mm)  where  the  feed  line  and  the 
associated  circuitry  is  fabricated.  The  second  and  third  layers  are  polyguide  (er  =  2.32,  h 

=  1.59  mm)  where  the  resonant  radiating  patches  are  fabricated  (see  Fig.  8  a  and  b). 
Bandwidths  ot  13%  and  18%  (return  loss  >  10  dB,  fg  ~  10.6  GHz)  were  measured  for  two- 

layer  and  three-layer  configurations  respectively.  For  comparison,  bandwidth  of  a  single 
layer  patch  on  alumina  substrate  is  only  1 .1%,  while  that  on  polyguide  is  about  6.6%.  This 
technique  of  capacitively-coupled  stick-on  antennas  is  well  suited  for  certain  applications 
where  the  antenna  may  be  regarded  as  an  additional  component  in  the 
transmitter/receiver  circuits. 

There  have  been  several  reports  of  two  layer  antennas  with  an  air-gap  in  between 
the  two  layers.  In  the  configuration  described  in  [19],  the  lower  patch  has  a  smaller  size 
and  acts  as  a  feeder  resonator.  Out  of  the  various  configurations  presented,  S-band 
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circular  disc  configuration  has  a  bandwidth  of  15%  (VSWR  <  2.0,  air-gap  thickness  not 
specified).  For  the  configuration  reported  in  [20],  a  low  er  dielectnc  foam  material  is  used 

to  separate  the  two  patches.  For  a  spacing  of  0.572  cm,  a  17.3%  BW  (VSWR  <  1 .92, 
3.85-4.58  GHz)  is  reported.  The  lower  patch  is  smaller  in  size  and  is  fed  by  a  probe  or  a 
microstrip  hne.  The  element  has  been  used  in  4x4  array  of  circularly  polarized  patches 
(two  feed  points  on  each  patch).  An  investigation  of  the  two  layer  antenna  performance 
as  a  function  of  the  air-gap  spacing  is  discussed  in  [21  ] 

An  experimental  study  of  vertically  stacked  triangular  patches  is  reported  in  [22], 
Both  dual  frequency  and  broadband  operation  are  described.  For  two  identical  equilateral 
tnangular  patches  fabricated  on  dielectric  layers  (with  er  =  2.55  and  thickness  1.6  mm) 
separated  by  a  foam  thickness  of  5.0  mm,  a  bandwidth  of  17.46%  (VSWR  <  2,  fg  =  3  61 
GHz)  has  been  reported.  For  comparison,  the  bandwidth  of  the  same  two-patch 
configuration  with  zero  foam  thickness  is  only  2.5%. 

All  of  the  vertically-stacked  multiple  resonator  configurations  discussed  above  have 
been  designed  by  experimental  iterations.  However,  analyses  of  such  configurations 
based  on  equivalent  transmission  line  model  [23]  as  well  as  based  on  full-wave  spectral 
domain  approach  [24,25]  have  been  reported  recently.  In  [23],  the  transmission  line 
model  analysis  of  a  two-layer  rectangular  patch  is  described.  The  upper  patch  is  smaller, 
and  the  region  common  to  both  the  patches  is  modeled  as  a  coupled  line  section. 
Conceptually,  the  approach  is  similar  to  that  described  in  [26]  for  a  suspended  patch 
antenna  excited  by  an  electromagnetically  coupled  inverted  microstrip  feed.  For  the 
configuration  desenbed  in  [23],  a  BW  of  13.5%  (f0  =  13  GHz,  VSWR  <  2,  lower  and  upper 
:  jbstrates  have  h  =  0,8  mm  and  cr  =  2.55)  has  been  repoded.  For  comparison  the 
bandwidth  for  a  single  patch  is  only  4.46%. 


Hankel  transform  domain  analysis  of  two-layer  cncular  patch  contiguration  reported 
in  [24,25]  is  an  extension  of  the  analysis  of  a  single  layer  open  circular  microstrip  raoiatmg 
structure  discussed  in  [27],  Complex  resonant  frequencies  of  the  unloaded  two-patch 
configuration  are  found  by  solving  the  equivalent  transverse  transmission  line  circuits  for 
TE  and  TM  components  in  the  Hankel  transform  plane.  Basis  functions  chosen  for 
current  distribution  take  into  account  the  edge  condition.  Green's  functions  in  spectral 
(Hankei  transform)  domain  are  in  algebraic  form.  Coefficients  of  the  immittance  matrix  are 
obtained  via  complex  integration  in  trie  spectral  domain.  Equating  coefficients 
determinant  to  zero  yields  complex  resonance  frequencies.  Input  VSWR  characteristic  is 
obtained  by  using  natural  complex  frequencies  computed  and  an  equivalent  circuit 
model.  Radiation  patterns  are  computed  without  inverting  the  Hankel  transform  current 
distribution.  Numerical  resonant  frequency  and  VSWR  variation  are  shown  to  be  in  good 
agreement  with  the  experimental  data  (er  =  2.55,  f  -  2.4-2. 8  GHz,  dielectric  layers  1 .6  mm 

thick,  air  spacing  between  layers  10.0  mm,  lower  patch  radius  20.8  mm,  upper  patch 
radius  21 .0  mm). 

Vertically  stacked  multi-resonator  configurations  are  the  most  widely  used 
broadband  microstrip  antenna  elements.  But,  convenient  practical  design  and  analysis 
procedures  are  not  available  so  far.  Analytical  treatments  based  on  the  extension  ot  cavity 
and  multiport  network  models  are  likely  to  be  developed  in  the  near  future. 

5.  IMPEDANCE  MATCHING  NETWORKS  FOR  BROAD-BAND  MICROSTRIP 

ANTENNAS 

Approaches  for  widebanding  of  microstrip  antennas  discussed  so  far  are  based  on 
either  the  concept  of  decreasing  the  resonator  Q-factor,  or  on  the  concept  of  coupled 
multiresonator  circuits.  An  entirely  different  approach  makes  use  of  wideband  impedance 
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matching  networks  to  reduce  the  reflection  loss  at  the  input  of  the  resonant  patch 
antenna.  It  is  an  obviously  natural  solution  from  microwave  circuit  designers'  point  of  view. 
However,  a  very  small  number  [28-32]  of  microstrip  antenna  designers  have  looked  upon 
the  broadbanding  problem  from  this  angle. 

As  pointed  out  in  [28,29,32],  the  maximum  bandwidth  obtainable  by  using  an 
impedance  matching  network  can  be  calculated  by  using  Fano's  broadband  matching 
theory  [33].  Treating  the  microstrip  patch  antenna  as  a  parallel  resonant  circuit  (as 
discussed  in  Sec.  2),  it  can  be  shown  the  input  impedance  bandwidth  BW  is  given  by 


BW  =  —  J  {(- 
G  V  GZ,- 


1  )(s-l  )/s } 


where  G  is  the  shunt  conductance  in  the  parallel  resonant  network  and  Zq  is  the 
characteristic  impedance  of  the  input  transmission  line.  When  Zq  =1/G,  the  relation  (6) 
reduces  to  (2)  of  Section  i .  However,  one  can  find  a  value  of  Zq  that  will  maximize  the  BW 
in  (6)  and  this  is  given  by 


Z°lo^,mum=G(F7T7i) 


For  this  optimum  value  of  Zq  the  optimum  BW  is  given  by 
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For  ttie  usual  value  of  s=2,  BW0p(  given  by  (8)  is  about  1 .2  times  the  bandwidth  given  by 
(2).  It  implies  that  if  we  do  not  insist  on  a  perfect  match  at  the  center  frequency,  the 


frequency  range  over  which  the  input  VSWR  is  less  than  2  can  be  increased  by  about 
20%.  This  is  consistent  with  Fano's  theory  [33]  according  to  which  maximum  BW  over 
which  a  parallel  resonant  circuit  can  be  matched  for  a  VSWR  <  sr  is  given  by  [28], 


BWmax=— — - - - 

max  Q  £n{(s+1)/(s-1  )}  (9) 

Comparing  the  values  given  by  (9)  with  those  given  by  (2),  it  can  be  shown  that  a 
bandwidth  improvement  by  a  factor  of  4.04  can  be  obtained  (for  s  =  2)  by  designing  a 
perfect  impedance  matching  network. 

Experimental  results  reported  in  [28]  point  out  a  bandwidth  improvement  by  a  factor 
of  3.9  (3  GHz,  BW  =  10%,  $  =  1.9)  by  using  a  shielded  microstrip  matching  network 
consisting  of  two  resonator  elements  etched  on  the  same  substrate  as  the  patch. 
Configuration  reported  in  [29],  on  the  other  hand,  makes  use  of  the  feed  probe 
inductance  (occurnng  in  a  thick  substrate  coaxial-fed  microstrip  antenna  element)  as  a  part 
ot  impedance  matching  network.  A  series  capacitance  located  at  the  top  end  of  the  probe 
(Kapton  layer  shown  in  Fig.  9)  forms  a  series  resonant  circuit  which  is  used  tor  broadband 
matching  of  the  parallel  resonant  circuit  constituted  by  the  resonant  patch  antenna. 
Experimental  results  show  a  marked  improvement  in  the  return  loss  (~ii.3  GHz,  er1  = 
1.05,  hi  =1.6  mm,  er2  =  3.3  mm,  h2  =  0.15  mrn)  but  bandwidth  improvement  data  is  not 
reported  in  [29], 

Another  wideband  circular  disc  antenna  making  use  of  a  series  resonant  circuit 
(teed  probe  inductance  being  a  part  of  the  series  resonance)  is  described  in  [30].  The 
matching  network  (details  not  given)  was  optimized  using  CAD  techniques  and  interlaced 
directly  with  a  triplate  array  teed  network.  For  a  6  mm  radius  patch  (er  =  2.2,  h  =  3.2  mm)  a 
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bandwidth  ot  35%  of  VSWR  <  1 .5  has  been  reported.  Radiation  patterns  over  the  whole 
band  (7  to  10.5  GHz)  shows  little  variation.  Gain  is  about  6  dBi. 

A  forthcoming  paper  by  Pues  and  Van  de  Capelle  [32]  describes  in  detail  the 
transmission-line  matching  network  design  for  broadband  microstrip  antennas 

A  single  stub  impedance  matching  technique  for  increasing  the  bandwidth  of  an 
electromagnetically  coupled  microstrip  patch  (Fig.  10)  is  discussed  in  [31].  In  this  case, 
the  input  impedance  of  the  electromagnetically  coupled  patch  was  measured  by  a 
network  analyzer  and  component  values  for  a  network  model  (a  capacitance  in  series  with 
a  parallel  RLC)  were  obtained  by  a  computer-aided  model  fitting  procedure.  A  stub  circuit 
was  then  designed  to  match  to  this  load.  A  bandwidth  of  13%  (VSWR  <  2,  3.375  GHz  to 
3.855  GHz,  eri  =  er2  =  2.2,  and  di  =  dg  =  0.158  cm)  was  thus  obtained. 

As  evident  from  the  above  discussions,  broadband  impedance  matching  networks 
provide  a  really  effective  means  of  increasing  the  usable  bandwidth  of  microstrip 
antennas.  A  combination  of  this  technique  with  the  other  approaches  discussed  earlier 
woulq  perhaps  lead  to  the  optimum  configuration  for  wide  bandwidth  microstrip  patches. 
Electromagnetically  Coupled  Microstrio  Patches 

The  example  of  broadband  patch  discussed  above  (described  in  [31])  makes  use  of 
an  electromagnetic  coupled  patch  (Fig.  10)  consisting  of  two  dielectric  layers.  The 
feeding  microstrip  line  is  printed  on  the  lower  layer  and  the  patch  is  on  the  top  of  the 
upper  dielectric  layer.  Such  a  coupling  arrangement  is  frequently  used  in  printed 
microstrip  antennas  [18,  26,  35-37],  Electromagnetically  coupled  patches  are  preferred 
as  the  feeding  microstrip  line  network  and  other  circuitry  are  located  on  a  thinner  (and 
possibly  higher  er  if  desired)  substrate  and  the  radiating  resonant  patch  is  fabricated  on 
the  thicker  and  low  er  layer. 
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Electromagnetically  coupled  patches  are  very  convenient  from  the  broadbanding 
point  ot  view  also.  Impedance  matching  networks  may  be  located  on  the  lower  layo' 
underneath  the  patch  without  requiring  additional  substrate  area.  Also  the  capacitive 
coupling  between  the  patch  and  the  microstnp  line  constitutes  a  senes  capacitance  which 
can  be  utilized  as  a  component  of  the  impedance  matching  network. 

Some  of  these  considerations  are  implicit  in  the  bandwidth  enhancement  method 
for  electromagnetically  coupled  strip  dipoles  discussed  in  [36].  In  this  paper,  it  is 
demonstrated  that  'f  parasitic  metallic  strips  are  incorporated  in  the  structure  either  co- 
planar  and  parallel  to  the  embedded  inicrostrip  line  open-end,  or  between  the 
transmission  line  and  microstrip  dipole,  then  substantial  bandwidth  enhancement  can  be 

ObiBiflGu . 

6.  BANDWIDTH  CONSIDERATIONS  FOR  CIRCULARLY  POLARIZED  PATCHES 

Two  different  design  philosophies  are  used  for  circularly  polarized  (CP)  microstrip 
patches.  One  of  these  makes  use  of  the  two  spatially  orthogonal  modes  of  a  square  or  a 
circular  patch.  These  two  modes  (each  of  them  giving  rise  to  a  linearly  polarized  radiation 
when  excited  independently)  are  fed  through  two  orthogonally  located  feeds  by  two 
signals  which  are  equal  in  magnitude  but  are  in  phase  quadrature.  The  bandwidth  of 
these  two-feed  circularly  polarized  antennas  is  the  same  as  that  for  the  corresponding 
linearly  polarized  antennas,  and  all  the  broadbanding  techniques  discussed  above  are 
applicable  to  this  class  of  circularly  polarized  antennas  also.  Reference  [38]  is  an  example 
where  two  vertically  stacked  circular  resonator  configuration  is  used  to  obtain  a  wide-band 
circularly  polarized  patch.  A  branch  line  hybrid  is  used  to  split  the  incoming  signal  for  two 
orthogonally  located  feeds.  A  bandwidth  of  13%  (VSWR  <  1.5)  is  obtained  when  the  two 
patches  (fabricated  on  r.r  =  2.55,  1 .6  mm  thick  dielectric  layers)  are  separated  by  air  gap  of 
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0.08  -  C.09  Xq .  The  radius  of  the  upper  patch  is  1.0  to  1.1  times  that  of  lower  patch 
(frequency  -  2.4  to  2.8  GHz). 

The  second  design  approach  used  for  circularly  polarized  patches  does  not  require 
a  signal  splitting  network  and  makes  use  of  only  a  single  feed  patch.  There  are  several 
configurations  [39-42]  based  on  this  approach.  All  of  these  configurations  make  use  ot 
two  spatially  orthogonal  modes  of  the  patch.  The  resonance  frequencies  of  these  two 
modes  are  made  slightly  different  from  each  other.  The  patch  is  led  at  a  frequency  in 
between  these  two  resonances  and  the  feed  location  is  selected  such  that  the  two 
modes  are  excited  equally.  Since  the  resonance  frequency  ot  one  of  the  modes  is  lower 
than  the  signal  frequency  and  that  of  the  other  mode  is  higher  than  the  signal  frequency, 
a  proper  selection  of  frequencies  causes  the  two  mooes  to  be  exciieti  9Cr  uui  of  phase. 
Thus  the  conditions  for  a  CP  radiation  are  met. 

The  bandwidth  of  the  single  feed  CP  patches  is  limited  not  by  the  input  VSWR  (as  is 
the  case  for  linearly  polarized  patches),  but  by  the  degradation  of  the  axial  ratio  as  one 
moves  away  from  the  center  design  frequency.  Howaver,  an  analysis  based  on  the  cavity 
model  [43]  shows  that  the  required  separation  between  the  two  resonance  frequencies 
and  hence  the  axial  ratio  bandwidth  is  inversely  proportional  to  the  Q-factor  of  the  patch. 

A  comparison  of  the  axial-ratio  bandwidths  for  vanous  CP  patches,  given  in  [40], 
shows  that  the  square  ring  configuration  has  a  wider  bandwidth  (5.2%,  axial  ratio  <  6  oB, 
f  =  3.0  GHz,  er  -  2.50,  and  h  *»  0.159  cm)  as  compared  to  crossed  strip,  almost  square 

patch,  corner  chopped  square  patch,  and  a  square  patch  with  a  diagonal  slot 
configurations.  Another  wideband  single  feed  OP  configuration  is  annular  ring  with  an  ear 
configuration  discussed  in  [42].  In  this  case,  TM12  mode  is  used  and  higher  bandwidth  is 
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obtained  tor  small  ratios  of  the  outer  to  inner  radii  of  the  nng.  A  6%  bandwidth  (axial  ratio  < 

6  dB,  er  =  2.52,  h  =  0.159  cm,  f  -  4,5  GHz)  is  reported. 

As  discussed  in  Section  2,  rectangular  ring  and  circular  ring  patches  provide  the 
maximum  impedance  bandwidth  among  the  single-layer,  single-resonator  patch 
configurations.  "Thus  it  is  not  surprising  to  note  that  these  two  configurations  also  provide 
the  maximum  axial  ratio  bandwidth  when  used  as  single-teed  CP  configurations. 

7.  BROAD  BAND  MICROSTRIP  ARRAY  CONFIGURATIONS 

As  the  individual  microstrip  patches  exhibit  antenna  gains  of  the  order  of  4-6  dB 
only,  they  are  frequently  used  in  linear  and  two-dimensional  array  configurations.  As  a 
general  rule,  arrays  of  broadband  microstrip  patches  yield  wide  bandwidth  This  is 
illustrated  in  [38j  where  4x4  array  of  broadband  patches  is  shown  to  yield  8.5%  bandwidth 
(VSWR  <  1.5,  gain  >  18  8  dBi,  efficiency  >  62%,  SLL  <  -22  dB)  at  2.6  GHz  using  the 
broadband  CF  patches  discussed  earlier. 

Specially  designed  broadband  configurations  for  linear  arrays  have  also  been 
reported  [44,45],  The  configuration  discussed  in  [44]  is  a  senes  ted  linear  array  ot 
rectangular  elements.  The  individual  elements  of  this  array  are  two-port  rectangular 
patches.  For  two-port  patches,  the  VSWR  bandwidth  depends  not  only  on  the  patch  but 
also  on  the  power  transmitted  to  the  output  port.  Thus  much  wider  impedance 
bandwidths  are  available  in  series-fed  array.  The  design  reported  in  [44]  has  a  bandwidth 
or  40%  (VSWR  <  2)  and  a  gain  of  ~10  dB  at  x-band.  The  configuration  employs  a  double 
section  non-uniform  antenna  (Fig.  11)  having  a  length  equal  to  about  5.0  Xmax  and 
consisting  ot  two  rows,  spaced  about  1 .5  A.max,  each  of  seven  radiating  patches. 

The  application  of  the  log-periodic  technique  to  the  series-fed  electromagnetically 
coupled  overlaid-patch  array  has  been  discussed  in  [45],  A  k-f}  analysis  ot  microstrip 


arrays  is  reported  and  indicates  that  the  microstnp  patch  is  not  an  optimum  element  lor 
log-penodic  arrays  and  that  elements  connected  directly  to  the  teed  line  result  in  the 
arrays  having  a  limited  bandwidth.  The  addition  of  the  series  capacitance  to  the  patcti 
equivalent  circuit,  implemented  by  electromagnetic  coupling,  allows  an  optimum  to  be 
approached  Log-periodic  overlaid  patch  array  design  (Fig.  12)  has  been  reported  to 
have  a  4  i  bandwidth  with  an  input  return  loss  of  8  dB,  a  gain  of  8  dBi  and  a  30°  backfire 
beam  whose  beamwidth  varies  from  63  to  32  degrees  across  the  band.  In  spite  ot 
limitations  on  the  ultimate  BW  obtainable,  this  design  certainly  extends  the  useful 
application  areas  of  microstrip  antennas. 

8.  CONCLUDING  REMARKS 

Various  techniques  used  for  increasing  the  bandwidth  of  microstrip  patch  antennas 
have  been  reviewed,  it  is  pointed  out  that  the  consideration  of  the  impedance  bandwidth 
alone  is  not  sufficient  for  broadband  microstrip  antennas,  and  a  concept  of  radiation 
bandwidth  is  introduced. 

It  has  been  discussed  (in  Sec.  2)  that  the  Q-factors  of  the  resonant  patches  can  be 
decreased  and  hence  their  bandwidths  can  be  increased  by  increasing  the  height  of  the 
dielectric  substrates.  The  problems  associated  with  thick  substrates  are:  increased 
reactance  ot  the  feed  junction,  and  the  radiation  from  the  microstrip  circuits  included  on 
the  same  substrate.  These  difficulties  may  be  overcome  by  using  electromagnetically 
coupled  patches  with  feed  network  and  associated  circuitry  on  the  lower  (and  thinner) 
substrate.  There  is  a  need  to  develop  design  procedures  tor  electromagnetically  coupled 
patches  based  on  the  extension  of  cavity  mode!  oi  multiport  network  model. 

Among  the  three  other  approaches  used  for  obtaining  wider  bandwidth,  more 
attention  needs  to  be  directed  to  the  application  of  wideband  impedance  matching 
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netwoik  approach  discussed  in  Sec.  5.  The  only  difficulty  with  this  method  is  'he  need  tor 
larger  substrate  area  required  tor  incorporating  the  matching  network  Again,  the  use  of 
electromagneticaily  coupled  patches  snould  prove  very  useful  since  the  matching 
network  can  be  incorporated  at  the  lower  level,  possibly  underneath  the  palch  itself. 
Further  research  efforts  need  to  be  pointed  in  this  direction. 

Between  the  two  different  methods  ot  constructing  muitiresonalor  coupled  patch 
configurations,  the  two-layer  configuration  with  vertically-stacked  patches  requires  smaller 
area  and  does  not  suffer  much  from  pattern  degradation  with  frequency  This 
configuration  is  used  frequently  but  most  of  the  designs  are  carried  out  by  experimental 
iterations.  Convenient  models  and  analysis  procedures  need  to  be  developed  so  that 
the  design  and  optimization  ot  these  configurations  can  be  cam  i  out  more 
systematically  Among  the  various  configurations  with  multiple  patches  on  the  same  layer, 
the  configuration  with  the  six  narrow  patches  coupled  along  the  non-radiating  edges  [15] 
has  yielded  the  most  promising  results.  Further  research  leading  to  analysis  and  design 
ot  this  type  of  geometry  is  desirable.  A  multiple  coupled  rmcrostrip  line  model,  proposed 
recently  tor  two  patches  coupled  along  the  non-radiating  edges  [i],  appears  to  be  a 
reasonable  approach  tor  this  purpose. 

The  search  tor  an  "ideal"  wideband  printed  microstrip  antenna  is  still  on.  Perhaps  a 
combination  ot  various  approaches  discussed  in  this  review  would  lead  to  an  optimum 
broadband  configuration.  We  can  look  forward  to  continued  research  in  this  area. 
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^  TEFLON  FIBERGLASS  (3.2mm  THICK) 


Fig.  1  A  rectangular  patch  antenna  with  parasitically  coupled  strips  (10). 
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Fin  2  A  rectangle  wltn  wo  V4  snorted  pa.cnes  coupled  to  the  non-radiating  edges 
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Fig.  3  (a)  Radiatmg-edges  gap-coupled  microstrip  antenna  [12]; 

(b)  Non-radiating  edges  gap-coupled  microstrip  antennas  (13];  and 

(c)  Four-edges  gap-coupled  microstrip  antenna  [13]. 
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Fig.  4  (a)  Radiating-edges  directly-coupled  microstrip  antenna  [14]; 

(b)  Non-radiating  edges  directly-coupled  microstrip  antenna  [14]  and 

(c)  Four-edges  directly-coupled  microstrip  antenna  [14], 
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Fig.  7 
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(a)  Two-patch  two-layer  antenna  with  smaller  patch  located  on  the  top  layer. 

(b)  Two-patch  two-layer  antenna  with  larger  patch  located  on  the  top  layer. 
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Fig.  8  (a)  and  (b)  Examples  ot  wideband  microstnp  anlennas  lor  circuit  integration  [18] 


antenna  patch 


Fig.  9  A  microstrip  anienna  with  a  series  L-C  network  used  for  widebana  impedance 
matching  [29]. 


Fig.  10  A  single  stub  impedance  matching  technique  for  increasing  the  bandwidth  of  an 
electromagnetically  coupled  microstrip  patch  [31]. 
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rig.  1 1  Double-section  wideband  series-linear  array  of  microstrip  patches  [44], 


Fig.  12  A  log-periodic  overlaid  array  of  electromagnetically  coupled  microstrip  patches 


CHARACTERISTICS  OF  ELECTROMAGNETICALLY  COUPLED  RECTANGULAR  PATCH 


ANTENNAS  WITH  OVERLAYING  PARASITIC  ELEMENTS 


by 

R.  Q.  Lee 

NASA  Lewis  Research  C  nter 
Cleveland,  OH  44135 


K.  F.  Lee 

Department  of  Electrical  Engineering 
The  University  of  Toledo 
Toledo,  OH  43606 


ABSTRACT 


The  ciidi  dc  Lei  i  s  i  i  cs  of  el  ecLt  omaynei  i  cal  iy  coupled  patch  (EMCP) 
antennas  consisting  of  a  driven  rectangular  element  in  the  bottom 
and  identical  parasitic  elements  on  the  top  have  been  investigated 
experimentally.  The  antenna  was  excited  at  the  resonant  frequency 
of  the  TMqi  mode  with  teflon  spacers  to  vary  the  spacings  between 
the  antenna  layers.  Experimental  results  indicate  that,  for  the 
two-layer  EMCP  antenna,  the  characteristics  can  be  separatee  into 
three  regions  depending  on  the  spacings  between  the  driven  and 
parasitic  layers.  Significant  improvement  in  bandwidth  is  possible 
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only  In  region  1,  and  high  gain  with  about  one  percent  bandwidth  Is 
achievable  In  region  3.  Region  2  generally  has  poor  radiation 
characteristics  and  Is  therefore  not  suitable  for  antenna 
applications.  For  the  three-layer  EMC?  antenna  operated  In  the  high 
gain  region  3,  experimental  results  show  that  the  overlaying 
parasitic  elements  can  enhance  the  gain  (db)  and  reduce  the 
beamwldth  of  the  antenna  by  a  factor  of  two  or  better. 

INTRODUCTION 

Microstrip  antennas  are  known  to  have  narrow  bandwidth  and  low 
gain.  In  recent  years,  research  on  bandwidth  improvement  utilizing 
paiasitic  element  has  been  extensive  and  significant  results  have 
been  made.  By  coupling  coplanar  parasitic  elements  to  a  fed  patch 
antenna,  a  bandwidth  five  times  as  large  as  that  of  a  single  patch 
has  been  reported  [1].  The  problem  of  gain  enhancement  has  received 
disproportionately  little  attention  in  the  literature.  A  method  to 
Increase  the  gain  of  printed  circuit  antennas  using  superstates  has 
been  proposed  by  N.  G.  Alexopoulos  et  al  [2].  This  method  requlies 
the  use  of  a  suoerstrate  with  eithei  very  large  dielectric  constant 
e  or  permeability  p.  The  practicality  of  this  method  depends  on  the 
availability  of  such  materials.  In  a  recent  article,  it  was 
demonstrated  that  by  placing  closely  spaced  parasitic  elements  in 
the  same  plane  of  a  driven  paten,  substantial  enhancement  of  the 
gain  was  obtained  [3].  It  appears  that,  as  in  the  case  of 
bandwidth,  parasitic  elements  play  an  important  role  in  enhancing 
the  gain  of  a  microstrip  patch. 
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This  paper  reports  an  experimental  investigation  of  a  two-layer  and 
a  three-layer  EMCP  antennas.  Although  the  two-layer  geometry  has 
been  studied  by  several  investigators  [4,  53,  the  emphasis  of 
previous  studies  was  on  bandwidth  improve. .lent ,  and  the  spacing,  s, 
between  antenna  layers  seldom  exceeded  0.15  X  (X  =  dielectric 
wavelength) . 

In  a  recent  study  C 6 ]  of  a  two-layer  EMCP  antenna  with  air  spacing, 
It  has  been  established  that,  expending  on  the  spacing,  s,  the 
antenna  characteristics  can  be  divided  into  three  regions: 
significant  improvement  in  bandwidth  is  possible  only  in  region  1, 
and  high  gain  is  achievable  in  region  3.  Region  2  generally  has 
poor  radiation  characteristics  and  is  therefore  not  suitable  for 
antenna  applications.  The  high  gain  region  begins  at  approximately 
0.3  X  (.61cm)  and  is  accompanied  by  a  narrow  bandwidth  of  1-2%.  It 
is  therefore  reminiscent  of  a  Vagi  antenna.  In  this  paper,  results 
for  the  characteristics  of  a  two-layer  and  a  three-layer  EMCP 
antennas  with  teflon  spacers  separating  the  antenna  layers  are 
presented.  The  three-layer  EMCP  antenna  with  two  overlaying 
parasitic  directors  was  designed  to  operate  in  the  high  gain  region. 

EX PERIHEMAL  RESULTS 

The  microstrip  patch  used  in  our  experiment  is  shown  in  fig.  la. 

The  patch  with  dimensions  a=  1 . 5  cm  and  b-1.0  cm  is  fabricated  on 
Cuflon  substrate  with  £,-=2. 17  and  thickness  0.025  cm  (10  mil).  The 


two-layer  EMCP  antenna  consisting  of  an  overlaying  parasitic  element 
and  a  driven  patch  is  shown  in  Fig.  1b.  The  driven  patch  was 
excited  with  &  coaxial  probe  at  the  resonant  frequency  of  the  TMq^ 
mode.  In  the  experiment,  the  spacing  betweer  antenna  layers  was 
incremented  in  multiples  of  0.381  mn  (15  mil)  using  teflon  spacers. 
For  each  spacing,  s,  the  E  and  H- plane  patterns  and  the  swept 
frequency  response  were  measured.  From,  these,  the  resonant 
frequency  as  well  as  the  antenna  gain  can  be  estimated.  The 
estimated  gains  computed  from  the  formula  [7] 

26000 

Gain  (db)  ~  10  log  <  * 

HPE°  HPh° 


are  found  to  agree  to  within  ±  0  6  db  of  the  measured  values.  The 
characteristics  for  the  case  of  10  mil  substrates  for  both  driven 
and  parasitic  patches  are  summarized  in  Table-1.  As  in  the  case 
with  air  separation,  the  antenna  characteristics  can  be  divided  into 
three  regions:  region  1  (0  <  s  <  0.114  cm)  exhibits  a  double-humped 
response  charactei i s i tes  with  an  estimated  gain  of  about  7  db; 
region  2  (0.305  <  s  <^0.381  cm)  has  "aonormal"  radiation  patterns 
while  region  3  (s  >  0.61  cm)  exhibits  "normal"  pattern  shapes  arid 
has  an  estimated  gain  of  about  11  db.  Representative  patterns  for 
each  region  are  shown  in  Fig.  2. 


To  evaluate  the  concept  of  gain  enhancement  using  parasitic 
director*,  a  three-layer  EMCP  antenna  with  two  identical  rectang  iar 
parasitic  patches  was  tested.  The  EMCP  antenna  was  designed  with 
proper  spacings,  s,  for  operation  in  the  high  gain  region.  For 
comparison  purpose,  antenna  gain,  radiation  patterns  and  swept 
frequency  response  of  a  single  micrcstrip  patch  and  a  two-layer  EMCP 
antenna  were  also  measured.  The  experimental  results  are  given  in 
Table  2.  Based  on  these  results,  the  following  comments  can  be  made 

a)  The  presence  of  parasitic  directors  reduces  the  3  db 
beamwidth  of  both  the  E  and  FI-plane  patterns.  The  3  db 
beamwidth  of  the  E-plane  is  reduced  from  103°  to  45°  with 
one  parasitic  director  and  to  30°  with  two  parasitic 
directors,  while  the  3  db  beamwidth  of  the  H-plane  is 
reduced  from  70°  to  50°  to  35°  respectively.  The  E  and 
H-plane  patterns  are  shown  in  Figures  3  and  4. 

b)  As  indicated  in  the  swept  frequency  measurement  in  Fig.  5, 
the  parasitic  directors  reduce  the  resonant  frequency  from 
10.12  GHz  for  the  single  patch  to  10.06  GHz  for  the 

two- layer  and  to  10.02  GHz  for  the  three- layer  EMCP 
antennas . 

c)  The  spacings  for  operating  in  the  high  gain  region  are 
0.35  X  for  the  two-layer;  0.35  X  and  0.47  X  for  the 
three-layer  antennas. 


d)  Antenna  gains  were  measured  using  the  image  method  under 
mismatched  conditions  [8].  The  measured  gain  for  the 
single  patch  is  4.7  db  compared  to  8.4  db  for  the 
two-layer  and  10.6  db  for  the  three-layer  EMCP  antennas. 
It  appears  that  further  increase  in  gain  is  possible  with 
additional  directors. 


CONCLUSION 


In  conclusion,  experimental  results  for  a  two-layer  and  a 
three-layer  EMCP  antennas  have  been  presented.  In  addition  to  the 
large  bandwidth  region  at  s  <  .114  Cm,  a  high  gain  region  is  found 
0.1  s  >  0.3  X.  For  the  three-layer  EMCP  antenna  operated  in  the  high 
gain  region,  the  parasitic  directors  enhance  the  gain  of  a 
rectangular  patch  antenna  from  4.7  db  to  10.6  db,  and  reduce  the  3 
db  bearnwidth  from  103°  to  30°  for  the  E-plane  and  from  70°  to  35° 
for  the  H-plane.  The  three-layer  EMCP  antenna  exhibits  similar 
characteristics  as  a  Yagi  antenna  with  over  10  db  gain  and  with 
about  1-27.  bandwidth. 
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Table  1 . 

Characteristics  of 

a  rectangular 

EMCP  antenna  with 

teflon 

separation 

:  a=l . 5 

cm,  b=l .0cm,  er 

=2.17,  t] =  1 2=0 . 0254cm. 

Spacing 

foi 

Pattern 

3  dB 

Estimated 

Bandwi dth 

s 

shape 

beamwi dth 

gai  n 

(cm) 

(GHz) 

(deg . ) 

(dB) 

r/°) 

0.000 

10.3 

normal 

88  x  65 

6.6 

7.1 

0.038 

10.8 

normal 

78  x  65 

7.1 

3.0 

0.076 

10.5 

normal 

82  x  65 

6.9 

1  .7 

0.114 

10.5 

normal 

78  x  65 

7.1 

1.5 

0.305 

10.3 

abnormal 

_ 

_ 

0.343 

10.2 

abnormal 

- 

- 

- 

0.381 

10.2 

abnormal 

- 

- 

- 

0.610 

10.1 

normal 

40  x  50 

11.1 

1  .0 

r\  r  a  n 

V  .  UHO 

i  r.  i 

IV.! 

norma  1 

40  x  55 

10.7 

1  .0 

0.762 

10.0 

normal 

40  x  50 

11.1 

1  .0 

Single 

10.2 

normal 

110  x  70 

5.3 

2.3 

patch 

‘Range  of  frequencies  for  which  the 
greater  than  one-half  of  its  value 


real  part  of  the  impedance  is 
at  resonance. 
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Table  2.  Characteristics  of  a  rectangular  electro¬ 
magnet!  cal  ly  coupled  patch  antenna  with 
a=l . 5  crn,  b=l  .0  cm 


Antenna 

Single 

Types 

Patch 

Two- layer 

Three- 1 ayer 

3  do  beamwidth 

(E-plane) 

103° 

45° 

30° 

3  db  beamwidth 

(H-plane) 

O 

c 

50° 

35° 

Gain  (db) 

4.7 

8.4 

10.6 

Ereq.  (GHz) 

10.12 

10.06 

10.02 

Element 
Spac i ng 


0.35X 


0.35X 

0.47X 


Fig.  1  Geometry  of  <a)  rectangular  microstrii 
patch  and  (b)  two-layer  EMCP  antenna. 
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Patterns  of  a  rectangular  EMCP  antenna  with 
teflon  separation:  a=1.5  cm,  b=1.0cm,  Gr«2.17, 
t]=t2=0.0254  cm. 

_  single  patch 

xxxxxxxx  s=0 .038  cm  (region  1) 

oooooooo  s=0.305  cm  (region  2) 

.  $=0,648  cm  (region  3) 
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Fig.  3  E-plane  patterns  for  (a)  single  patch,  (b)  two-layer 
EMCP  antenna  and  (c)  three-layer  EMCP  antenna. 
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ABSTRACT: 


When  a  microstrip  antenna  is  mounted  on  a  conical 
surface,  the  planar  element  geometry  which  conforms  to 
that  of  the  cone  is  the  annular  sector  element.  This 
paper  presents  analytical  and  experimental  results  for 
the  annular  sector  microstrip  antenna  when  mounted  on 
both  planar  and  conical  surfaces.  The  cavity  model  is 
applied  to  analyze  the  resonant  frequency,  input 
impedance,  and  radiation  patterns  of  the  element  on  a 
planar  surface.  This  theory  is  then  expanded  upon  and 
applied  to  the  design  of  elements  mounted  on  conical 
surfaces.  Experimental  pattern  characteristics  of  the 
element  mounted  on  a  conical  surface  are  presented,  for 
both  the  TM01  and  TM^  modes.  A  set  of  design  curves 
has  been  developed  to  provide  a  practical  means  of 


specified  mode  and  resonant  frequency. 


«  1  **>  /->  VN 

wuocu 


1.0  INTRODUCTION: 


The  microstrip  antenna  has  proven  to  be  a  very 
useful  element  in  many  applications,  because  it  may  be 
designed  in  a  variety  of  geometric  shapes,  and  it  can 
be  conformally  mounted  on  flat  and  curved  surfaces. 
Extensive  analysis  has  been  done  for  the  more  commonly 
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used  microstrip  antenna  geometries,  including  the 
rectangle,  3quare,  and  the  circular  disk1"3.  However, 
when  a  patch  element  is  conformally  mounted  on  a 
conical  surface,  the  geometrical  constraints  require 
that  the  annular  sector  shape  be  used.  On  a  planar 
surface,  the  annular  sector  element  is  defined  as  shown 
in  Figure  1.  Although  some  analysis  has  been  done  for 
this  element  in  the  planar  configuration4'5,  very 
little  work  has  been  reported  in  the  literature 
regarding  analysis  and  practical  considerations  for  the 
design  of  such  an  element  mounted  on  conical  surfaces. 

This  paper  applies  the  cavity  model  theory  to 
determine  the  resonant  frequency,  polarization,  input 
impedance,  and  radiation  patterns  for  the  planar 
annular  sector  element.  This  theory  is  extended  to  the 
design  of  elements  on  conical  surfaces  by  projecting 
the  dimensions  of  the  conical  structure  onto  a  planar 
surface,  using  the  cavity  model  to  design  elements  on 
this  equivalent  planar  surface,  and  computing  the 
element  dimensions  as  projected  onto  the  cone. 

2.0  THEORY l 

The  cavity  model  theory  was  originally  developed 
by  Y.  T.  lx?  et  al1  for  application  to  a  wide  variety  of 


338 


canonical  shaped  elements.  Basically,  this  theory 
models  the  microstrip  patch  element  as  a  cavity  bounded 
by  electric  conductors  on  the  top  and  bottom,  and 
magnetic  conductor  walls  on  all  sides.  Restricting  the 
analysis  to  that  of  a  "thin"  substrate  allows  the  field 
components  to  be  modelled  as  constant  along  the  z-axis, 
with  the  E-field  in  the  cavity  oriented  along  the  z- 
axis  (TMZ) ,  By  applying  the  appropriate  boundary 
conditions  in  a  cylindrical  coordinate  system,  this 
theory  may  be  applied  specifically  to  the  annular 
sector  element. 

Using  such  a  cavity  model,  a  series  of  (TM)  modal 

solutions  may  be  obtained.  Since  this  element  is 

analyzed  in  the  cylindrical  coordinate  system,  the 

solutions  are  expressed  in  terms  of  circumferential 

field  variations  as  TM _ _  modes,  with  the  subscripts 

f  H 

defined  ass 

p  ■  number  of  circumferential  field  variations 

(P  m  0,  2,  ...) 

q  “  the  qth  radial  root  for  a  given  p 
(q**  1,  2,  3,  •••) 

To  focus  on  the  practical  applications  for  this 
•lament,  the  analysis  presented  here  shall  concentrate 
on  the  design  of  two  specific  modes,  the  TM01  and  TM-^ 
modes.  The  fundamental  difference  between  these  two 
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modes  is  that  their  internal  field  distributions,  and 
thus  their  radiation  patterns,  are  polarized  on 
orthogonal  axes.  As  shown  in  Figure  ?,  the  E-field  has 
a  radial  component  (Ep  )  for  the  TM01  mode,  and  a 
circumferential  component  (E^ )  for  the  TM-^  mode.  This 
provides  a  great  deal  of  flexibility  when  implementing 
the  annular  sector  element  in  specific  design 
applications . 

2 . 1  Resonant  Frequency 

The  resonant  frequency  is  computed  by  analyzing 
the  internal  field  structure  of  the  element,  subject  to 
the  appropriate  boundary  conditions.  The  solution  is 
derived  in  terms  of  a  wave  function  ^pq:  (E2  =~jufxip0q) 

tf'pq  -  t  AlJm(*P>  +  BlYra(kP>  3  cos(m0)  (1) 

where  Jm  and  are  Bessel  functions  of  the  first  and 
second  kind,  respectively,  of  order  m.  Applying  the 
boundary  conditions  on  the  electric  and  magnetic  fields 
leads  to  a  characteristic  equation: 

JA(Ka)*A<*b)  “  JA(kb)Y;(ka)  =  0  (2) 

m  =  p7 T/pQ,  p  =  o,  i,  2,  ... 

Any  annular  sector  element  design  must  satisfy  (2)  . 


Thus,  given  the  physical  dimensions  (a,  b,  and  0O)  of 
the  element,  the  resonant  frequency  may  be  computed  by 
solving  (2)  for  k,  where  k  =  2n£j  u0€0er  -  However, 
many  design  problems  specify  the  frequency  of 
operation,  and  desire  to  determine  the  element 
dimensions.  This  can  be  accomplished  best  by  casting 
(2)  in  another  form.  Let  us  express  the  characteristic 
equation  (2)  as  follows: 

YA<ka)  =  *£<“>_>.  o) 

J^(ka)  J^(kb) 

In  this  form,  the  solution  is  obtained  for  a  given 
frequency  (which  fixes  k)  and  angle  fbQ  (which  fixes  m) 
by  determining  the  inner  (a)  and  outer  (b)  radii  that 
satisfy  the  equation. 

The  practical  application  of  this  may  be  seen  by 
plotting  x^(kr)/J^(kr)  vs  kr  for  the  desired  modes. 
The  resultant  graphical  solutions  provide  a  set  of 
design  tools  for  determining  the  dimensions  of  an 
annular  sector  element.  The  characteristic  curves  for 
the  TMQ1  mode  are  shown  in  Figure  3 .  Note  that  the 
solutions  for  the  TM01  mode  are  independent  of  the 
sector  angle  A  set  of  TM^  mode  solutions  are 

shown  in  Figure  4  for  a  number  of  sector  angles  ranging 
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ANNULAR  SECTOR  PATCH  ANTENNA 

7 M  Ot  MODE  SOLUTIONS,  FOR  ANY  ANGLE  % 


from  10  to  60  degrees.  These  design  curves  are  applied 
by  choosing  the  desired  sector  angle  pQ,  and  finding 
two  points,  ka  and  kb,  which  have  equal  values  on  the 
vertical  axis. 

The  dimensions  calculated  with  this  method  do  not 
account  for  the  fringing  fields  which  extend  beyond  the 
physical  dimensions  on  the  radiating  edges.  Thus,  to 
maintain  resonance  at  a  given  frequency,  the  element 
dimensions  must  be  adjusted  to  compensate  for  the 
fringing  effects.  Formulas  have  been  developed  for 
various  element  shapes2' 6,7 .  For  the  annular  sector 
element  the  microstrip  line  open  circuit  formula8  was 
modified  to  fit  more  closely  with  the  measured  data: 

£1  =  0. 412h  (fr  +  0 « 300)  (Weff/h  +  0.264)  (4) 
(£r  -  0.258)  (Weff/h  +  0.800) 

Thus,  (4)  has  modified  the  original  formulation  by 
using  the  relative  dielectric  constant  er,  rather  than 
an  effective  dielectric  constant.  This  method  was 
found  to  provide  a  more  accurate  prediction  of  the 
resonant  frequency,  similar  to  observations  made  by 
others  using  different  element  shapes6.  weff  an 
effective  element  width,  defined  below. 

Using  Al,  the  physical  dimensions  of  the  element 
may  be  adjusted  to  account  for  the  fringing  effects. 


The  resonant  length  of  the  element  depends  upon  the 
selected  mode.  Shown  in  Figure  2,  the  TM01  resonant 
length  is  on  the  radial  axis,  whereas  the  TM1;L  resonant 
length  is  on  the  angular  axis.  Thus,  the  element 
corrections  to  account  for  fringing  in  these  two  modes 


are  derived  as: 

TMg1  Mode: 

weff  =  9>o  <5a> 

2 

a'  =  a  +  A1  (5b) 

b'  =  b  -  Al  (5c) 

P'o  *  0O  (5d) 

TMi:l  Mode: 


weff  "  b  "  a 
a '  -  a 


b' 

V'o 


b 


4A1 
a  +  b 


(6a) 

(6c) 

(6d) 


The  primed  variables  in  (5)  and  (6)  represent  the 
corrected  physical  dimensions  of  the  element  after 
accounting  for  the  fringing.  (Note:  angular  variables 
are  expressed  in  radians.) 
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2.2  Input  Impedance 


The  cavity  model  may  be  used  to  calculate  the 
patch  input  impe  ance  using  a  coax  probe  feed. 
However,  the  initial  model  must  be  modified  to  account 
for  the  feed  location  in  the  cavity,  and  to  properly 
include  the  power  loss  due  to  radiation  and  dissipation 
mechanisms.  Basically,  the  input  impedance  is 
proportional  to  the  E-field ,  as  given  by: 


'in 


1=1 


amp 


(?) 


Thus,  the  probe  feed  is  modelled  as  a  1  amp  current 
source  located  at  (  pf ,  0j)  .  To  account  for  this 
discontinuity,  the  analysis  was  modified  to  an 
inhomogeneous  problem  using  a  technique  of  mode 
matching.  This  technique  allows  the  patch  element  to 
be  divided  into  two  homogeneous  regions,  each  of  which 
satisfies  the  source-free  boundary  conditions.  The  sum 
of  these  modes  are  then  matched  to  the  source 
discontinuity.  The  wave  function  given  by  (1)  is 

expanded  about  the  boundary  at  p  =  p^: 

CD 

^pq  =  J  [AiJm(*p)  +  BpYjjj ( kp )  ] cos (m£S )  (8a) 

p=0  (a  <  p  <  pf) 

CD 

tf'pq  “  y  [ApJm(^P)  +  BpYro(kp) ]cos(mp)  (8b) 

^  p^O  (Pf  <  P  <  b) 
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The  wave  functions  in  (8)  were  expanded  in  a  series 
form  to  include  the  reactive  impedance  contributions 
from  all  higher  crder  modes.  These  formulations  were 
solved  by  enforcing  field  continuity  between  the  two 
regions.  Thus,  the  input  impedance  solution  was 
derived  and  found  to  agree  with  that  by  Ou^ : 


'in 


CO 

2 

p=0 


[Ja(kpf)Yi(ka)  -  Jffi(ka) Ym(kpf ) ] 


[ Jm(kb) Y^(ka)  - 


Ym(kb) J^(ka) ] 


•  CYi(kb)Jn(kPf)  -  j;(kb)Ym(kpf)]cgs2j^.gfl 


1  +  <5 


Op 


(9) 


where  <50p  ■»  1  for  p  =  0,  and  0  for  p  +  0.  However,  (9) 
does  not  yet  account  for  the  power  losses  due  to 
radiation  (Prad)  ,  and  dissipation  in  the  dielectric 
(Pd)  and  copper  plating  (Pc) .  Richards,  et  al9 
developed  an  effective  loss  tangent  deff?  equal  to  1/Q 
of  the  cavity; 


fieff  =  1  *  prad  *  pd  *’  pc  (10) 
Q  2  we 

Thus,  the  power  terms1  and  the  stored  electric  energy 
We10  are  computed  and  used  to  develop  this  effective 
loss  tangent.  This  term  for  0eff  is  then  used  to 
express  the  wave  number  k  as  a  complex  quantity: 
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(11) 


keff  “  27rf  /  (1  “  ^eff^ocofr 

Finally,  the  complex  wave  number  in  (11)  is  substituted 
into  (9)  to  provide  an  accurate  model  for  computing  the 
input  impedance  of  the  annular  sector  element. 

2.3  Radiation  Patterns 

The  radiation  characteristics  of  this  element  were 
computed  using  the  equivalence  principle.  The  electric 
fields  at  all  edges  of  the  element  may  be  derived  from 
(8)  where  Ez  =  -j  Huygen's  principle  was  then 
applied  to  obtain  the  equivalant  magnetic  current 
densities  on  each  of  the  cavity  walls: 

H  =  -2n  x  5ZEZ  (12) 

By  integrating  these  magnetic  currents  along  the  four 
edges  of  the  element,  the  far-field  radiation  patterns 
were  computed.  However,  to  account  for  the  edge 
effects  of  the  finite  substrate  ground  plane,  the 
Geometric  Theory  of  Diffraction  (GTD)11  was  applied. 

The  radiation  pattern  theory  was  again  applied 
only  to  the  annular  sector  element  mounted  on  a  planar 
surface,  although  recent  work  has  been  done  to  develop 
theory  for  evaluation  of  pattern  characteristics  on 
more  complex  conformal  surfaces12. 
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3.0  EXPERIMENTAL  RESULTS: 


The  theory  developed  for  the  planar  annular  sector 
patch  was  applied  to  design  both  TM01  and  TM^  planar 
elements.  Furthermore,  to  verify  the  application  of 
this  theory  to  antenna  design  on  conformal  surfaces, 
the  planar  dimensions  were  used  to  compute  equivalent 
element  dimensions  as  projected  onto  a  cone.  The  cone 
structure  fabricated  for  these  tests  had  a  30°  apex 
angle  and  was  truncated  at  both  ends,  as  indicated  in 
Figure  2.  Th®  antennas  were  fabricated  on  a  0.062" 
thick  substrate  with  a  dielectric  constant  € r  =  2.2. 
The  design  frequency  was  chosen  as  fQ  «  2.5  GHz . 

The  characteristic  equation  (3)  and  fringing  terms 
(4) -(6)  were  used  to  compute  the  element  dimensions. 
However,  when  optimizing  the  input  impedance,  it  was 
observed  that  the  impedance  formula  (9)  provided  a  more 
accurate  prediction  of  the  resonant  frequency. 

The  resonant  frequency  measurements  are  compared 
with  theory  in  Table  l  for  each  of  the  element  designs. 
The  resonant  frequency  was  defined  as  the  frequency 
which  has  the  minimum  VSWR,  or  optimum  impedance  match. 


TABLE  1:  RESONANT  FREQUENCY: 
MEASURED  VS.  THEORETICAL 


Element 

Measured 

Theoretical 

Theoretical 

Description 

(GHz) 

Using  Eq.  (3) 

Using  Eq.  (9) 

(%  error) 

(%  error) 

TM01  Planar 

2.539 

2.500 

2.522 

(1.5%) 

(0.7%) 

TM01  Conical 

2.527 

2.500 

2.522 

(1.1%) 

(0.2%) 

TM^  Planar 

2.535 

2.500 

2.521 

(1.4%) 

(0.6%) 

TM1-^  Conical 

2.452 

2.500 

2.521 

(1.9%) 

(2.8%) 

The  input  impedance  was  optimized  by  applying  (9) , 
and  locating  the  feed  point  to  maximize  the  2 : 1  VSWR 
bandwidth.  The  measured  and  theoretical  impedance  of 
the  four  antenna  elements  are  shown  in  Figures  5-8 .  As 

ran  hp  coon  fho  aaraanont1  4e  nonoral  lw  frrrirl  uith  fhn 

exception  of  the  TM01  conical  element  in  Figure  6. 
This  element  has  good  agreement  for  the  resonant 
frequency,  but  shows  a  discrepancy  in  the  impedance 
magnitude,  which  is  believed  to  be  caused  by  the  close 
proximity  of  the  patch  radiation  edges  to  the  truncated 
cone  edges. 

The  radiation  patterns  were  measured  for  the  TM01 
and  TM^  conical  annular  sector  elements,  and  are  shown 
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FIGURE  8: 

INPUT  IMPEDANCE  FOR  TMjj  CONICAL  ELEMENT 
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in  Figures  9  and  10,  respectively.  The  pattern 
characteristics  for  these  two  elements  are  summarized 
in  Table  2: 

TABLE  2:  MEASURED  RADIATION  PATTERN  CHARACTERISTICS 
FOR  CONICAL  MOUNTED  ANNULAR  SECTOR  ELEMENTS 

Element  E-Plane  K-Plane  Peak  Gain 
Description  HPBW  HPBW  (dBi) 


™01 

Mode 

90  Deg. 

94  Deg. 

6.8  dBi 

™11 

Mode 

86  Deg. 

7  5  Deg . 

5.5  dBi 

One  discrepancy  to  be  noted  is  that  the  TM01  patterns 
have  an  H-plane  half-power  beamwidth  (HPBW)  which  is 
widgr  than  the  E-plane  HPBW.  Although  this  is  not 
typical  for  a  planar  microstrip  antenna,  it  may  be 
explained  by  referring  to  Figure  2(a),  observing  that 
the  radiating  edges  of  this  element  are  conformally 
curved  around  the  circumfex-ence  of  the  cone,  and  thus 
radiating  over  a  wider  azimuth  range  than  a  planar 
element.  Overall,  the  radiation  patterns  indicate  that 
the  annular  sector  element  may  provide  good  coverage 
when  mounted  on  a  conical  surface. 


FIGURE  9 


(a)  ELEVATION  PLANE  (E-PLANE) 


(b)  AZIMUTH  PLANE  (H-PLANE) 

TM  CONICAL  ELEMENT  RADIATION  PATTERNS 
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(a)  ELEVATION  PLANE  (H-PLANE) 


(b)  AZIMUTH  PLANE  (E-PLANE) 

'  I  CUKE  10:  TMn  CONICAL  ELEMENT  RADIATION  PATTERNS 
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4.0  CONCLUSIONS!  AND  FURTHER  APPLICATIONS: 


This  paper  has  applied  the  cavity  model  theory  to 
analyse  the  annular  sector  patch  on  planar  and  conical 
surfaces.  To  provide  a  practical  design  method,  a  set 
of  design  curves  was  presented  for  the  TMQ1  and  TM13_ 
modes.  Accurate  formulas  were  developed  to  properly 
account  for  the  edge  fringing.  Experimental  results 
confirmed  the  accuracy  of  the  planar  element  theory 
when  applied  to  the  design  of  elements  on  conical 
surfaces. 

Further  applications  being  pursued  are  the  use  of 
annular  sector  elements  in  the  design  of  conical 
microstrip  arrays,  and  the  extension  of  this  theory  to 
the  design  of  circular  polarized  (CP)  elements. 
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MODE  PURITY  OF  RECTANGULAR  PATCH  ANTENNAS 
WITH  POST  AND  APERTURE  EXCITATIONS 
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Abstract 


micrustnp  putcii  uiucium^>  oiici  n’^ny  <iitr active  lCatures  WuiCu  iiiCiUuC  icw  proii»c, 
low  cost,  confomiability,  compactness,  and  light  weight.  Even  the  inherent  high  quality 
factor  of  these  antennas  can  be  advantageous  if  a  narrow  band  of  operation  is  desired. 

The  microstrip  patch  antenna  can  be  excited  by  various  methods  such  as  slot  line, 
coplanar  waveguide,  microstrip  line,  apertures,  and  posts.  Each  of  the  methods  above  has 
advantages  and  disadvantages  in  the  areas  of  fabrication,  spurious  radiation,  feed 
network/antenna  design,  and  input  impedance.  It  is  the  focus  of  this  paper  to  investigate 
the  post  feed  and  aperture  coupling  methods  of  excitation  for  identical  rectangular 
microstrip  patches  of  various  aspect  ratios.  Comparison  between  the  two  feeds  will  center 
around  the  issue  of  mode  purity,  in  this  case,  the  0,1  mode. 


To  experimentally  detect  the  excitation  of  modes  other  than  the  0,1  mode,  properties 
of  the  far-field  radiation,  in  particular,  cross-polarization  in  the  principal  planes,  will  be 
used. 
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Introduction 

For  linear  polarization,  it  is  usually  desirable  that  the  microstrip  patch  antenna  be 
excited  in  a  single  mode.  Associated  with  the  desired  mode  is  a  particular  radiation  pattern. 
Any  additional  modes  which  are  excited  by  the  feed  mechanism  may  result  in  additional 
radiation  which  is  likely  to  be  undesirable.  Hence,  some  types  of  degradation  of  the 
dominant  mode  radiation  pattern  can  be  attributed  to  spurious  modes  of  the  patch. 

One  detrimental  effect  oi  mode  contamination  is  a  nse  in  the  cross-polarization 
levels.  F'0‘  t.  rectangular  microstrip  patch  antenna  operating  in  the  0,1  mode,  the  cross- 
polarization  in  the  principal  planes  is,  theoretically,  zero.  Here  we  use  Ludwig's  second 
definition  ;  ij  :V  r  -ross-polarization  and  defu-e  E$  as  the  •  o-pol  component  and  Lq  as  cross- 
pol  component  for  the  H-plane.  In  the  E-plane,  the  reverse  is  tr  .  Eg  becomes  the  co-pol 
component  and  E<j>  becomes  the  cross-pol  component.  For  the  0,1  mode,  cross¬ 
polarization  in  either  principal  plane  can  be  attributed  to  the  excitation  of  additional  modes. 

Oberhart  and  Lo  previously  studied  the  cross-polarization  of  rectangular  patch 
antennas.  In  their  study,  the  cavity  model  theory  [4]  for  the  rectangular  microstrip  patch 
antenna  was  employed  to  determine  a  lECp0xl/IExpol'  ratio  for  a  given  feed  location  and 
various  aspect  ratios  ranging  from  1.00  to  2.1 1.  R.  C.  Hansen  also  discusses  the  cross¬ 
polarization  of  square,  rectangular,  and  circular  microstrip  patches  [5].  Rectangular 
patches  were  the  subject  of  investigation  by  Huynh,  Lee,  and  Lee  [6]  in  which  the  cavity 
model  theory  was  used  to  determine  relative  co-pol  to  cross-pol  levels  as  a  function  of  not 
only  aspect  ratios  but  also  feed  position,  substrate  thickness,  and  relative  permitivity  of  the 
substrate  supporting  the  patch. 

The  scope  of  this  paper  is  to  compare  the  purity  of  the  0,1  mode  for  identical 
rectangular  microstrip  patch  antennas  where  one  is  excited  by  a  post  and  the  other  is  excited 
by  an  rectangular  aperture.  This  comparison  will  be  conducted  for  various  aspect  ratios 
ra:  ging  from  1.5  to  3.5  for  both  feed  types. 
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Construction  and  Measurements 

The  geometries  of  the  two  patch  antennas  are  shown  in  Figures  1  and  2.  The  post- 
fed  patch  (Figure  1)  consists  of  a  copper  rectangle  of  appropriate  dimensions  which  is 
supported  2mm  above  a  ground  plane  by  Rohacell ,  a  low  loss,  low  permitivity  foam  (£r 
=1,071  and  a  loss  tangent  of  ,00156).  The  post,  which  is  the  center  conductor  of  ,085 
coaxial  cable,  is  located  at  (x=0.0mm,  y=l0.0mm).  The  post  at  this  position  will  couple 
strongly  into  the  electric  field  of  the  0,1  mode,  less  strongly  to  that  of  other  modes. 

The  aperture-fed  patch  (Figure  2)  is  configured  identically.  However,  instead  of 
the  post  feed,  a  rc  Hang-.-lar  aperture  (,5mm  by  20mm)  is  cut  in  the  ground  plane.  This 
anenute  <'j  located  0mm,  y=25.0mm)  and  will  cou.  V  strongly  to  the  magnetic  field 

of  the  3,1  mode  but  will  ot  couple  to  as  many  higher  order  modes  as  the  post  because  of 
its  symmetric  location.  A  microstrip  line  which  runs  along  the  underside  of  the  ground 
plane  separated  by  a  1/16"  of  Rexolite  2200  dielectric  is  used  to  excite  the  siot. 

Scattering  parameters  were  measured  on  a  HP8510  Network  Analyzer  to  determine 
the  exact  resonant  frequency  of  the  0,1  mode.  Both  Through-Line-Reflect  calibration  and  a 
standard  3-term  calibration  were  used  for  error  correction. 

At  the  appropriate  resonant  frequencies,  far-field  iadiation  patterns  were  taken  for 
l>oth  polarizations  in  both  principal  planes  and  one  diagonal  plane.  The  measurements  were 
conducted  on  a  20'  by  20'  ground  plane  range.  The  receiver  used  was  a  Scientific  Atlanta 
model  which  has  a  40  dB  dynamic  range.  Great  care  was  taken  to  ensure  proper 
orientation  of  the  microstrip  patch  antenna  and  the  alignment  of  the  receiving  horn  antenna. 

Results  and  Discussion 

The  original  dimensions  of  the  rectangular  patches  were  a  =  175mm  by  b  =  50mm 
which  gives  an  aspect  ratio  of  3.5.  The  aspect  ratios  of  the  patches  were  varied  from  3.5  to 
1.5  in  0.5  increments  by  merely  trimming  the  "a"  dimension  (while  keeping  the  "b" 
dimension  constant).  Although  the  cavity  model  theory  predicts  the  0,1  resonance  to  be 
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3.00  GHz,  due  to  the  presence  of  the  fringing  fields  at  the  edge  of  the  patch,  the  actual 
resonance  frequency  is  depressed.  The  amount  by  which  the  "resonant  frequency"  is 
perturbed  is  also  a  function  of  aspect  ratio.  In  addition  to  the  effect  of  the  fringing  fields, 
the  resonant  frequency  of  the  aperture-  fed  patch  is  decreased  because  of  the  presence  of  the 
slot.  Scattering  parameters  were  measured  to  determine  the  exact  0,1  resonant  frequency. 
Far-field  radiation  patterns  were  measured  at  the  exact  resonant  frequency.  Although  far- 
fieid  data  was  measured  in  both  principal  planes  and  in  one  diagonal  plane,  the  most 
meaningful  data  was  gathered  in  the.  H-plane. 

Figure  3  shows  the  H-plane  far-field  radiation  patterns  for  1.5  aspect  rario 
rectangular  patches  fed  both  by  post  and  aperture.  The  cross-polarization  component  b.0  is 
relatively  non-existent  for  the  aperture-fed  patch.  However  for  the  post-fed  patch,  cross 
polarization  levels  as  high  as  -23  dB  are  observed.  The  shape  ol  the  Ho  pattern  suggests 
that  the  2,0  mode  is  being  slightly  excited.  Nothing  unusual  is  noted  for  the  component 
except  for  a  ripple  in  the  pattern  at  zenith  which  is  believed  to  be  the  result  of  multiple 
scattering. 

Data  for  the  2.0  aspect  ratio  patches  are  presented  in  Figure  4.  In  the  post-fed 
patch,  the  2,0  mode  presence  can  clearly  be  detected  in  the  pattern  for  the  Hq  component. 
This  is  expected  since  the  cavity  model  theory  predicts  that  the  0,1  and  2,0  modes  have  the 
same  resonant  frequency  at  this  aspect  ratio.  Again  the  cross-polarization  level  for  the 
aperture-fed  patch  is  considerably  lower,  with  peaks  in  the  neighborhood  of  -35  dB. 

As  the  aspect  ratio  is  increased  to  2.5  (Figure  5),  the  aperture-fed  patch  is  again 
well  behaved  with  respect  to  cross-polarization  performance.  However  in  the  case  of  the 
post-fed  patch,  the  cross-polarization,  with  peaks  at  the  -23  dB  level,  still  indicate  the 
presence  of  the  2,0  mode.  The  co-polarization  component,  E<j>,  begins  to  show  side  lobes 
at  0  =  70’  that  are  33  dB  down.  These  side  lobes  result  from  the  extreme  length  of  the 
principal  radiating  magnetic  current. 
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At  aspect  ratio  of  3.0  (Figure  6),  curiously  enough,  both  post  and  aperture  feeds 
produce  almost  non-existent  cross-polarization.  The  E^  component  of  the  post-fed  patch 
shows  that  the  side  lobe  is  now  well  defined  and  levels  are  now  peaked  at  0  =  60  and  have 
risen  to  the  -20  dB  level.  Side  lobes  also  begin  to  appear  in  the  apertuie-fed  patch  at  the 
-31  dB  level  for  6  =  70'. 

The  final  aspect  ratio  measured  was  3.5  (Figure  7).  The  component  for  both  the 
aperture  and  post  feeds  show  that  the  side  lobe  level  is  around  -23  dB  and  that  the  side 
lobes  are  not  as  well  defined  as  in  the  case  of  the  post-fed  patch  of  aspect  ratio  3.0.  The  Ee 
component  of  the  post-fed  patch  reveals  that  the  4,0  mode  is  being  excited.  The  cross- 
polarized  field  for  the  aperture-fed  patch  peaks  around  the  -35  dB  level. 

ivFijtyu 

It  has  been  shown  that  one  aperture  feed  produces  a  purer  0,1  mode  excitation  than 
a  particular  post  feed  over  a  variety  of  aspect  ratios  for  rectangular  patch  antennas.  Only 
one  aperture  dimension  was  investigated  and  no  attempt  was  made  to  vary  the  input 
impedance  by  manipulation  of  the  feeds.  It  is  quite  certain  that  the  mode  purity  resulting 
from  the  aperture  excitation  is  a  function  of  both  aperture  feed  dimension  and  position. 
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A  p  e  r  t  u  r  c  -  f  e  o!  patch 


Figure  2  Geometry  of  the  aperture-fed  rectangular  patch. 
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Figure  4  H-plane  of  both  the  post-fed  and  aperture-fed  rectangular  patches  of  aspect  ratio  2.0. 
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Figure  5  H-plane  of  both  the  post-fed  and  aperture-fed  rectangular  patches  of  aspect  ratio  2.5. 
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CIRCULAR  POLARIZATION  FROM  STACKED  MICROSTRIP  PATCHES 


by 

Andrew  Adrian 
Texas  Instruments 
McKinney,  Texas 

Daniel  H.  Schaubert 
University  of  Massachusetts 
Amherst,  Massachusetts 

ABSTRACT 

A  novel  method  of  achieving  circularly  polarised 
radiation  with  rectangular  laicrostrip  patches  is 
described.  In  this  geometrical  arrangement,  two 
stacked  microstrip  patches  are  orthogonally  polarised. 
The  bottom  patch  acts  as  a  partial  ground  plane  for 
the  top  patch,  and  the  top  patch  does  not 
significantly  affect  the  radiation  from  the  bottom 
patch.  Energy  is  coupled  to  the  top  patch  via  a 
diagonal  aperture  in  the  center  of  the  bottom  patch 
while  the  entire  system  is  fed  by  a  single  excitation 
of  the  bottom  patch.  This  excitation  can  consist  of 
any  suitable  feed  for  a  microstrip  patch.  In  this 
case,  a  probe  feed  was  utilized.  Circularly  polarized 
radiation  then  occurs  when  the  aperture  is  tuned  for 
3  dB  coupling  wi'h  a  90  degree  phase  shift. 
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(Fig.  la).  Consequently,  the  top  patch  does  not 
significantly  affect  the  radiation  of  the  bottom 
patch,  and  the  bottom  patch  acts  as  a  partial  ground 
plane  for  the  top  patch. 

2.2  Excitation 

Only  the  bottom  patch  needs  an  external 
excitation.  This  can  be  implemented  as  any  suitable 
microstrip  patch  antenna  feed.  in  the  case  that  is 
being  described  here,  a  probe  was  utilized. 
Consequently,  the  top  patch  is  excited  by  energy  that 
is  coupled  through  the  aperture  from  the  bottom  patch 
( Fig.  2 ) . 

The  aperture  itself  is  centered  on  the  bottom 
patch.  However,  its  major  axis  is  oriented  on  a  45 
degree  angle  to  both  patches.  In  order  to  achieve 
circularly  polarized  radiation,  the  aperture  must  be 
tuned  for  3  dB  coupling  with  a  90  degree  phase  shift. 
This  can  be  accomplished  by  adjusting  the  aperture 
size  and  .its  exact  shape. 

3.0  EXPERIMENTAL  RESULTS 

3.1  Geometry 

An  experimental  prototype  of  the  ACSCFP  was 
fabricated  with  two  layers  of  0,062  inch  thick  Rogers 
R/T  Duroid  5880  which  has  a  dielectric  permittivity  of 
of  2.20.  Dimensions  of  the  exact  geometry  are 
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displayed  in  Fig.  3.  Amplitude  and  phase  tuning  of 
the  aperture  resulted  in  its  unique  shape. 

3.2  Center  Frequency 

The  operating  center  frequency  of  the  element  is 
3.470  GHz  where  its  return  loss  is  6.8  dB .  This  is 
evident  from  the  cusp  in  the  Smith  Chart  and  return 
loss  plots  of  the  ACSCPP  (Fig.  4)  as  well  as  the 
measurements  of  axial  ratio  over  frequency  (Table  1). 

A  spin  -linear  radiation  pattern  that  was  taker,  in  the 
E-plane  of  the  bottom  patch  at  3.470  GHz  is 
illustrated  in  Fig.  5.  The  axial  ratio  at  this 
frequency  is  1.8  dB . 

3.3  Bandwidth 

The  operating  bandwidth  of  this  element  is 
limited  by  its  axial  ratio.  Below  3.450  GHz  and  above 
3,490  GHz,  the  axial  ratio  is  greater  than  3  dB  (Table 
1).  This  resulted  in  a  3  dB  axial  ratio  bandwidth  of 
1.01  %.  From  measurements  of  the  relative  radiated 
power  of  each  patch  in  the  stack  and  the  axial  ratio 
of  the  combined  patches,  the  phase  shift  through  the 
aperture  was  calculated  (Table  1).  From  Table  1  it  is 
evident  that  at  the  low  end  of  the  element's  frequency 
band,  a  poor  phase  shift  through  the  aperture 
contributes  to  the  poor  axial  ratio  there. 

Conversely,  on  the  uihet  end  of  the  band,  a  poor  power 
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split  limits  the  axial  ratio  bandwidth. 

The  antenna  tested  here  has  been  compared  to  a 
single  patch  with  a  diagonal  slot  for  circular 
polarization  (Kerr)  [1].  The  Kerr  antenna  was 
fabricated  on  a  substrate  of  the  same  dielectric 
permittivity  as  the  ACSCPP.  Conversely,  the  substrate 
was  0.128  inch  thick  in  order  to  have  the  same  overall 
thickness  as  the  ACSCPP.  This  was  felt  to  be  the  most 
fair  comparison  of  antenna  bandwidths.  However,  the 
ACSCPP  require  that  only  the  bottom  patch  be  fed, 
which  means  the  feed  probe  or  microstrip  feedline 
corresponds  to  a  substrate  one-half  the  total  antenna 
thickness.  This  can  be  very  important  for  wide 
bandwidth  antennas  on  thick  substrates  where  feed 
inductance  and/or  radiation  can  cause  severe  problems. 

The  3  dB  axial  ratio  bandwidth  of  the  Kerr 

antenna  was  optimized  to  a  value  of  1.46  %,  which  is 

greater  than  the  prototype  ACSCPP.  However,  the 

ACSCPP  may  be  capable  of  greater  bandwidth  if  an 

optimized  design  is  developed.  This  development  can 

utilize  the  following  interdependent  parameters  to 

obtain  the  most  desirable  performance: 

Bottom  patch  width 
Aperture  size  and  shape 
Top  patch  width 

Substrate  thicknesses  and  permittivities. 
During  the  development  of  the  prototype  ACSCPP,  it  was 
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found  that  the  resonant  frequencies  of  the  top  and 
bottom  patches  are  affected  differently  by  changes  in 
the  aperture  size.  Therefore,  trial  and  error  design 
is  tedious  and  an  analysis  of  the  element  would  be 
very  benificial. 

4.0  CONCLUSION 

This  paper  has  documented  a  proof  of  concept 
design  for  ACSCPP.  The  input  return  loss  was  6.8  dB, 
and  the  3  dB  axial  ratio  bandwidth  was  1.01  %. 
Improving  the  input  match  in  this  case  can  probably  be 
accomplished  by  moving  the  feed  probe  closer  to  a 
•radiating  edge  of  the  bottom  patch.  The  additional 
distance  between  the  aperture  and  the  probe  may  also 
result  in  an  easier  tuning  procedure.  Gains  in  axial 
ratio  bandwidth  may  be  possible  by  optimizing  the 
aperture  for  an  improved  amplitude  and  phase  coupling 
profile.  Furthermore,  use  of  the  ACSCPP  requires 
feeding  only  the  bottom  patch.  This  may  be 
advantageous  for  antennas  on  very  thick  substrates  (as 
when  wide  bandwidth  is  desired). 

Finally,  the  ACSCPP  may  have  a  potential  which 
has  yet  to  be  investigated.  Since  the  radiators  are 
stacked,  materials  with  different  dielectric 
permittivities  and/or  thicknesses  can  be  used  for  each 
layer  of  the  fabrication.  This  could  provide  a  means 
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for  obtaining  dual  frequency  circular  polarization 
with  a  single  element  or  circular  polarization  at  one 
frequency  and  linear  at  one  or  two  others.  This  would 
require  resonating  the  patches  in  both  linear 
directions. 
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Figure  1 


(a)  Top  and  (b)  side  views  of  stacked, 
orthogonally  polarised  microstrip  patches. 
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Figure  2  (a)  Top  and  (b)  side  views  of  stacked, 

orthogonally  polarized  microstrip  patches 
with  an  aperture  coupled  top  patch. 
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Figure  3  Dimensions  of  the  prototype  ACSCPP.  The 

(a)  patches  are  dimensioned  in  cm  and  the 

(b)  aperture  is  labeled  in  mm. 


Figure  4  (a)  Input  im; 
ACSCPP. 
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Abstract 


A  rigorous  analysis  of  the  natural  resonance  frequencies  and  input  impe¬ 
dance  characteristics  of  an  annular-ring-loaded  (ARL)  circular-disk  micro¬ 
strip  antenna  is  presented.  Using  vector  Hankel  transforms  (VHT),  the 
problem  is  formulated  in  terms  of  vector  dual-im,egral  equations.  Galerkin’s 
method  is  then  used  to  solve  the  equations  to  obtain  the  resonance  frequen¬ 
cies  and  the  current  distributions  on  the  conductive  patches  arising  from  a 
probe  excitation.  Due  to  the  singular  nature  of  the  current  distribution,  the 
singularity  subtraction  method  has  been  used  to  accelerate  the  convergence  of 
basis  function  expansions.  Experiments  for  determining  resonance  frequen¬ 
cies  and  input  impedance  characteristics  of  an  ARL  circular-disk  micros^rip 
antenna  with  various  substrate  thicknesses  have  been  made.  The  theoreti¬ 
cal  results  are  in  good  agreement  with  the  experiment  data  even  when  the 
electrical  thickness  of  the  substrate  is  0.1  wavelength.  It  is  shown  that  this 
theory  can  be  used  to  analyze  some  microstrip  antennas  with  an  electrically 


thick  substrate  and  mutual  coupling  between  conductive  patches  or  between 
the  patch  and  the  feed  of  a  microstrip  antenna. 


I.  Introduction 

In  recent  years,  microstrip  antennas  have  aroused  great  interest  in  both 
theoretical  research  and  engineering  applications  due  to  their  low  profile, 
light  weight,  conformal  structure,  and  ease  in  fabrication  and  integration 
with  solid-state  devices.  Because  of  the  inherent  narrow  bandwidth  of  mi¬ 
crostrip  antennas,  many  attempts  have  been  made  to  improve  their  frequency 
characteristics.  In  general,  there  are  two  efficient  approaches  to  broaden  the 
bandwidth  of  the  microstrip'- antennas,  i.e.,  by  increasing  the  substrate  thick- 
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ness  and  by  combining  several  patches  with  similar  resonance  frequencies  on 
the  same  or  different  layers  of  an  antenna  structure. 

Unfortunately,  some  popular  analysis  methods  for  microstrip  antennas 
in  the  past,  such  as  the  cavity  model  method,  are  no  longer  valid  in  the 
above  cases  because  these  methods  do  not  consider  the  effects  of  surface 
waves  and  fringing  fields  at  the  t  Ige  of  the  patch.  They  cannot  be  used 
to  model  easily  the  mutual  coupling  between  conductive  patches  or  between 
the  patch  and  the  feed.  Some  efforts  have  been  made  to  analyze  the  above 
problems  [1]- [3] ;  most  of  these  used  the  exact  Green’s  functions  and  spectral 
domain  moment  formulation.  Due  to  the  limitations  of  computer  resources, 
the  numerical  results  for  the  antennas  with  the  thickness  greater  than  0.054 
wavelength  have  not  been  given  in  [1].  In  [2],  an  artificial  approximation  for 
edge  extension  has  been  used  as  a  correction  because  the  analysis  in  [2]  was 
still  based  on  a  cavity  model  even  though  the  substrate  was  electrically  thick. 
The  analysis  in  [3]  did  not  consider  the  singular  nature  of  the  current  at  the 
feed  region. 


In  this  paper,  the  annular-ring  loaded  (ARL)  circular-disk  microstrip 
antenna  with  mutual  coupling  between  the  ring  and  the  disk  has  been  ana¬ 
lyzed  in  the  spectral  domain  by  using  the  vector  Hankel  transform  (VHT). 
First  we  study  the  natural  resonance  frequencies  of  the  ARL  circular- disk 
antenna.  The  problem  is  formulated  in  terms  of  a  set  of  vector  dual  integral 
equations.  GaJerkin’s  method  is  then  used  to  convert  the  integral  equation 
into  a  system  of  linear  algebraic  equations.  The  natural  resonances  can  be 
found  by  searching  for  the  zeroes  of  the  coefficient  matrix.  Good  agreement 
between  theoretical  results  and  experiment  data  has  been  obtained.  Next, 
the  problem  of  an  ARL  circular -disk  microstrip  antenna  excited  by  a  probe 
is  studied.  Due  to  the  singular  nature  of  the  current  distribution  on  the  con¬ 
ductive  patch,  the  singularity  subtraction  method  has  been  used  to  accelerate 
tne  convergence  of  the  basis  function  expansion.  The  induced  EN1F  method 
is  used  to  calculate  the  input  impedance  of  the  probe-excited  antenna  at  the 
ground  plane  of  the  antenna.  Good  agreement  has  been  found  between  the 
numerical  results  and  measured  data  even  for  substrate  thicknesses  up  to  0.1 
wavelength. 


II.  Integral  Equations  for  the  ARL  Disk  Antenna 

Using  the  vector  Hankel  transform  (VHT)  formulation  [1,4,5,7,S],  the 
vector  dual  integral  equations  for  the  ARL  disk  antenna  shown  in  Figure  1 
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are  giver,  by 


e„(p)  =  J dkpkpHn(kpp)  ■  G(kp)  •  K „(kp)  =  0,  p  <  ax,  a  <  p  <  b,  (1) 
o 


OO 

e„(p)  =  J  dkpkpHn(kpp)  •  Kn{kp)  —  0,  aj  <  p  <  a,  p  >  b,  (2) 


where 


e»(p)  = 


EPn(p)/ cos  n<j> 
is  related  to  the  n-th  harmonic  of  the  electric  field, 


-Kpn(p)/ cos  n<f> 
■h',in(p)/sinn<j!> 


*n(/»)  = 

is  related  to  the  n-th  harmonic  of  the  current  on  the  ring  and  the  disk, 


(3) 


(4) 


K“l  ''  “  1  A'n2((.,) 


(5) 


is  the  VHT  of  the  current  Kn(p )  on  the  ring  and  the  disk, 

^■n(kpp')  — 


KM 

jr^Jn{kpp)  Jn(kpp) 


(6) 


is  the  kernel  of  the  vector  Hankel  transform,  and 

rv;  \&(l-R™)  0 

0  f£(l+RTE) 


(7) 


is  the  dyadic  Green’s  function  for  the  p  and  <j>  components  of  the  current  in 
the  VHT  spectral  space. 

In  the  above,  R™  and  RTE  are  generalized  reflection  coefficients  for 
TM  and  TE  waves,  respectively.  They  can  be  expresed  as 


R 


TE.TM 


RTE,TM  +  RTE,TMc2,klld 
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For  a  dielectrical  substrate  backed  by  a  PEC  ground  plane, 


(8) 


rjTE  _  Ml  kz  —  MfcU  _  kz  —  k\; 
1  Mi  kz  +  pkiz  kj  +  fcji 


RTE 

nX7 


-1, 


(9) 
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TM  __ 
01  “ 


f ikx  -  ikit 
eikt  +  ekit  ’ 


Therefore,  we  obtain 


R 


TM 

12 


=  1. 


R™  = 


icikx  cos  kizd  —  ek  j.  sinking 
2«i kz  cos  kizd  +  efc i,  sin  kltd' 


(10) 


(11) 


where 


rTE  _  ikx  sin  kltd  +  kix  cos  kixd 
ik x  sin  kixd—  klx  cos  klxd' 

kx  =  (k2-ki)K  klx=--(ki-ki)li, 

k2  =  k]  =  u>Vici  =  u.'2/iei, 


(12) 


(13) 


and  d  is  the  substrate  thickness.  Equations  (1)  and  (2)  are  the  consequence 
of  imposing  the  mixed  boundary  conditions  on  the  conductive  patches  at  the 
plane  2=0,  namely, 


e„(p)  =  0,  p  <  a,  and  a<p<b  (2  =  0),  (14) 

*n(p)  =  0,  outside  the  patches  (2  =  0)  (15) 

where  e„  and  ic„  are  given  by  equations  (1)  and  (2),  respectively. 


III.  Selection  of  Basis  Funtions 

In  equations  (1)  and  (2),  the  current  distribution  Kn(p )  or  its  VHT 
K n{kp)  is  unknown.  We  can  solve  for  the  unknowns  using  Galerkin’s  method. 
To  do  this,  we  expand  the  unknown  current  in  terms  of  a  basis  set.  An  orthog¬ 
onal  and  complete  set  of  eigenfunctions  derived  from  the  current  distributions 
of  the  modes  of  a  magnetic-wall  cavity  is  used  to  expand  the  currents  on  the 
patches. 

The  current  on  the  annular  ring  corresponding  to  the  field  of  a  TMnm 
mode  of  a  magnetic  wall  cavity  is  [7] 


*n(p)  = 


Am^nm{p)  a<  p  <b 
0  elsewhere 


(16) 


where 


and 


,  (  \  —  *l,n{0nmP/a-) 

^  ~  l/&Vn(/W/a)J  ’ 


— )  —  >7n(/?„m  — )Arn(/?ram)  ~  Jn(0nm)^  n{0ntn  ) 
a  a  (i 


(17) 


594 


and 


fini^ripP/ &)  —  Jni&npP/ &)  ^n(^np)  ^n(^np)^n(®«  ?p/~)  (38) 

The  eigenvalue  anp  is  determined  by  <j>n(anpc )  =  0.  The  VHT  of  «„(/>)  in 
equation  (26)  can  be  derived  to  be: 


where 


and 


K  n(fcp)  =  Bp<f>np(kp) 


<f>np{kp) 


0 

.  kl-(onr/a^ZnP^p)  . 


Znp(kp)  —  ^^>n(o!npC)  Jn(t:p6)  <}>n(o(nps)Jn(kpa). 


(29) 


(30) 

(31) 


For  the  circular  disk,  the  current  corresponding  to  the  field  of  the  TEn, 
mode  is  [8] 


,  ,  f  Dqrj nq(p)  p  <  «i 

"  ^  0  elsewhere 


where 


l 


Vnq{p)  = 


f^Jn(XnqP/a  l) 
L  J'niXnqp/a  0 


(32) 


(33) 


The  eigenvalue  \nq  is  determined  by  Jn{Xnq )  =  0,  and  is  the  g-th  root  of  the 
equation.  The  VHT  of  K„(p)  in  equation  (32)  can  be  written  as 


hn(l'p)  —  DqlJnq(k  p) 


(34) 


i»f  ore 
*t  livi 


_  r  o 


(35) 


A  complete  set  of  vector  basis  functions  approaching  the  current  distri¬ 
bution  on  the  antenna  patches  then  can  be  written  as 


*n{p)  =  < 


M  P 

^  ^  ^m^nmlp)  4"  ,  B p<f>nP(f')  O.  ^  p  <  b 

p=l 


m=r- 1 


l  Q 

C£nl{p)  +  ^2  DqVnq{p) 

1=1  q=l 

0 


P  <  a  i 
elsewhere 


(36) 
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and  the  VHT  of  (p)  in  equation  (36)  is 


M  P  L  Q 

Kn(fcp)  =  ^  Aml}>nm{kp)  +  ^  Bp4>np(kp)  +  C£nt{kp)  +  nq{kp). 

m=nl  p=  1  /sal  y=l 

(37) 

Substituting  equation  (3V)  into  (1),  we  obtain 
M  00 

f>«  /  MM*n(*v>)  •  G(kp)-^nm{kp) 

P  00 

E  B>Jk  pdkpUn(kpp)G{kp)  ■  ^np(kp) 

kpdkpHn(kpP)-G(kp)-(nl(kp) 


+ 


i  °° 

+X>/ 

j=i  ' 

«  7  _  _ 

~h  ^  Dq  J  kpdkpHn(kpp)  ■  G(kp)  ■  f{nq{kP)  ~  0. 


(38) 


It  can  be  shown  from  the  derivation  of  K ri{kp)  that  equation  (2)  is  satisfied 
automatically.  Therefore,  we  need  only  to  solve  equation  (38)  to  obtain  the 
natural  resonance  frequen.  ies  of  the  microstrip  antennas. 


IV.  Characteristic  Equation  for  Eigenmodes 

Equation  (38)  can  be  converted  to  a  matrix  equation  by  applying  the 
testing  functions  in  Galerkin’s  method.  Multiplying  equation  (38)  by  p0|„(p), 
p&njip),  pCfc(p)  and  pf/‘n,(p),  respectively,  and  integrating  with  respect  to  p 

from  0  to  oo  for  i  =  1, . ,M,  j  —  1, . ,P,  k  =  1, . ,L,  t  =  I, . ,Q , 

and  using  Parseval’s  Theorem  for  the  vector  Hankel  transform  [4],  equation 
(38)  became  a  system  of  M  4-  P  -f  L  +  <3  linear  algebraic  equations  and  can 
be  written  in  the  following  matrix  form 


■  a'*"7  " 

'  A' 

A**  Aw  A**  A4’ 

B 

AS‘  AH  A*  A<’ 

C 

_,v  T, 

.D. 

where 

OO 

—  J  dkpk^m(kp)-G(kp).f)n}(kPl 

o 


(39) 


(40) 
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where  a  and  (5  can  be  either  t/>,  <j>,  f  or  tj.  Due  t.o  reciprocity,  the  matrix 
in  equation  (39)  is  symmetric  or  A“  =  .  In  the  above,  the  superscript 

t  implies  transpose.  In  order  for  the  system  of  equations  (39)  to  have  a 
nontrivial  solution,  its  determinant  must  vanish.  Denoting  the  determinant 
in  equation  (39)  by  D(u?),  we  can  conclude  easily  that  the  natural  resonance 
frequencies  u u  of  the  microstrip  antenna  must  satisfy  the  following  equation 

D(tni)  ■=  0.  (40a) 


In  general,  the  roots  of  equation  (40)  are  complex  with  a  negative  imag¬ 
inary  part,  indicating  damping  due  to  radiation  loss. 


V.  Input  Impedance  of  the  Probe  Excited  Antenna 

In  order  to  find  the  input  impedance  of  the  ARL  circular-disk  microstrip 
antenna  excited  by  a  probe  as  shown  in  Figure  2,  we  have  to  consider  the 
excitation  probl  n  [1] .  In  this  case,  the  electric  field  on  the  patches  is  not 
only  due  to  the  current  on  the  patch,  but  also  due  to  the  current  on  the 
probe.  To  determine  the  transverse  electric  field  on  the  patches  due  to  a 
vertical  probe  wuh  a  finite  radius  embedded  in  the  aielecL  lie  suusLate,  we 
have  to  find  the  field  in  the  substrate  due  to  the  probe  first.  Using  the  dyadic 
Green’s  function  formalism,  we  can  express  the  z-component  of  the  electric 
field  in  the  substrate  due  to  the  probe  as  (1,  5] 


<  =  - 
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^  Jo{kpp)J0[kh,R) 


d — l  d+l 

> r—  z  < x — 


(41) 


where  R,  l  and  I  are  radius,  length,  and  assumed  uniform  current  of  the 
probe  in  the  2-direction,  respectively,  as  shown  in  Figure  2.  The  function 
sinc(i)  =  sin  (i)/x, 

Using  the  transmission  and  reflection  properties  of  the  electromagnetic 
waves  in  a  stratified  medium  [Gj,  we  can  find  the  field  in  the  upper-half  space 
from  equation  (41).  For  a  probe  location  at  p  —  xb\  using  the  addition 
theorem  we  can  express  the  field  in  terms  of  a  series  of  Dessel  functions. 
Then,  the  transverse  components  of  the  electric  field  in  the  upper-half  space 
due  to  the  probe  can  be  derived  from  Ez  yielding 


e£„  cos  n<p 

l^nj 

sin _ 

3yy 


(42) 


and 


where 
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where  (5„m  is  the  Kronecker  delta  function  and  X™  =  14-  Rjo4  ■ 


(45) 


When  the  ARL  circular-disk  antenna  is  being  excited  by  a  probe,  the 
current  distribution  on  the  disk  in  the  vicinity  of  the  probe  has  a  singularity 
due  to  the  singular  nature  of  the  field  around  the  probe,  which  results  in 
slow  convergence.  Therefore,  the  current  J§{p)  of  a  lossy  magnetic  wall 
cavity  under  a  probe  excitation  is  added  to  the  regular  current  distribution 
j£(p)  on  the  patch  to  accelerate  the  convergence  [5]  and  the  VHT  of  J§(p) 

i  p  i  •  _  .i j  r _ _  -  f  i  cl 

Ctin  ut;  louiiu  in  ciuscu  luim  oo  [i,uj 
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Jn(fel.y) 
J'n(k  L«) 
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ki -i*  dn(kpa) 
_  at.Mhisl] 

ki,k,a  J 


}■ 


(46) 


where  kl3  is  the  wavenumber  of  the  dielectric  medium  in  the  magnetic  wall 
cavity  and  can  be  written  as 

*.,  =  (l  +  ia)—  (47) 

a  i 

where  q  is  a  positive  loss  factor.  The  expression  in  (4C)  consists  of  both 
the  singular  term  and  the  reflected  wave  term.  Since  we  axe  only  interested 
in  the  singular  term  in  (46),  we  can  choose  a  to  increase  the  loss  such  that 
the  reflected  wave  term  is  unimportant.  Considering  the  excitation  condition 
arid  using  the  singularity  subtraction,  we  can  rewrite  equation  (1)  as 


J  dkpkplin(kpp)  ■  G (*„)  •  [K *{kp)  +  K* 

o 


d^v.H„(V)S„(^)  =0 


(48a) 
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3yy 


or 


oo 

J  dkpkpHn(kpP)  •  G (fcp)  •  K«(fcp) 

0 

oo 

-  [ dkpkpUn(kpP ) .  [Sn(fc„)  +  G(fcp) .  K*(fe„)]. 

*  (48b) 

Using  Galerkin’s  method,  performing  operations  similar  to  those  given  in 
Section  IV,  we  obtain  a  system  of  linear  algebraic  equations.  Solving  the 
matrix  equation,  the  current  distribution  on  the  patches  and  then  the  z- 
component  of  the  total  field  E  in  the  substrate  can  be  determined.  Therefore, 
the  input  impedance  of  the  antenna  under  a  probe  excitation  can  be  derived 
by  using  the  induced  EMF  method  [9] 

Zin  =  -  j  E  •  3dV/I2  (49) 

v 

The  above  is  a  variational  expression  for  Z,n.  The  input  impedance  expres¬ 
sion  in  terms  of  the  current  distribution  both  on  the  patches  and  on  the 
probe  can  be  derived  from  equation  (49)  as 
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where  the  subscript  nl  denotes  the  first  elements  of  vectors  Kjj,  K^f.  The 
vector  K^,  which  corresponds  to  the  regular  part  of  the  current  distribu¬ 
tion  on  the  disk,  is  obtained  by  solving  the  integral  equation  (48b)  using 
Galerkin’s  method.  It  is  given  by 


At  i 

^  Anrnil>nm{kp)  +  Cn£nl{kp).  (51) 

m=  1  1=1 
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VI.  Results  and  Discussions 


In  this  paper  the  vector  integral  equations  are  solved  using  Galerkin;s 
method.  As  in  other  problems  solved  with  this  method,  the  convergence  of 
the  solution  has  to  be  considered.  Furthermore,  to  enhance  the  accuracy 
of  the  numerical  calculation,  the  integral  paths  are  deformed  to  avoid  the 
singularities  near  the  real  axis.  The  oscillatory  properly  of  the  integrand 
and  the  convergence  of  the  infinite  integral  have  to  be  evaluated  a  priority 
to  select  a  suitable  numerical  integral  method  and  appropriate  truncation 
points  for  the  infinite  integrals.  The  singularity  subtraction  method  has  been 
used  to  accelerate  the  convergence  of  the  integral  with  a  singular  current 
distribution.  The  slowly  decaying  parts  of  the  integral  in  equation  (50)  have 
been  subtracted  out  and  evaluated  in  close  form. 

The  numerical  analysis  for  natural  resonance  frequencies  in  this  paper 
are  the  same  as  those  in  Figure  3  of  [7]  when  L  —  0  and  Q  —  0  (for  the  case 
of  an  annular-ring  antenna)  or  as  those  in  Figure  3  of  [8]  when  M  =  0  and 
P  —  0  (for  the  case  of  a  circular-disk  antenna).  This  is  because  the  above 
are  the  special  cases  of  ARL  circular-disk  microstrip  antennas. 

As  shown  in  the  theoretical  analysis  and  the  experiment,  ARL  circu¬ 
lar  disk  microstrin  antennas  have  dual  resonance  frequencies,  which  can  be 
used  to  broaden  the  bandwidth  of  the  antenna.  In  general,  one  resonance  is 
higher  than  the  natural  resonance  of  the  annular  ring  antenna,  and  the  other 
resonance  is  lower  than  the  natural  resonance  of  the  circular-disk  antenna. 
The  difference  depends  on  the  mutual  coupling  (odd  and  even)  between  the 
two  patches;  therefore,  it  depends  on  the  geometrical  configuration  of  the  mi¬ 
crostrip  antenna  and  the  dielectric  constant  of  the  substrate.  Figure  3  shows 
the  dual  resonance  frequencies  as  the  function  of  the  substrate  thickness.  An 
appropriate  excitation  condition  is  important  to  obtain  the  dual  resonance 
frequences. 

Figure  4  shows  the  resonance  frequencies  and  the  bandwidth  of  two  ARL 
circular-disk  microstrip  antennas  with  different  patch  sizes  as  a  function  of 
substrate  thickness.  It  can  be  shown  that  the  resonance  frequencies  decrease 
and  the  bandwidths  increase  as  the  thickness  of  the  substrate  increases.  The 
numerical  results  are  in  good  agreement  with  experiment  data  even  when  the 
substrate  thickness  increases  to  0.1  wavelength. 

Figures  5-7  show  the  input  impedance  characteristics  of  the  ARL  cir¬ 
cular-disk  microstrip  antenna  with  various  electrical  thicknesses  as  a  function 
of  frequencies.  In  Figure  5  the  electrical  thickness  of  the  antenna  substrate 
is  d  =  0.054A,;  M,  P,  L,  Q ,  the  number  of  the  basis  functions  for  current 
distribution  on  the  patches,  are  2, 1,2,1,  respectively.  The  agreement  between 
theoretical  results  and  experimental  data  is  excellent. 
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Figure  6  shows  the  input  impedance  of  the  antenna  with  electrical  thick¬ 
ness  d  =  0.07A<,  The  numbers  of  the  basis  functions  for  the  current  on  the 
patches  are  the  same  as  in  Figure  5.  The  results  of  the  induced  EMF  method 
still  agree  well  with  the  experiment  data.  The  impedance  locus  shift  between 
the  induced  EMF  method  and  experiment  data  is  the  effect  of  the  edge  ex¬ 
tension  caused  by  the  fringing  field  at  the  edge  of  the  patches.  This  fringing 
field  effect  is  caused  by  singular  charge  distribution  near  the  edge;  hence, 
more  basis  functions  are  required  to  approximate  this  fringing  field  effect. 
Therefore,  the  resonances  obtained  by  numerical  calculation  with  an  insuf¬ 
ficient  number  of  basis  functions  tend  to  be  higher  than  that  obtained  from 
the  experiment.  The  thicker  the  substrate,  the  more  significant  the  fringing 
field  effect,  therefore,  the  greater  the  shift  between  the  locus  obtained  from 
the  theory  and  experiment  for  a  given  number  of  basis  functions. 

Figure  7  shows  the  input  impedance  of  the  antenna  with  electrical  thick¬ 
ness  d  =  O.lA*.  The  number  of  the  basis  functions  of  the  TM  mode  for  current 
distribution  on  the  circular  disk  is  increased  to  8  (L  =  8)  to  account  for  the 
fringing  field  effect.  A  very  good  agreement  between  the  numerical  and  the 
experimental  results  has  been  achieved.  Increasing  the  number  of  basis  func¬ 
tions  for  the  circular  disk  rather  than  for  the  annular  ring  is  because  the 
dominant  resonance  is  coming  from  the  disk  rather  than  the  ring. 

Figure  8  shows  the  input  impedance  locus  for  L— 2  and  L=S  when  the 
electrical  thickness  d  —  0.1  A,.  The  edge  extension  caused  by  the  fringing 
field  at  the  edge  of  the  patch  cannot  be  neglected,  and  more  basis  functions 
of  the  TM  modes  have  to  be  chosen  to  account  for  the  fringing  field  effect 
when  the  substrate  is  electrically  thick. 

It  can  be  shown  from  the  above  discussion  that  the  theory  in  this  paper 
can  be  used  to  analyse  the  characteristics  of  microstrip  antennas,  especially 
the  antennas  with  mutual  coupling  between  patches,  which  are  difficult  to 
analyse  with  some  popular  methods  such  as  the  cavity  model  method.  Fur¬ 
thermore,  this  theory,  when  applied  to  the  analysis  of  the  microstrip  anten¬ 
nas  with  an  electrically  thick  substrate,  gives  results  of  resonances  and  input 
impedances  in  good  agreement  with  experimental  data  if  we  use  more  basis 
functions  of  the  TM  modes  tc  account  for  the  fringing  field  effect.  How¬ 
ever,  much  more  computer  resources  are  required  in  such  a  case.  It  is  also 
shown  from  the  numerical  analysis  that  the  effect  of  the  probe  on  the  res¬ 
onances  is  much  less  than  the  effect  or.  the  fringing  field  of  the  patch  when 
the  microstrip  antenna  is  electrically  thick.  We  have  not  used  any  artifical 
correction  on  the  theory  for  the  good  agreement  with  measurement  results. 
The  uniform-probe-current  assumption  does  not  affect  the  input  impedance 
appreciably  because  the  input  impedance  is  dominated  by  the  disk  current, 
even  for  antennas  with  electrically  thick  substrates. 


References 

[1]  W.  C.  Chew  and  J.  A.  Kong,  “Analysis  of  circular  micro-  trip  disk  an¬ 
tenna  with  thick  dielectric  substrate”,  IEEE  Tran?.  Artcrnas  Propagat., 
v.  AP-29,  pp.  68-76,  Jan.  1981. 

[2]  M.  Davidovitz  and  Y.  T.  Lo,  “Input  impedance  of  a  probe-fed  circu¬ 
lar  microstrip  antenna  with  thick  substrate,”  IEEE  Trans.  Antennas 
Propagat .,  v.  AP-34,  pp.  905-911, July  1986. 

[3]  S.  M.  Wright  and  Y.  T.  Lo,  “Efficient  analysis  of  infinite  microstrip 
antenna  arrays  on  electrically  thick  substrates”,  RADC-TR-86-8  Final 
Technical  Report,  Rome  Air  Development  Center,  University  of  Illinois, 
March  1986.  ADA169011 

[4]  W.  C.  Chew  and  T.  M.  Habashy,  “The  use  of  vector  transform  in  solv¬ 
ing  some  electromagnetic  scattering  problems”,  IEEE  Trans.  Antennas 
Propagat .,  v.  AP-34,  pp.  871-879,  July  1986. 

[5]  W.  C.  Chew,  Mixed  Boundary  Value  Problems  in  Microstrip  and  Geo¬ 
physical  Probing  Applications,  Ph.  D.  dissertation,  Massachusetts  Inst. 
Technol.,  Cambridge,  1980. 

[6]  J.  A.  Kong,  Theory  of  Electromagnetic  Waves.  New  York:  Wiley,  1975. 

[7]  S.  M.  Ali,  W.  C.  Chew  and  J.  A.  Kong,  “Vector  Hankel  transform  anal¬ 
ysis  of  annular-ring  microstrip  antenna,”  IEEE  Trans.  Antennas  Prop¬ 
agat.,  v.  AP-30,  pp.  637-644,  July  1982. 

[8]  W.  C.  Chew  and  J.  A.  Kong,  “Resonance  of  nonaxial  symmetric  modes 
in  circular  microstrip  disk  antenna”,  J.  Math.  Phys.,  v.  21,  n.  10,  pp. 
2590-2598,  Oct.  1980. 

[9]  E.  C.  Jordan  and  K.  G.  Balmain,  Electromagnetic  Waves  and  Radiating 
Systems.  Prentic  Hall:  New  Jersey,  1968. 


13 


4U3 


Fig.l  Geometrical  configuration  of 
annular-ring' loaded  circular  disk 
antenna 
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Fig. 2  ARL  circular-disk  microstrip 
antenna  excited  by  a  probe 
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Fig. 5.  Comparison  of  experiment  and  theory  for  input 
impedance,  where  d  =  0.054  wavelength.  full  symbols  —  theory, 
open  symbols  —  experiment. 
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DESIGN  AND  ECONOMIC  CONSIDERATIONS  FOR  BROADBAND  PHASED-ARRAY  ELEMENTS 

George  J.  Monser 


1.0  INTRODUCTION 

Today's  technology  and  economy  demand  that  the  arrays  function  over 
multi-octaves  and  that  the  cost  per  element  be  small.  This  talk  addresses 
both  needs.  Three  stages  of  linear-array  development  are  presented.  Stage  1, 
where  the  consideration  was  for  a  linearly-polarized  array.  Stage  2,  where 
the  consideration  was  to  extend  the  array  to  include  an  orthogonal  polariza¬ 
tion.  The  final  stage  describes  areas  of  consideration  for  maximizing  the 
array  bandwidth. 

2 . 0  STAGE  1  DESIGN  CONSIDERATIONS 

In  forming  a  linear  array,  as  depicted  in  Figure  1,  two  dimensions  are 
important.  The  element  dimension  d_,  which  is  aligned  with  the  array,  is  fixed 
by  the  highest  frequency  of  intended  operation.  If  d_  exceeds  1/2  wavelength, 
secondary  lobes  will  appear  when  the  beam  is  scanned  beyond  about  30°  from 
boresight.  For  a  variety  of  reasons,  system  requirements  generally  Impose  a 
limit  on  the  formation  of  secondary  lobes  so  that  the  dimension  is  approxi¬ 
mately  ■ /2  where  .A  is  the  wavelength  corresponding  to  the  uppermost  frequency. 
The  dimension  H,  perpendicular  to  the  array,  represents  a  free  dimension  in 
linear  array  design  limited  only  by  the  minimum  f ield-of-view  (FOV)  orthogonal 
to  the  array.  The  FOV  in  degrees  is  given  approximately  by  51  (A/H),  where 
A,  is  the  free-space  wavelength  and  H  is  the  height  both  in  consistent  measure¬ 
ment  units  (i.e.,  inches,  CMs ,  etc.). 

Once  11  is  chosen,  the  3pace  per  array  element  for  interior  element 
design  is  set  for  the  depth  or  length  of  the  array  element  are  generally 
limited  by  the  installation  space.  Where  space  is  a  premium,  such  as  in 
tactical  aircraft,  sufficient  depth  for  an  optimum  design  may  not  exist. 

Figure  2  shows  an  array  element  cell  of  a  phased-array  with  dimensions  H 
and  ji.  Both  11  and  _d  are  determined  by  the  ’  '^hest  frequency  of  operation.  For 
example,  if  the  highest  frequency  corresponds  to  0.65  inches  in  wavelength  and 


Element:  Dimensions 


a  FOV  of  30®  is  required,  then  from  equation  (1),  H  *  1.105  inches,  and  from 
the  (1/2)X  requirement,  d  *  0.325  inches.  Thus,  the  array  element  dimensions 
for  this  application  are  fixed  (1.105  x  0.325  inches). 

If  a  bit  of  the  grating  lobe  is  allowed  to  exist,  a  slightly  larger  d_ 
can  be  used.  If  the  FOV  is  re-defined,  from  say  the  half-power  limits  to  the 
1/3  power  limits,  then  a  slightly  larger  is  permitted.  In  any  case,  both 
dimensions  are  established. 

3 . 0  STAGE  2  DESIGN  CONSIDERATIONS 

With  the  array-element  real  estate  established,  it  Is  necessary  to 
look  for  array  elements  chat  can  be  fit  within  these  constraints  of  Figure  2, 
and  which  are  economical  to  produce.  A  horn  within  which  ridges  can  be  placed 
to  yield  a  3/1  bandwidth  was  selected  as  one  candidate  element.  Mechanical 
tolerance  and  performance  studies  were  conducted  to  support  the  use  of  cast 
elements  for  the  final  design.  For  the  other  candidate  element,  it  was 
visualized  char  an  axray  of  printed  notches,  with  one  board  per  array  ele¬ 
ment,  offered  an  economical  and  broadband  solution.  In  addition,  the  simplest 
solution  for  integrating  the  printed  notches  In  the  horn  array  would  consist 
of  replacing  a  section  of  each  horn  wall  with  the  printed  notch  board.  The 
notches  radiate  primarily  from  the  step  discontinuity  near  the  aperture  plane 
so  that  they  would  function  effectively  and  so  would  the  horns  for  the  walls 
necessary  for  propagating  the  TEj^q  mode  would  be  nearly  intact  except  for  the 
small  notches. 

Photographs  of  these  elements  are  shown  in  Figure  3.  Since  gain  versus 
frequency  is  most  important,  each  element  was  first  tested  as  an  isolated  ele¬ 
ment,  Results  are  given  in  Figure  4.  In  this  figure,  the  results  are  plotted 
versus  the  critical  design  dimension  expressed  as  a  fraction  of  a  wavelength. 

In  both  cases,  the  critical  dimension  is  for  the  notches  aligned  with  the  a_ 
dimension  of  the  horn  (see  Figure  2).  To  achieve  the  extended  band  for  the 
horn  element  it  is  necessary  to  place  ridges,  as  shown,  inside  the  array  ele¬ 
ment.  The  ridges-  Lower  the  cut-off  frequency  of  the  horn.  Without  ridges, 
the  cut-off  frequency  corresponds  to  a  wavelength  equal  to  (2  x  a).  With 
ridges,  the  cut-off  wavelength  can  be  increased  significantly,  to  at  least 
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Iiotographs  of  Elements 


GAIN  RELATIVE  TO  GAIN  STANDARD  (dB) 


(3  x  ) ,  reaulting  in  •  lower  useable  frequency  as  frequency  and  wavelength 
are  Inversely  related.  A  further  point  should  be  made  and  that  is,  while 
the  H.  or  £  dimension  is  fixed  at  the  aperture,  by  the  minimum  FOV,  inside 
the  horn  the  _a  dimension  is  selected  to  achieve  the  low  frequency  limit  of 
operation  (i.e.,  4  GHz)  and  at  the  same  time  avoiding  higher  order  waveguide 
modes.  Should  this  occur  a  lowering  of  the  uppermost-  frequency  would  result. 

Heving  selected  the  elements,  an  array  was  constructed  consisting  of 
both  types  of  elements  interlaced  so  that  two  polarizations  could  be  trans¬ 
mitted  £r  received  by  the  array.  One  polarization  is  intercepted  or  formed 
by  the  horn  elements,  parallel  to  the  array,  and  the  other  polarization,  per¬ 
pendicular  to  the  array  is  intercepted  or  formed  by  the  notch  elements.  A 
photograph  of  this  array  is  shown  in  Figure  5_. 

4.0  PERFORMANCE  TESTS 

Having  constructed  the  array,  array  element  VSWR  was  first  measured, 
as  it  determines  element  efficiency.  Here  it  is  important  to  note  that  the 
VSWR  of  concern  is  that  realized  when  all  array  elements  ate  simultaneously 
driven.  This  is  referred  to  as  the  active  VSWR  as  contrasted  to  the  passive 
VSWR  measured  on  a  single  element  in  the  array  with  all  other  elements  termi¬ 
nated  or  loaded.  In  this  design  it  was  found  that  at  the  lower  frequencies 
the  passive  VSWR  for  the  horn  elements  predicted  efficiencies  on  the  order 
of  50  percent.  However,  from  active  VSWR  values  closer  to  80  percent  result. 
Active  match  can  be  found  by  measuring  the  passive  reflection  coefficient  and 
the  coupling  coefficients  from  ail  elements  to  element  under  teat,  with  all 
quantities  expressed  as  complex  values  or  phasors.  The  inset  box  (see 
Figure  6)  shows  how  active  match  is  found.  There,  it  may  be  observed  that 
for  the  horn  element  the  passive  reflection  partially  cancels  the  coupled 
energy  resulting  in  a  low  active  reflection  coefficient  and  good  active  VSWR. 
However,  for  the  notch  elements  the  vectors  do  not  cancel  to  the  same  extent. 
For  this  element,  the  lowest  passive  reflection  is  achieved  by  design  and  the 
resulting  coupling  and  active  VSWR  tolerated. 
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Figure  6a.  Phasoc  Relationships  (Horns) 
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CASE  NO.  2:  PHASORS  FOR  S/L 
NO.  5,  0°  SCAN,  F=6.5  GHz 


Once  chess  VSWR  determinations  were  complete,  gains  for  all  elements 
were  checked  over  the  full  frequency  band  to  assure  continuous  coverage. 
Typical  swept-gain  plots  are  given  in  Figure  7.  There,  it  can  be  seen  mat 

none  of  the  elements  had  drop-outs  in  gain  over  the  full  band  Prom  5  to  m  wnz 

(i.e.,  continuous  coverage). 

Pattern  coverage  for  each  type  of  element  is  given  in  Figures  8  and  9. 

In  Figure  8,  coverage  in  the  array-plane  at  mid-band  is  shown  to  exceed  100 

between  half-power  points  (i.e.,  —3  dB  points  on  the  patterns).  Thus,  when 
all  elements  of  the  array  are  driven  with  phasing  set  to  steer  the  beam,  good 
scan  characteristics  should  result  over  a  100°  sector  (i.e.,  _^50n  from  bore- 
sight).  Figure  9  shows  the  coverage  orthogonal  to  the  array  exceed  30°  for 
half  power  limits  (i.e.,  -3  dB  points  on  the  patterns). 

5.0  COST  ANALYSIS 

Production  costs  for  each  type  of  array  element  were  estimated. 
Excluding  one-time  tooling  costs  for  each  element,  it  was  determined  that 
each  array  element  (both  polarizations)  could  be  manufactured  for  well 
under  $100.  For  example,  the  cost  of  cast  horn  elements  was  $10.  Adding 
$10  per  element  for  "touch  labor"  in  assembly  and  $25  for  production  test 
yields  a  total  cost  under  $50. 

6.0  CONCLUDING  REMARKS 

This  talk  has  described  the  design  of  an  economically-producible 
dual-polarized  array -capable  of  4/1  bandwidth  for  both  polarizations  and 
6/1  bandwidth  for  the  printed-notch  array. 

Recent  investigative  efforts  have  indicated  that  the  bandwidth  for 
the  horn  elements  can  possibly  be  extended  to  cover  a  6/1  band. 
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Figure  7b.  Swept  Gain  Full  Band 

(Vertical  and  Horizontal  Polarizations) 


- H  (HORN) 

(GAINS  RELATIVE) 
4/7/81 


Figure  8.  Coverage  in  Array  Plane  at  Mid-Band 
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A  LARGE  MULTIBEAM  LENS  ANTENNA  FOR  EHF  SATCOM 
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Electromagnetic  Sciences,  Inc. 

Norcross,  Georgia  30092 

ABSTRACT 

Development  of  very  large,  satellite  based,  multiple  beam  antennas  (MBAs)  using  ex¬ 
tremely  high  frequencies  (EHF)  for  satellite  communications  (SATCOM)  is  the  subject  of 
this  paper.  Very  large  MBAs  having  several  hundred  beams  is  a  significant  Increase  in  the 
complexity  of  an  MBA  by  today's  standards.  The  main  issue  to  be  resolved  is  whether  the 
beamforming  network  (BFN)  required  for  such  an  MBA  can  be  implemented  in  a  practical 
manner.  The  paper  is  extracted  from  a  more  complete  report  (Montgomery,  Runyon,  and 
Fuller,  1988]  recently  published  and  presents  strong  support  for  the  feasibility  of  these 
large  MBAs. 

Two  primary  approaches  to  realizing  the  BFN  are  examined:  (1)  a  constrained  BFN  using 
interleaved  switch  networks  with  an  adaptive  combining  network,  and  (2)  a  space  fed 
array  consisting  of  modules  capable  of  precision  phase  and  amplitude  control. 
Consideration  is  given  to  the  addition  of  solid  state  low  noise  amplifiers  for  uplink  an¬ 
tennas,  however,  nulling  requirements  place  stringent  requirements  on  these  components. 
Both  approaches  are  compared  with  regard  to:  gain  and  noise  temperature,  adaptive  null¬ 
ing,  power  requirements,  weight  and  size,  complexity  and  cost,  reliability,  and  robustness. 

Common  to  either  MBA  approach  is  the  use  of  a  large,  high  performance,  zoned  dielectric 
lens  and  the  use  of  feed  dusters  to  obtain  improved  scanning  and  sideiobe  performance. 

h  INTRODUCTION 

The  use  of  EHF  for  SATCOM  is  expected  to  grow  substantially  over  the  next  decade. 
EHF  frequencies  allow  more  user  channels  while  accommodating  wide  band  users.  EHF 
frequencies  also  permit  high  efficiency  antennas  with  narrow  beamwidths,  high  gain,  and 
low  sidelobes  with  minimal  size  and  weight  penalties  —  ideal  for  satellite  based  systems. 
These  antenna  characteristics  are  well  suited  for  minimizing  the  adverse  effects  of  spa¬ 
tially  distributed  jammers  both  with  and  without  the  use  of  adaptive  nulling  networks. 

Electronic  scanning  of  an  MBA  is  accomplished  by  the  electronic  control  of  th©  BFN.  The 
critical  EHF  components  such  as  switches  and  phase  shifters  which  comprise  the  BFN 
typically  utilize  ferrite  technology  which  is  capable  of  achieving  low  insertion  losses  with 
low  switching  powers  while  maintaining  the  reliability  required  for  satellite  applications. 
Typically,  the  BFN  uses  power  from  one  or  mere  antenna  feed  ports,  by  controlling  switch 
matrices  and  variable  power  combiners,  to  move  the  active  port  centroid  in  the  focal  region 
of  either  a  reflector  or  lens  to  provide  an  electronically  agile  beam. 

This  paper  considers  antennas  with  component  beams  having  a  0.5*  half  power 
beamwidth.  Depending  on  the  beam  spacings,  this  beamwidth  implies  approximately  700- 
900  beams  within  the  earth  fieid-of-view  (FOV).  By  utilizing  a  high  performance  dielectric 
lens  with  component  beam  clusters,  the  MBA  has  considerable  flexibility  in  controlling  the 
secondary  pattern  characteristics  over  the  entire  earth  FOV  while  also  permitting  high 
performance  adaptive  nulling. 
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2.  LEAM  FORMING  METHODS 


Two  techniques  for  the  realization  of  the  large  MBA  are  examined:  (11  a  traditional, 
constrained  BFN  using  seven  interlaced  switch  trees  (refer  Figure  11  with  a  seven  port 
variable  power  combiner  (VPC)  and  phasing  network;  and  (2|  a  space  fed  array  comprised 
of  precision  phase  and  amplitude  control  modules  at  each  port  -  and  denoted  as  a  beam 
port  array  (refer  Figure  2).  The  interlaced  switch  tree  is  only  used  with  the  constrained 
BFN,  The  beam  port  array  achieves  the  same  excitation  flexibility  with  a  fundamentally 
different  technique. 

Beam  Geometry 

A  hexagonal  beam  geometry  is  assumed  for  both  techniques.  Figure  3  illustrates  the 
hexagonal  geometry  and  defines  the  two  beam  crossover  and  the  triple  crossovers  (i.e., 
three  beam  crossovers).  A  fundamental  design  parameter  for  either  method  is  the  number 
of  beams  required.  Table  1  illustrates  the  number  of  beams  within  the  circular  earth  FOV 


TABLE  1.  Hexagonal  Beam  Parameters 


nber  of  Beams 

Beam  Spacing* 

7 

7.44782° 

19 

4.13131 

37 

2.97913 

61 

2.27156 

91 

1.86195 

127 

1.56149 

265 

1.06397 

685 

0.64764 

769 

0.60761 

859 

0.57291 

1  Assumes  earth  radius  =  8.6® 


and  their  spacing  when  the  earth  radius  passes  through  the  triple  crossover  of  the  hexago¬ 
nal  grid.  In  order  to  determine  the  beamwidths  of  the  component  beams,  the  concept  of 
minimum  gain  for  a  switched  beam  system  is  used.  The  minimum  gain  of  a  switched  beam 
system  occurs  at  the  triple  crossover  point  and  is  related  to  the  two  beam  crossover  angle 
as  shown  in  Figure  3.  The  minimum  gain  (relative  to  peak)  approximately  determines  the 
two  beam  crossover  (relative  to  peak)  and  the  beamspacing  to  beamwidth  ratio  by  using  a 
Gaussian  beam  model.  Table  2  illustrates  these  quantities  for  a  parameter  range  of  inter¬ 
est.  The  final  crossover  levels  are  chosen  by  maximizing  the  minimum  gain. 
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Figure  2.  Beam  Port  Array  Concept 
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TABLE  2.  Beam  Crossover  Parameters 


Two  Beam 

Tripie 

Beam  Spacing/ 

Crossover 

Crossover 

Beam  width 

-2.26  dB 

-3.00  dB 

0.866 

-3.00 

-4.00 

1.000 

-4.00 

-6.33 

1.166 

-6.00 

-6.87 

1.291 

-6.60 

-7,33 

1.354 

-6.00 

-8.00 

1.414 

A  proper  determination  cf  the  minimum  gain  requires  the  consideration  of  the  effects  of 
the  finite  feed  horn  spacing  on  the  spillover  loss  and  taper  efficiency.  Figure  4  illustrates 
the  minimum  gain  as  a  function  of  the  two  beam  crossover  level  for  the  case  of  685  and 
869  beams.  The  solid  and  dotted  curves  are  the  result  of  an  approximate  analysis  while  the 
discrete  data  represent  the  results  of  a  more  complete  lens  ray  tracing  analysis  using  a 
conical  horn  pattern  with  finite  wall  thickness.  The  reduced  gain  of  the  ray  tracing  solution 
in  comparison  to  the  approximate  result  is  a  consequence  of:  (11  the  nonsymmetric  conical 
horn  pattern,  (2)  a  finite  feed  hom  wall  thickness,  and  (3)  lens  ohmic  loss  and  reflection. 
The  minimum  gain  derived  with  either  analysis  is  maximum  when  the  crossover  level  is 
approximately  5.6  dB  below  the  peak.  This  value  is  lower  than  the  3  dB  which  is  com¬ 
monly  used.  However,  note  that  the  gain  in  the  vicinity  is  slowly  varying;  for  example,  a  4 
dB  crossover  results  in  approximately  a  0.4  dB  gain  reduction.  Furthermore,  there  are 
frequently  added  considerations  (e.g.  nulling  or  low  sidelobe  cluster  excitations)  which  can 
favor  the  use  of  higher  crossover  levels.  Nevertheless,  a  5.6  dB  two  beam  crossover  level 
is  used  for  both  beam  forming  methods  studied. 

Constrained  BFN 

Figure  6  illustrates  the  port  designation  pattern  of  the  seven  interlaced  switch  networks 
in  the  central  region  of  a  large  constrained  BFN.  Each  of  the  switch  networks  is  observed 
to  have  a  hexagonal  periodicity  which  is  related  to  the  fundamental  port-to-port  spacing. 
Mnr*  imnortjinHv.  this  cwtrrmt.rv  np.miit.s  nlanar  switch  networks  for  each  of  the  seven 

- r  - - J  »  -  a -  V  M -  *■  — 

interlaced  networks.  This  is  more  easily  seen  in  Figure  8  (also  Figure  1),  which  illustrates 
the  port  designations  for  an  868  port  constrained  BFN  and  the  planar  row  and  column 
connections  for  the  number  4  port,  The  row-to-row  spacing  of  the  planar  networks  is  found 
to  be:  At  =  0.32733  s,  where  s  is  the  port-to-port  spacing.  A  typical  row-to-row  spacing  is 
0.190  inch  and  accommodates  EHF  switch  heights  and  permits  the  use  of  standard  height 
waveguide.  The  port-to-port  spacing  (typically  1.638  inch)  within  each  planar  network 
easily  permits  the  use  of  added  switches  at  each  port  for  increased  isolation.  Figure  7 
illustrates  one  of  the  row  networks  for  the  869  port  feed.  Because  of  the  periodicity,  each 
network  is  similar  to  all  others  and  only  differ  by  which  ports  are  deleted.  For  the  869  port 
feed,  the  basic  building  block  network  is  a  18:1  binary  divider.  If  a  port  of  the  binary 
network  is  not  required  for  a  particular  row,  it  is  omitted.  It  is  assumed  that  the  number 
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Minimum  Gair  vs.  Crossover  Loss 


Figure  5.  Seven  Port  Constrained  BFN  Layout  Pattern 


Figure  0.  Row  end  Column  Switch  Tree  Connections  for  869  Port  MBA 
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Figure  7.  Typical  Row  Switch  Matrix  for  869  Port  MBA 


of  switches  in  any  path  is  a  constant  in  order  to  facilitate  port  phase  tracking  over  the 
frequency  and  temperature  ranges  of  interest.  This  requires  the  use  of  dummy  switches  as 
shown  in  Figure  7.  Figure  7  also  illustrates  that  the  number  of  unique  part  types  is  small 
even  with  a  discrete  implementation,  although,  the  quantities  of  each  type  may  be  large. 
Techniques  for  an  integrated  network  assembly  are  available  and  include:  flangeless  switch 
housings,  drop-in  switches,  or  an  integrated  housing  such  as  shown  in  Figure  8.  Indeed, 
the  fabrication  of  large  BFNs  depends  on  using  techniques  such  as  these  to  render  the 
method  practical. 

The  element  of  the  constrained  BFN  which  ultimately  permits  the  flexibility  afforded  by 
feed  clusters  is  the  variable  power  combiner  and  phasing  network  as  shown  in  Figure  9a. 
The  network  is  comprised  of  two  types  of  components:  VPCs  and  360°  phase  shifters.  The 
VPC  is  realized  by  utilizing  hybrids  and  phase  shifters  as  shown  in  Figure  9b.  The  phase 
shifters  for  the  VPC  need  only  have  a  differential  electrical  length  of  90®  (for  the  mixed 
hybrids)  or  180*  (for  like  hybrids).  The  requirements  of  the  network  are  normally  dictated 
by  the  adaptive  nulling  properties  and  would  indicate  phase  errors  should  be  less  1  ®  and 
amplitude  errors  less  than  0.1  dB.  These  type  of  errors  certainly  indicate  high 
performance  components. 

Beam  Port  Array 

The  beam  port  array  (BPA)  is  a  fundamentally  different  feed  from  the  constrained  feed. 
The  BPA  illustrated  in  Figure  2  is  space  fed  by  a  focused  phased  array  which  is  capable  of 
creating  a  localized  illumination  in  the  the  BPA  aperture.  The  centroid  of  the  illumination 
can  be  varied  and  hence  results  in  beam  scanning.  In  order  to  provide  the  same  cluster 
control  as  the  constrained  BFN,  the  BPA  uses  phase  and  amplitude  modules  (PAMs)  at 
each  element.  The  PAMs  are  realized  by  using  a  variable  power  combiner  (VPC)  with  an 
internal  termination  of  the  isolation  port.  The  phase  control  is  achieved  by  using  360® 
phase  shifters  in  contrast  tc  the  90*  or  180®  phase  shifters  normally  used  in  VPC  net¬ 
works.  The  BPA  offers  a  potential  advantage  for  eliminating  the  waveguide  switch  net¬ 
works  and  the  VPC  networks  of  the  constrained  BFN. 


Beam  port  array  concepts  have  been  examined  previously  by  Tang  and  Winter  (1973). 
Tang  and  Winter  used  a  focused  phased  array  in  conjunction  with  a  constrained  lens  in 
both  the  feed  region  and  the  radiating  aperture.  In  contrast  to  the  present  study,  Tang 
and  Winter  used  half  wavelength  element  spacings  for  both  the  focused  phased  array  and 
the  feed  constrained  ions;  furthermore,  Tang  and  Winter  utilized  nonplanar  surfaces  in 
order  to  maximize  the  illumination  efficiencies.  Thi9  previous  antenna  design  should  be 
contrasted  with  the  present  embodiment,  where  the  port-to-port  spacing  of  the  BPA  is 
approximately  two  wavelengths.  Furthermore,  the  goal  of  the  present  design  is  to  use 
planar  interfaces  on  either  side  of  the  BPA  with  compensating  cluster  amplitude  and 
phase  weights.  Additionally,  a  phase  and  amplitude  module  is  utilized  at  each  element 
instead  of  a  simple  constraining  transmission  line  connection.  These  nontrivial  differences 
are  intended  ultimately  to  reduce  the  system  cost. 


Figure  10  illustrates  a  single  phase  and  amplitude  module  utilizing  ferrite  components 
and  includes  the  feed  horns  and  a  circular  polarizer.  Figure  11  illustrates  a  37  port  subar¬ 
ray  of  these  phase  and  amplitude  modules  using  planar  interfaces  and  indicates  the  appar¬ 
ent  simplicity  of  the  method. 
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Integrated  Planar  Switch  Matrix 
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(ti>  Network  Schematic 


(b)  Variable  Power  Combiner  Configurations 


Figure  9.  Seven  Port  Amplitude  and  Phase  Network 
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Figure  10.  Beam  Port  Phase  and  Amplitude  Element 


Figure  U.  37  Element  Grouping  of  Phase  and  Amplitude  Modules 


3.  LENS  DESIGN 


Several  uniform  dielectric  lenses  were  considered  and  a  variant  of  the  lens  originally 
considered  by  Shinn  [1956]  and  more  recently  investigated  by  Rotman  [1983,19841  was 
selected.  This  type  of  lens  approximates  the  performance  of  an  ideal  thin  lens  which 
exhibits  the  property  that  the  coma  error  is  small  in  comparison  to  the  spherical  error 
[Cheston  and  Shinn,  1952).  The  lens  considered  by  Shinn  has  the  property  that  the  outer 
surface  is  spherical  with  the  center  coinciding  with  the  focal  point.  Referring  to  Figure  12, 
the  lens  is  designed  by  using  Snell’s  law  at  the  spherical  surface  (surface  2)  and  using  the 
equipath  condition  for  all  rays.  The  more  difficult  problem  is  to  determine  the  zone  step 
locations.  This  is  accomplished  by  the  simultaneous  solution  of  two  coupled  nonlinear 
equations. 

A  more  generalized  lens  is  realized  by  permitting  the  outer  surface  radius  to  assume  other 
values.  This  generalized  lens  provides  designs  which  are  quite  similar  to  more  classical 
lenses  such  as  the  meniscus  lens  and  the  plano-convex  lens.  The  generalized  lens  design  is 
an  exact  duplication  of  the  hyperbolic  lens  when  the  radius  approaches  infinity.  Figure  13 
illustrates  the  lens  contours  of  several  lenses  considered  in  this  study.  It  is  clear  that  the 
lens  approaches  a  thin  spherical  lens  as  the  number  of  zone  steps  increases.  However,  an 
increase  of  the  number  of  zone  steps  results  in  a  decrease  of  the  bandwidth  of  the  lens 
since  the  zoning  is  exact  only  at  a  single  frequency.  For  the  present  study,  seven  zones  are 
normally  used.  An  added  benefit  of  zoning  the  lens  is  a  decrease  in  the  lens  weight. 
Table  3  illustrates  the  weight  and  on  axis  thickness  of  a  34  inch  generalized  lens  for 
various  number  of  zones.  A  lens  with  seven  zones  is  approximately  a  third  of  the  unzoned 
weight  and  represents  a  significant  weight  advantage. 


TABLE  3.  Generalized  Lens  Properties  vs.  Number  of  Zones 


Number  of  Zones 

On  Axis  Thickness 

Weight 

1 

5.353  in 

93.34  lbs 

2 

4,906 

83.02 

3 

4.459 

70.94 

4 

4.013 

60.11 

a 

3.566 

r.n  &4 

a 

3.119 

42.25 

7 

2.673 

35.26 

8 

2.226 

29.59 

9 

1.779 

25.26 

10 

1.333 

22.28 

D  =  34.0  in.,  f/D=1.5,  Minimum  Thickness^ 
R0=0.85f,  Rexolite  dielectric 

=  0.375  in., 
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4.  PERFORMANCE  ANALYSIS 


Analysis  Method 

The  analysis  of  the  radiation  properties  of  a  large  MBA  involves:  the  prediction  of  a  feed 
pattern,  ray  tracing  from  each  feed  through  a  large  zoned  dielectric  lens  to  a  planar 
aperture,  and  computation  of  the  far  field  pattern  using  an  FFT  algorithm  Perhaps  the 
most  important  element  of  the  analysis  is  the  ray  tracing  analysis.  Referring  to  Figure  1 4, 
a  projected  aperture  plane  is  defined  with  a  square  cell  periodic  gTid.  The  ray  tracing 
solution  seeks  the  intercepts  on  both  the  inner  and  outer  lens  surfaces  such  that  Snell’s 
law  is  obeyed  with  the  ray  passing  through  the  selected  grid  coordinates  and  the  source 
point.  This  solution  involves  the  simultaneous  solution  of  two  coupled  nonlinear  equations. 
The  phase  of  the  ray  follows  directly  from  :  (1)  the  total  path  length  with  a  consideration 
of  the  dielectric  media  along  ray,  and  (2)  the  phase  of  the  feed  pattern.  The  amplitude  is 
evaluated  by  the  absolute  feed  pattern  magnitude  in  the  direction  of  the  ray  path  and  the 
proper  treatment  of  the  divergence  factors  [Lee,  et.  al.,  1982;  Deschamps,  1972]  at  each 
interface  which  includes  the  transformation  of  the  two  principal  radii  of  curvature  of  the 
wavefront  as  it  crosses  dielectric  boundaries.  A  scalar  average  of  the  Fresnel  transmission 
coefficients  for  TE  and  TM  polarizations  is  included  at  each  interface.  The  loss  tangent  of 
the  dielectric  is  included  to  properly  account  for  dissipation  for  rays  within  the  lens. 

Feed  Clusters 


The  use  of  a  cluster  of  feeds  with  appropriate  amplitude  and  phase  excitations  permits 
considerable  flexibility  by  allowing:  (1)  alternate  feed  loci,  12)  tailoring  the  secondary  pat¬ 
tern  characteristics  such  as  sidelobes,  and  13)  adaptive  nulling. 


Three  feed  loci  are  considered:  <1)  a  spherical  surface  having  a  radius  equal  to  the  lens 
focal  length,  (2)  a  compromise  focal  surface,  and  (3)  a  planar  surface.  A  goal  of  the  present 
study  is  to  utilize  planar  surfaces  where  possible  to  ease  interface  difficulties  associated 
with  feed  loci  having  a  finite  curvature.  Figure  16  illustrates  single  feed  patterns,  for  a 
planar  feed  locus  and  a  compromise  feed  locus  at  the  edge  of  earth  coverage  using  a 
32.6  inch  lens.  (The  feeds  are  pointed  toward  the  lens  apex.)  The  defocusing  of  the  planar 
locus  is  quite  evident  while  the  compromise  focal  surface  produces  well  focused  patterns. 
Tiie  gain  of  the  component  beam  on  a  planar  focal  surface  is  -1.1  dB  below  the  compro¬ 
mise  locus.  (An  additional  gain  reduction  of  0.62  dB  is  incurred  when  the  feeds  are 
directed  parallel  with  to  the  lens  axis.)  In  order  to  overcome  the  planar  surface  defocusing, 
the  performance  of  this  same  antenna  with  seven  element  hexagonal  clusters  is 
considered.  Figure  16a  illustrates  the  far  field  patterns  at  8.6°  scan  for  a  seven  element 
feed  with  a  uniform  amplitude  and  phase  illumination.  Figure  16b  illustrates  this  same 
case  except  that  the  feed  phases  have  been  adjusted  for  cophased  component  beams.  The 
compensation  provided  by  the  cluster  is  obvious, 


Figure  17  illustrates  the  far  field  patterns  of  this  cophased,  seven  element  cluster  with  the 
outer  six  elements  attenuated  relative  to  the  center  element  by  various  amounts.  Sidelobe 
levels  of  better  than  -30  dB  are  achieved  when  the  edge  elements  are  approximately 
-12  dB  relative  to  the  center  element. 
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Figure  16.  Lens 


<b)  Phase  Compensation 


Patterns  of  a  Uniform  Amplitude  Seven  Element  Cluster 


Feed  cluster  control  can  also  be  used  for  adaptive  nulling.  The  steady  state  weights  in  the 
presence  of  jammers  are  related  to  the  quiescent  (i.e.  without  jammers)  weights  via  the 
channel  covariance  matrix  (Fenn,  1982).  The  eigenvectors  (or  eigenbeams)  and  eigenvalues 
of  the  covariance  matrix  are  useful  in  interpreting  the  adaptive  algorithm.  Only  the  larger 
eigenvalues  will  produce  adapted  weights  which  vary  significantly  from  the  quiescent 
weights.  The  algorithm  removes  the  projections  of  the  interfering  source  eigenbeams  from 
the  quiescent  weight  vector,  which  results  in  a  null  in  the  antenna  angular  response  in  the 
direction  of  the  interference.  Alternately,  the  adaptive  process  may  be  viewed  as  forming 
beams  in  the  direction  of  the  jammers  and  subtracting  the  eigenbeams  from  the  quiescent 
beam  to  form  a  null  in  the  jammer  direction.  Figure  18  illustrates  this  process  for  the  a 
large  MBA  with  a  seven  element  hexagonal  cluster  using  a  34  inch  diameter  lens  at 
44.5  GHz.  A  user  is  assumed  at  the  origin  while  a  jammer  is  assumed  at  -0.2897*.  The 
quiescent  beam,  the  dominant  eigenbeam  (i.e.  the  largest  eigenvalue),  and  the  adapted 
beam  ate  shown.  It  is  clear  that  the  dominant  eigenbeam  has  a  maximum  in  the  direction 
of  the  jammer  and  that  the  adapted  beam  has  a  well  formed  null  in  the  direction  of  the 
jammer. 


5.  BFN  COMPARISON 

The  constrained  BFN  and  the  beam  port  array  are  compared  with  a  number  of 
considerations,  including;  (1)  weight  and  size,  (2)  gain  and  noise  temperature,  (3)  power 
requirements,  (4)  complexity  and  cost,  (5)  adaptive  nulling,  (6)  reliability,  and  (7) 
robustness. 

Weight  and  Size  i  Power  Requirements 

The  weight  of  a  satellite  based  antenna  is  one  of  its  most  critical  properties,  but  even  a 
lightweight  antenna  can  have  hidden  weight  penalties  associated  with  the  generation  of 
power  or  the  support  structure  of  a  large  antenna.  For  this  reason  it  is  important  to 
consider  these  parameters  together. 

Table  4  illustrates  the  estimated  weight,  power,  and  size  of  an  859  port  antenna  system. 
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Approach 

DP 

Weight 

Elect. 

Weight 

Total 

Power 

Added  Total 

System  Wt*  Weight 

Size 

Constrained 

119.56 

75.55 

29.34 

20.54 

215.65 

33  X  73 

Beam  Port 

76.36 

11.42 

219.30 

153.50 

241.28 

33  X  73 

1  Weight  (lbs),  Power  (Watts),  Size  (inches) 

1  Weight  attributed  to  prime  power  source  at  0.7  Ib/Watt 
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QUIESCENT  BEAM 


These  estimates  are  based  on  a  comprehensive  examination  of  these  antenna  systems 
utilizing  advanced  state-of-the-crt  techniques.  For  the  constrained  BFN  the  estimate  in¬ 
cludes:  waveguide;  horns,  polarizers,  transitions,  active  switches,  dummy  switches,  isola¬ 
tion  switches,  flanges,  a  32.6  inch  Rexolite  dielectric  lens,  a  seven  port  VPC  network  with 
330°  phase  shifters,  support  structure,  switch  and  VPC  drivers  (with  full  standby  redun¬ 
dancy),  central  electronics  (with  full  redundancy),  and  electronic  enclosure  weight  (includ¬ 
ing  EMP  shielding),  For  the  beam  port  array  the  estimate  includes:  horns,  polarizers, 
transitions,  phase  and  amplitude  modules,  a  32.6  inch  Rexolite  dielectric  lens,  a  511 
element  focused  phased  array  feed,  support  structure,  phase  and  amplitude  module 
electronic  drivers  (with  full  redundancy),  phased  array  drivers  (with  full  redundancy),  cen¬ 
tral  electronics  (with  full  redundancy),  and  electronic  enclosure  weight  (including  EMP 
shielding).  The  total  hardware  weight  is  195  lbs  for  the  constrained  BFN  and  only  88  lbs 
for  the  beam  port  array.  Hence,  there  initially  appears  to  be  a  significant  weight  savings 
for  the  beam  port  concept.  However,  the  power  requirements  of  the  beam  port  array  are 
considerable  in  comparison  to  the  constrained  BFN,  primarily  due  to  the  power  associated 
with  the  phased  array  feed.  This  results  in  an  increase  in  the  estimated  satellite  weight  by 
approximately  0.7  lb  per  Watt  of  the  required  power  (Cummings,  1987],  When  this  total 
system  weight  is  considered,  the  total  hardware  weight  is  216  lbs  for  the  constrained  BFN 
and  241  lbs  for  the  beam  port  array  -  a  different  conclusion  than  originally  anticipated 
This  conclusion  may  vary  aramatieally  if  the  possibility  of  solid  state  components  are 
considered,  however,  the  impact  on  the  power  supplies  will  still  be  significant. 

Gain  I  Noise  T emperature 

The  two  antenna  systems  are  identical  in  the  selection  of  the  dielectric  lens,  feed  horns, 
and  feed  locus.  Hence,  the  antenna  gain  relative  to  this  radiating  portion  of  the  antenna 
does  not  distinguish  the  antennas.  However,  the  insertion  losses  of  the  remaining  beat., 
forming  network  provide  one  basis  for  a  quantitative  tradeoff.  The  estimated  insertion  loss 
of  the  constrained  BFN  is  5.38  dB.  For  the  beam  port  array,  the  estimated  insertion  loss 
is  9.4  dB  near  broadside  and  10.6  dB  at  the  edge  of  scan.  (A  large  portion  of  the  loss  is 
associated  with  the  coupling  of  the  beam  port  array  to  the  focused  phased  array  —  which 
is  aggravated  by  the  size  of  the  beam  port  horn  elements.)  An  alternative  beam  port  array 
with  a  single  horn  feed  was  also  examined  —  but  exhibits  an  insertion  loss  of  22.1  dB  at 
broadside  and  25.5  dB  at  the  edge  of  scan. 

The  use  of  low  noise  amplifiers  (LNAs)  can  also  be  considered.  The  natural  location  for  the 
LNAs  for  the  constrained  BFN  is  at  the  input  to  the  seven  channel  VPC  network  and 
results  in  a  loss  of  1.93  dB  ahead  of  the  amplification.  For  either  of  the  beam  port 
concepts,  the  LNAs  are  considered  ahead  of  the  phase  and  amplitude  control  modules  with 
only  a  loss  of  0.2  dB  ahead  of  the  amplification,  The  properties  of  the  LNA  are  assumed 
to  be:  10  dB  gain  per  stage  and  a  4  dB  noise  figure.  Table  5  illustrates  the  noise  tem¬ 
perature  of  the  antenna  systems  based  on  the  beam  forming  network  loss  and  neglecting 
the  antenna  temperature.  The  system  noise  temperature  is  shown  as  a  function  of  the 
number  of  LNA  gain  stages.  Clearly,  the  beam  port  concept  relies  on  the  use  of  LNAs  to 
achieve  an  acceptable  result  and  is  contrasted  with  the  constrained  BFN  which  can  be 
used  without  LNAs.  Note  that  single  stage  amplifiers  will  suffice  for  either  the  constrained 
BFN  or  the  beam  port  array  with  the  phased  array  feed.  However,  the  beam  port  array 
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with  a  horn  teed  requires  a  two  or  three  stage  amplifier  for  satisfactory  operation.  (These 
choices  are  indicated  in  Table  6.)  The  prospects  for  the  development  of  a  multistage  LNA 
at  EHF'  with  precise  amplitude  and  phase  tracking  over  the  required  temperature  and 
frequency  band  are  remote  in  the  near  future. 


TABLE  5.  Noise  Temperature  Tradeoff 


Antenna  System 

No.  Gain 
Stages 

Temp  (“K)‘ 

Constrained  BFN 

0 

710.9 

1 

900.0  <- 

2 

851.5 

3 

846.6 

BPA  (Phased  Array 

0 

3039.6 

Feed) 

1 

776.4  <- 

2 

503.0 

3 

476.8 

BPA  (Horn  Feed) 

0 

102606.9 

1 

10732.0 

2 

1498.7 

3 

575.4  <- 

4 

483.0 

1  <-  indicates  selected  number  of  gain  stages 


A  measure  frequently  used  for  defining  the  receiving  properties  of  a  system  is  the  gain  to 
temperature  ratio,  G/T,  where  T  is  the  effective  system  temperature  at  some  reference 
point  and  O  is  the  gain  prior  to  this  reference.  A  nominal  gain  when  a  cluster  is  excited  is 
44  dB  for  the  863  port  case  when  the  outer  ring  is  -10  dB  relative  to  the  center  and  a 
planar  surface  is  used.  Thus  typical  values  of  O/T  are.  14.5  dB  for  the  constrained  BFN, 
15.1  dB  for  the  beam  port  array  with  a  phased  array  feed.  Without  the  use  of  LNAs,  the 
resulting  O/T  ratios  are  16.6  dB  and  9.2  dB  for  the  constrained  BFN  and  beam  port 
array,  respectively. 

Complexity  I  Cost 

While  the  cost  of  any  future  system  is  subjective  at  best,  the  complexity  is  quite  objective. 

The  859  port  constrained  BFN  has  1986  switches  (1711  active  switches  including  isolation 
switches),  19  phase  shifters  associated  with  the  VPC  amplitude  and  phasing  network,  and 
7  LNAs  at  the  input  to  the  7  port  VPC  network.  This  results  in  a  total  of  1730  control 
functions. 
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The  859  port  beam  port  array  has  1718  high  precision  phase  shifters  associated  with  the 
phase  and  amplitude  control  modules,  611  modest  precision  phase  shifters  associated  with 
the  focused  phased  array,  and  859  bias  switches  for  on/off  control  of  the  LNAs  for  added 
port  isolation.  This  results  in  a  total  of  3088  control  functions. 

In  addition,  to  the  increased  number  of  controls  for  the  beam  port  array  concept,  the  type 
of  controls  for  this  approach  are  also  considerably  more  complex  (i.e.  an  8  bit  phase  shifter 
driver  in  comparison  to  a  switch  driver). 

Assuming  the  hardware  cost  is  approximately  proportional  to  the  number  of  controls,  the 
cost  of  the  beam  port  array  is  at  least  twice  that  of  the  constrained  BFN.  In  a  practical 
consideration,  the  cost  differential  could  be  even  greater  because  of  the  added  complexity 
of  the  phrse  shifters  compared  to  the  switches. 

Adaptive  Nulling 

Adaptive  nulling  properties  of  the  constrained  BFN  versus  the  beam  port  antenna  are 
identical  when  the  losses  of  the  networks  are  not  considered.  However,  there  can  be  a 
significant  difference  in  the  implementation  of  the  adaptive  algorithm  as  well  as 
differences  due  to  the  configuration  of  and  the  practical  errors  of  the  components. 

F'or  a  seven  port  nulling  algorithm,  the  constrained  BFN  requires  the  control  of  seven 
360°  phase  shifters  and  six  VPCs  (another  12  phase  shifters)  -  a  total  of  19  controls.  Tho 
beam  port  array  requires  the  control  of  14  phase  shifters  and  hence  does  offer  some 
advantage. 


The  means  uf  extracting  an  error  signal  are  significantly  different.  Typically,  directional 
couplers  are  used  at  the  input  of  the  VPC  nulling  network  to  extract  the  channel  samples 
for  a  constrained  BFN.  The  beam  port  array  does  not  permit  this  method  of  obtaining  a 
channel  sample.  Instead  a  dither  of  the  desired  channel  can  be  used  to  modulate  specific 
ports  in  amplitude  and  phase.  The  inherent  simplicity  of  the  dither  method  appears  advan¬ 
tageous  compared  with  the  more  complex  structure  of  the  VPC  network  used  with  a 
constrained  BFN. 


Reliability 

The  complexity  of  the  large  MBA  is  considerable  in  comparison  to  even  the  most  complex 
MBA  by  today’s  standards.  Hence,  the  issue  of  reliability  is  an  important  concern  for 
either  approach.  Total  standby  redundancy  has  been  assumed  for  ail  electronics  for  either 
approach  —  however,  triple  redundancy  may  be  required  for  key  components. 

An  interesting  aspect  associated  with  the  reliability  of  a  beam  for  a  very  high  gain  MBA  is 
that  finite  beam  failures  can  be  tolerable.  Since  the  antenna  has  a  large  number  of  beams 
with  very  close  spaemgs,  it  mey  be  possible  to  repoint  the  satellite  only  slightly  and  utilise 
other  operational  beams  with  only  a  slight  degradation  of  antenna  performance. 
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Robustness 


In  the  present  application,  robustness  refers  to  the  ability  of  the  antenna  system  to  be 
designed  for  a  wide  variety  of  applications.  In  this  sense,  both  antenna  systems  under 
consideration  are  quite  robust.  However,  the  beam  port  antenna  allows  additional  robust¬ 
ness  with  regard  to  reconfiguration  of  the  beam  patterns  since  the  central  electronics  can 
be  sufficiently  general  as  to  allow  increasing  the  cluster  size  and  shape  beyond  c  seven 
port  hexagonal  layout.  This  is  not  possible  with  the  constrained  BFN. 

8.  CONCLUSIONS 

This  paper  has  examined  the  general  goal  of  achieving  practical  designs  for  large  multiple 
beam  antennas  for  EHF  MILSATCOM.  The  two  basic  antennas  under  consideration  are:  a 
constrained  BFN  and  a  space  feed  antenna.  The  constrained  BFN  antenna  represents  a 
more  traditional  antenna  concept  using  mature  ferrite  waveguide  switches  and  waveguide 
feed  networks.  The  beam  port  antenna  represents  an  approach  which  potentially  can  offer 
weight  and  layout  advantages.  However,  recent  technology  has  made  the  constrained  BFN 
attractive,  even  for  a  large  MBA.  In  fact,  a  constrained  BFN  with  several  hundred  feeds 
and  seven  port  feed  clusters  is  attainable  with  today’s  technology;  and  in  particular,  in 
terms  of  complexity,  weight,  and  power  consumption,  a  large  constrained  BFN  is  competi¬ 
tive  with  less  capable  concepts  presently  under  consideration  for  EHF  MILSATCOM. 

The  beam  port  concept  with  a  focused  phased  array  feed  can  offer  equally  acceptable 
microwave  performance  with  less  microwave  hardware.  However,  with  presently  available 
electronic  designs,  the  power  consumption  of  the  ferrite  scanned  phased  array  tends  to 
make  this  approach  less  attractive  than  the  constrained  feed.  With  substantial  advances  in 
EHF  L,NA  technology,  the  phased  array  could  be  replaced  with  a  single  horn,  or  the  ferrite 
array  replaced  with  a  solid  state  MMIC  array  design  that  consumes  less  power.  Given 
these  advances,  which  may  result  from  the  current  industry  wide  emphasis  on  MMIC,  the 
space  feed  could  offer  significant  advantages  in  terms  of  reduced  weight  and  increased 
beam  forming  flexibility. 
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MEASURED  ANTENNA  PERFORMANCE  OF 


AN  EHF  HYBRID  SCAN  ARRAY  ANTENNA  FOR 
MILSATCOM  APPLICATIONS 


J.  F.  Pedersen 

Hazeltine  Corporation 
Wheeler  Laboratory 
Greenlawn,  NY  11740 

ABSTRACT 

Hazeltine  has  completed  a  contract  to  design  an  Advanced 
Development  Model  of  an  EHF  Hybrid-Scan  Array  for  potent ial 
application  to  airborne  MILSATCOM  terminals.  This  paper  will 
present  the  design  and  electrical  tests  on  the  antenna  assembly. 

The  antenna  assembly  .includes:  a  dual-band  aperture  consist¬ 
ing  of  waveguide  radiators  with  isolating  filters  to  prevent  coup¬ 
ling  between  receive-band  and  transmi t-band  radiators;  a  meander¬ 
line  circular-wave  polarizer;  a  conformal  dual-band  radome  capable 
of  passing  circular  pomcized  waves  over  a  70°  range  of  incidence 
angles;  a  dual-band  rotary  joint  and  the  associated  electronics 
and  phase  shifters  to  stabilize  and  steer  the  antenna  beam  over 


nearly  a  hemisphere. 


SECTION  1.0 


INTRODUCTION 

The  task  of  providing  high  performance  airborne  platforms  with 
SATCOM  capability  has  involved  one  major  conflict.  Specifically! 
the  user  would  like  to  have  this  communication  capability  with  no 
antenna  protrusion,  for  reasons  ’.’elating  to  aircraft  signature, 
drag,  and  blockage  of  other  antennas.  However,  the  need  for  near- 
hemispheric  high  gain  coverage  has  required  significant  antenna 
protrusion  and  large  radomes.  Although  conformal  fully  phased 
arrays  are  bcrng  developed,  they  are  not  presently  viable  candi¬ 
dates  for  a  practical  terminal. 

Iiazcltine  has  pursued  the  technology  of  hybrid  scanning 
(electronic  scanning  in  elevation  and  mechanical  scanning  in 
azimuth)  as  a  cost-effective  solution  to  the  protrusion/gain 
trade-off^- . 

hi  SHF,  an  exploratory  model  of  a  receive-only  antenna2  and  an 
advanced  development  mode,  of  a  fly-able,  full  duplex  antenna-®  have 
employed  a  horizontal  "turntable"  aperture  with  electronic 
scanning  from  broadside  to  end-fire. 

At.  EKF,  however,  end-fire  scanning  is  not  practical  for  the 
gain  levels  required.  The  EHF  design  concept  uses  a  tilted  array, 
which  is  located  below  the  level  of  the  fuselage  (figure  1-1). 

The  antenna  array  elec.troni.cal!  y  scans  the  beam  ±33°  in  elevation, 
while  it  rotates  mechanically  360°  in  azimuth.  It  is  protected  by 
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Figure  1-1.  EHF  Hybrid  Scan  Antenna  Configuration  With 

Flush  Radome 

a  radome  that  is  contoured  to  the  aircraft  surface.  The  installa¬ 
tion  is  flush,  meeting  performance  criteria  from  zenith  to  20° 
above  the  horizon4. 

In  addition  to  addressing  the  major  antenna  design  issues, 
this  program  has  also  considered  installation  of  the  antenna  on  an 
aircraft  such  as  the  Bl-B  including  preliminary  structural  analy¬ 
sis.  This  latter  effort  has  been  conducted  by  Electrospace 
Systems,  Inc.  (ESI)  as  a  team  member. 

The  realisation  of  the  design  involves  the  development  of 
several  key  components  Including: 

•  A  dual  frequency  band  radome,  which  must  operate  at  20  GHz 
and  44  GHz  with  circularly-polarized  waves  and  with 
incidence  angles  ranging  from  0°  to  70° 
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•  A  dual  frequency  band  aperture  of  interleaved  20  GHz  and 
44  GHz  slot  elements 

•  A  dual  frequency  band  polarizer,  which  must  operate  over  a 
range  of  incidence  angles  from  0°  to  35° 

SECTION  2.0 

DESIGN  OF  THE  KEY  COMPONENTS 

2 . 1  RADOME 

If  the  phased  array  were  horizontal  and  flush  with  the  air¬ 
craft  surface,  it  could  be  protected  from  the  environment  by 
dielectric  windows  in  the  radiating  elements  or  with  a  radome 
attached  to  the  rotating  radiating  surface.  However,  with  the 
internal  location  of  the  antenna,  an  independent  radome  is 
required  that  can  be  made  flush  with  the  aircraft  surface. 

The  flush  radome  must  operate  at  a  small  band  of  frequencies, 
about  20  GHz  and  44  GHz,  with  circularly  polarized  waves  and  with 
incidence  angles  ranging  front  0°  to  70°.  The  radome  must  provide 
high  transmission  efficiency  and  minimum  degradation  of  the  axial 
ratio  of  circular  polarization.  Standard  radome  designs  do  not 
appear  to  be  satisfactory.  Electrically  acceptable  designs  have 
poor  mechanical  strength,  while  radomt  that  are  adequate 
mechanically  have  degraded  electrical  performance. 

Acceptable  electrical  performance  and  adequate  mechanical 
strength  has  been  achieved  by  using  a  double-quarter  wavelength 


(at  20  GHz)  dielectric  sandwich  with  a  shunt  susceptance  between 
the  walls.  The  susceptance  produces  an  effect  that  is  large  at  44 
GHz  and  small  at  20  GHz,  giving  independent  control  of  the  radome 
match  angles.  The  computed  performance  of  the  radome  is  shown  in 
figures  2-la,  2-lb,  2-lc,  and  2-ld.  Figure  2-la  shows  the  reflec¬ 
tion  coefficient  for  H-plane  incidence  at  20  GHz  and  44  GHz.  Note 
that  the  match  angles  are  nearly  equal,  providing  low  reflection 
at  20  GHz  and  44  GHz  over  a  wide  range  of  incidence  angles. 

Figure  2-lb  shows  the  transmission  phase  difference  between 
orthogonally  polarized  waves  at  20  GHz  and  44  GHz.  Figure  2-lc 
shows  a  plot  of  axial  ratio  as  a  function  of  incidence  angle.  At 
70°  incidence  and  44  GHz,  the  axial  ratio  is  3.5  dB;  this  is  a 
result  of  an  increase  in  the  transmission  phase  difference  at  70° 
shown  in  figure  2-lb.  At  all  angles  less  than  70°,  the  axial 
ratio  is  significantly  lower. 

As  shown  in  figure  2-2,  the  measured  insertion  loss  of  the 
2-foot  by  2-foot  radome  panel  is  less  than  0.5  dB  at  20  GHz  and 
less  than  0.75  dB  at  44  GHz  over  the  full  range  of  incidence 
angles.  The  measured  performance  should  be  compared  to  the 
computed  loss  shown  in  figure  2-ld. 

2.2  INTERLEAVED  20  GHZ,  44  GHZ  APERTURE 

Linearly  polarized  elements  and  a  wave  polarizer  are  used  to 
achieve  circular  polarization.  Separate  elements  are  used  for  the 
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transmission 


Ficure  2-1.  Computed  Radome  Performance 


ANGLE  OF  INCIDENCE  (DEGREES) 

20  40  60  80 

I - 1 - 1 - 1 - 1 


of  Incidence 

receive  band  and  the  transmit  band.  The  elements  are  interleaved, 
forming  a  single  aperture  assembly. 

The  design  of  a  dual-band  array  requires  a  technique  for  pre¬ 
venting  the  high  band  power  from  coupling  into  the  low  band  radia¬ 
tor  This  can  be  accomplished  in  two  ways. 

One  method  uses  cross-polarized  elements.  Two  possible 
arrangements  of  cross-polarized  slots  are  shown  in  figure  2-3.  In 
either  arrangement,  the  net  component  of  high  band  current  flews 
along  the  length  of  the  low  band  element.  If  the  low  band  element 
is  thin  and  below  cut-off,  it  is  not  excited  at  the  nigh  band 
frequency. 

To  scan  the  array,  a  phase  shift  is  introduced  to  each  row  of 
slots.  A  component  of  current  can  flow  between  rows  of  slots 
caused  by  this  phase  difference.  The  arrangement  of  slots  shown 
in  figure  2-3a  does  not  couple  to  this  current,  and,  therefore,  it 
is  the  preferred  arrangement  of  cross-polarized  elements. 
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b. 

Figure  2-3.  Two  Arrangements  of  Cross-Polarized  Slots  Showing 
How  High  Band  Currents  Cancel  Transverse  to  the  Long  Dimension 

of  Low  Band  Slot 

The  second  method  ot  decoupling  the  low  band  radiator  from  the 
high  band  power  is  to  use  a  filter  that  is  internal  to  the  low 
band  element.  Either  cross-polarized  or  co-polarized  elements  may 
be  used  with  this  approach.  The  filter  would  be  designed  to  put 
an  effective  short-circuit  plane  at  the  face  of  the  low  band 
element  in  both  the  fundamental  TE^g  and  TF^g  modes  at  the  high 


bond  fiequency.  The.  filter  would  also  pass  the  lew  hand  frequency 
with  small  reflection  loss. 

The  element  spacing  between  rows  of  low  hand  elements  is 
determined  by  the  scan  requirements  and  the  tiiyh  band  frequency. 
The  preferred  arrangement  of  cross-polar ized  elements  has  the  long 
dimension  of  the  low  band  element  in  the  plane  that  has  the  small¬ 
est  interelen.ont  spacing.  Consequently,  the  low  band  element  will 
be  electrically  small'’.  The  low  band  element  is  also  scanned  in 
the  H-plane,  which  reduces  its  radiation  power  factor  further*’. 

To  increase  the  -adiation  power  factor  and  maintain  a  low 
reflection  coetf icient ,  the  low  band  element  area  must  be 
increased.  This  results  in  a  nearly  square  element  rather  than  a 
thin  element,  which  is  desirable  for  minima r  coupling.  For  these 
reasons,  c  co-polarized  arrangement  of  slots  wirh  decoupling  fil¬ 
ters  in  the  low  band  elements  was  selected  for  the  antenna 
ape* 1 1 u  r c . 

The  computed  filte,.  perform  nee  is  shown  in  figure  2-4.  Note 
the  pass  band  and  two  reject  bands  in  figure  2-4;  one  reject  band 
cor  responds  to  the  fundamental  high  band  waveguide  mode  while  the 
second  reject  band  represents  the  TE2q  mode  response. 

Figure  2-4  shows  the  phase  slope  of  the  filter  across  its  pass 
bard.  Since  the  filter's  phase  slope  is  not  flat  across  the  pass 
band,  a  shift  in  the  center  frequency  of  the  filter  caused  by 


tolerances  can  introduce  a  phase  error  across  the  aperture.  The 


Figure  2-4.  Reflection  Loss  vs  Frequency  of  Isolating  Filter 
amount  of  phase  error  was  evaluated  as  a  random  function  of 
expected  machining  tolerances.  The  result  is  shown  in  figure  2-5 
as  a  histogram.  Each  bin  is  approximately  1.5°  wide;  a  total  of 
2000  trials  were  made  with  the  number  of  trials  indicated  within 
each  bin.  The  phase  error  introduced  by  the  filter  can  produce 
error  sidelobes  at  a  -27  dB  Level.  This  error  is  acceptable  since 
it  is  5  dB  below  the  design  sidelobe  level. 

The  computed  filter  performance  was  verified  by  fabricating  a 
small  section  of  the  array.  Coupling  measurements  and  radiation 
pattern  measurements  were  recorded  with  and  without  the  filter  in 
the  low  band  element.  Samples  of  these  results  are  shown  in  the 
next  four  figures. 

Figure  2-6a  shows  the  coupled  high  band  signal  into  the  low 
band  waveguide.  It  is  at  a  maximum  level  of  -5  dB  in  the  TE1Q 
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Figure  2-5.  Histogram  of  Filter  Phase 
mode.  Figure  2-6b  shows  the  coupled  signal  with  the  filter.  The 
coupled  signal  level  has  been  reduced  to  -27  dB. 

Radiation  pattern  measurements  were  also  recorded.  Sample 
results  are  shown  in  figures  2--7  and  2-8.  Note  the  significant 
pattern  degradation  in  figure  2-7  without  the  filter  and  the 
improvement  shown  in  figure  2-8  with  the  filter.  A  control  pat¬ 
tern  is  shown  on  each  figure;  it  corresponds  to  a  pattern  of  the 
high  band  array  without  the  low  band  elements  present. 

2.3  POLARIZER 


The  circular  polarizer  converts  the  polarization  of  an  inci¬ 
dent  wave  for  linear  to  circular  at  both  20  GHz  and  44  GHz,  over  a 
range  of  incidence  angles  from  0°  to  ±35°. 

An  approach  to  designing  a  polarizer7  uses  printed-circuit 
elements  that  introduce  inductive  and  capacitive  suseeptance  to  an 
incident  wave.  A  variation  of  this  design^  uses  inductive 
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MEASURED  COUPLING  INTO  LOW  BAND  ELEMENT-NO  FILTER 


measured  coupling  to  low  band  element  with  filter 


Figure  2-6.  Coupled  Signal 

elements,  which  are  serpentine  rather  than  straight.  This  latter 
approach  (a  meander-line  polarizer)  was  selected  for  the  £HF 
polarizer. 

The  meander-line  polarizer  produces  good  circular  polarization 
over  a  wide  frequency  range.  This  is  accomplished  by  using 
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polarization-selective  elements  that  have  opposite  frequency  char¬ 
acteristics,  thereby  maintaining  a  nearly  constant  90°  phase  shift 
over  the  frequency  band. 

There  are  several  problems  involved  with  the  design  of  a  cir¬ 
cular  polarizer  for  operation  over  a  range  of  incidence  angles. 
These  include:  the  change  of  effective  wavelength  with  incidence 
angle;  the  angle-sensitive  behavior  of  the  polarizer  elements;  and 
the  geometric  situation  that  precludes  maintaining  the  desired  45° 
orientation  of  the  polarizer  elements  relative  to  the  polarization 
of  the  incident  wave  over  the  range  of  incidence  angles. 

The  last  problem  can  be  partially  overcome  by  orienting  the 
polarizer  elements  at  an  angle  that  deviates  from  45°  and  that 
represents  the  best  compromise  over  the  range  of  incidence  angle. 

The  requirement  for  a  combination  of  both  wide-band  and  wide- 
angle  performance  can  be  approached  by  designing  the  polarizer  for 
operation  at  only  the  two  5-percent  frequency  bands,  about  20  GHz 
and  44  GHz,  rather  than  over  the  complete  band  between  these* 
frequencies . 

Figures  2-9a  through  2-9c  show  the  computed  polarizer  perform¬ 
ance.  Three  curves  are  shown  on  each  graph  for  three  angles  of 
incidence.  The  frequency  bands  of  interest  are  noted  on  each 
figure.  Figures  2-9a  and  2-9b  show  the  reflection  coefficient  of 
the  principal  polarizations.  Figure  2-9c  shows  the  resulting 
axial  ratio  of  the  polarizer.  The  measured  axial  ratio  of  the 
polarizer  vs  angle  of  incidence  is  shown  in  figure  2-10. 
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Figure  2-9.  Computed  Polarizer  Performance 
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Fimjie  2-10.  Measured  Axial  Ratio  ot  Polarizer  vs  Angle 
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2 . 4  POWER  COMBINERS 

For  an  array  that  scans  in  only  one  piano,  the  elements  may  be 
arranged  in  a  rectangular  grid.  'i’r.e  clement-  spacing  in  v he  plane 
of  scan  is  smaller  to  prevent  a  radiating  giatin]  lobe.  Even  for 
a  moderate  scanning  range  (SI0),  the  scar,  plant;  e ,  orient  spacing 
will  preclude  the  use  oi  standard  vj.j--r.gu d »  as  tat  feed.  St  rip¬ 
line,  suspeneer.  strip!  i»»e,  riego  va-eg aid  ?,  and  t  educed  height 
waveguide  are  alternatives. 

The  antenna  consists  of  many  coirpc-'eM  3  connected  in  series. 
The  dissipative  loss  is  additive  and  it  can  be  apprt.-cic.bj.e  in 

t  ^  A-  _  1  - -  :  .1-..  -  4-  TPT  T  IT*  P  4  /-t^  L’ t  $  -  4  «r/^  m  r-  -  1  ^4-'  <r 
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n  K.  u.*U  *-  4.  V. 


total  antenna  loss  is  required  to  achieve  prac.-ical  ERP  and  C/T 
levels.  Since  reduced  height  waveguide  yields  the  lowest  dissipa¬ 
tive  loss,  all  power  combiners  use  this  medium. 

Sidelobe  control  was  required  for  each  principal  plane  of  the 
array.  In  azimuth,  the  amplitude  taper  ;s  achieved  by  using  addi  - 
tional  Icuels  of  equal  (3  dll)  splitters.  The  extra  levers  quan¬ 
tize  the  illumination  in  coarse  3  dll  steps;  however,  the  modest 
sidelobe  requirement  of  -16  dB  maximum  in  1 ransmit  and  -22  dB  in 
receive  can  be  achieved  bv  using  this  technique. 

In  elevation,  a  coarse  approximation  to  the  desired  illumina¬ 
tion  will  not  yield  these  same  sidelobe  levels-  this  is  because 
the  array  contains  fewer  rows  of  elements  than  elements  per  row. 
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As  a  result,  nonequal  E-plane  waveguide  dividers  are  used  to 
derive  the  amplitude  illumination  in  elevation. 

Figure  2-11  shows  measured  azimuth  patterns  (nonscanning 
plane)  of  the  array  including  radome  at  35°  elevation.  Sidelobe 
levels  as  a  function  of  elevation  are  shown  tabulated  in  figure 
2-11.  Figure  2-12  shows  typical  measured  patterns  in  the  scanning 
of  the  array. 


Figure  2-11.  Typical  Measured  Azimuth  Plane  Radiation  Patterns 
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SECTION  3.0 


ARRAY  CONSTRUCTION 

Figure  3-1  is  a  photograph  of  the  antenna  pedestal.  Tire 
p  .'Stal  consists  of  two  housings  -  an  inner  housing  that  rotates 
and  a  stationary  outer  housing.  The  millimeter  wave  components, 
the  electronics  to  control  and  stabilise  the  beam,  and  the  phase 
shifter  driver  cards  are  assembled  in  the  inner  housing. 

The  stationary  housing  contains  the  motor  stator  and  all.  power 
supplies  needed  to  drive  the  electronics  and  phase  shifters. 

A  blower,  attached  to  the  top  of  the  inner  housing,  draws  air 
up  through  the  pedestal  and  exhausts  the  hot  air  through  the  array 
of  holes  at  the  periphery  of  the  pedestal  as  seen  in  figure  3-1. 

A  slip  ring  routes  control  signals  and  power  supply  voltages 
between  housings.  A  dual-band  20  GHz/44  GHz  waveguide  rotary 
joint  (manufactured  by  Kevlin,  Inc.)  mounts  inside  the  slip  ring. 

The  beam  is  scanned  electronically  by  using  ierrite  phase 
shifters,  supplied  by  Electromagnetic  Sciences,  Inc. 

A  phase  shifter  assembly  contains  16  phase  shifters.  Three 
assemblies  (two  44-GHz  and  one  20-GHz)  aie  arranged  across  the 
azimuth  power  divider,  feeding  the  array  in  elevation.  The  driver 
cards  are  located  in  two  bays  and  connect  to  the  KF  assembly  bv  a 
flexible  balanced  strip  transmission  line. 
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Finn  re  3-1.  EHF  Hybrid  Scan  Array 
Reduced  height  waveguide  power  dividers  are  used  Cor  power 
distribution  to  the  phase  shifters  and  to  the  elements.  The 
polarizer  is  mounted  in  a  frame  0.25  inch  in  front  of  the  array. 

The  cadome  for  this  program  is  a  flat  2  foot  by  2  foot  panel 
mounted  to  a  B-52  fuselage  via  an  adaptor  plate  that  conforms  to 
the  fuselage's  cylindrical  shape. 


SECTION  4.0 
CONCLUSIONS 


Under  this  program,  all  the  technologies  required  to  provide  a 
flush  antenna  system  for  an  airborne  MILSATCOM  terminal  have  been 
reduced  to  practice.  The  array  is  a  cost-effective  solution  for 
aircraft  where  a  flush  conformal  array  is  needed. 
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ABSTRACT: 

The  measurement  of  one  dimensional  planar  arrays  that  are  electrically  long 
is  currently  a  difficult  task.  The  accurate  measurement  of  the  antenna  patterns 
is  limited  by  large  antenna  range  requirements  for  far  field  ranges,  by  the 
number  of  near  field  points  that  must  be  measured  for  near  field  probes,  and  by 
the  great  size  and  quality  of  the  reflector/feed  for  compact  ranges.  This  paper 
describes  the  use  of  quadratic  and  higher  phase  correcticn  applied  to  the  phased 
array. 

The  errors  that  result  from  this  approach  are  analyzed  as  a  function  of 
range,  antenna  electrical  length,  sidelobe  level  and  degree  of  correction. 

Both  computer  simulations  and  analysis  are  used  to  generate  the  limitations 
of  this  approach.  Far  field  measurements  made  in  the  near  field  are  currently 
planned . 

INTRODUCTION: 

The  Air  Force  is  currently  investigating  methods  of  mounting  large  lightweight 
phased  arrays  on  the  fuselages  of  long  aircraft  for  long  range  airborne  surveillance. 
Such  arrays  tend  to  be  very  long  in  one  dimension,  but  moderate  in  size  in  the 
other.  Because  of  the  size  of  these  antennas  the  measurement  of  their  basic 
patterns  is  very  difficult.  For  example  the  far  field  distance  for  a  10'  x  100' 
aperture  at  C-band  would  be  about  23  miles  away,  by  the  2D/A  criteria  for 
a  very  low  sidelobe  array  the  distance  would  be  much  greater  if  very  little 


pattern  distortion  is  desired.  Compact  ranges  for  this  size  antenna  require  a 
massive  high  tolerance  reflector  and  near  field  ranges  would  be  extremely  expensive. 
One  technique  to  measure  the  far  field  patterns  while  using  a  conventional 
antenna  range  is  based  on  the  fact  that  in  one  dimension,  the  smaller  one,  the 
far  field  range  is  at  a  reasonable  distance.  For  the  other  dimension  the 
antenna  is  well  within  the  Fresnel  field,  so  an  accurate  pattern  can  only  be 
found  by  putting  in  a  phase  correction  over  the  aperture-  Since  this  correction 
is  only  precisely  correct  at  one  angle,  repeated  corrections  are  therefore 
required  to  get  a  good  approximation  to  the  far  field  at  infinity  pattern.  In 
this  report  the  first  two  higher  order  correction  factors  are  examined  quadratic 
alone  and  quadratic  and  cubic.  An  analysis  of  the  residual  error  is  given  and 
from  that  analysis,  an  estimate  of  the  range  of  angles  is  over  which  a  particular 
correction  gives  small  addd  error  is  found  That  range  of  angles  over  which  the 
pattern  of  infinity  determines  the  number  of  resets  of  the  phase  shifters. 


SIMULATION: 


For  the  analysis  a  computer  simulation  was  developed.  This  involved  writing 
three  individual  computer  programs.  The  first  program  generates  the  actual  far 
field  pattern,  the  result  is  shown  in  fig  (1).  Using  the  pattern  as  a  reference 
a  second  program  was  written,  using  a  modified  expression  which  Includes  the 
correction  factors.  First  the  program  was  run  including  only  the  quadratic 
correction.  This  was  performed  for  several  different  cases.  First  the  distance 
R0  are  repeated.  The  program  was  then  run  including  the  third  order  correction 
and  as  with  the  quadratic  correction  the  parameters  R0  end  Gfe  are  varied.  To 
determine  the  degree  of  difference  between  the  actual  far  field  pattern  and  the 
pattern  which  include  the  correction  terms  another  simulation  was  run.  This 
program  subtracts  the  actual  far  field  values.  This  gives  a  measure  of  the  error 
in  adding  the  correction  terms. 

The  results  of  the  simulation  for  the  quadratic  term  correction  alone  are 
presented  in  fig  (2-5).  These  are  for  a  96>antenna.  For  the  addition  of  the 
third  order  correction  the  results  are  shown  in  fig  (6-9). 

With  the  results  of  the  computer  runs  a  3et  of  graphs  have  been  developed. 
These  graphs  are  a  measure  of  how  much  error  one  can  expect  In  applying  this 
technique  to  a  design. 

Fig  (10)  shows  how  the  error  will  vary  with  distance  in  applying  the 
quadratic  correction.  ( FI 1 )  shows  the  effect  of  adding  the  third  order  correction. 
The  graph  represents  how  many  times  the  phase  shifter  will  have  to  be  set  in  order 
to  achieve  a  certain  error  margin. 

ANALYSIS: 

In  order  to  estimate  the  residual  error  in  using  phase  correction 
in  the  Fresnel  field,  several  assumptions  are  made  to  derive  the  dominant  error 
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forms.  These  assumptions  include: 

1)  The  1/R  amplitude  term  can  be  approximated  by  1/RO,  that  the  Tresnel 
distance  chosen  is  large  enough  so  that  only  the  phase  term  contributes 
significantly  to  the  residual  error. 

2)  The  type  of  patterns  that  concern  us  are  low  sidelobe  patterns  with 
Taylor  illiminations  or  similar  and  beyond  a  small  angular  range  these 
patterns  have  sin  x/x  behavior, 

3)  The  quadratic  phase  error  at  the  aperture  edge  is  less  than  a  radian. 
The  actual  distance  from  the  Fresnal  field  point  to  a  point  on  the  aperture 

is  R.  Since  the  array  is  in  the  near  field  in  one  dimension  only  the  aperture 
can  be  treated  as  a  linear  array. 


?  Vi, 

Rj(Rp+x-2xR  cos  of  )  (1) 
where  R  is  the  distance  from  the  near  field  point  to  the  center  of  the  array 
and  x  is  the  distance  from  the  center  of  the  array  to  any  point  on  the  array. 


For  R°>x. 


x2-2xR  cosa  Ro-l(x2-2xRo  cos  u  )  2Ro  +  1  (x2-2xRo  cos  <*)8Ro 


p=  b  + 


2Ro2 


8Ro4 


12Ro6 


SU-2X  ftp  IQS  Q  + 

128  Ro8 


etc 


(2) 


after  combining  terms 

R=Ro-xcosb<  +  x2(l/2Ro-cos  39)  +  x8(l/2Ro2cos  w  -  1/2  Ro2  cos8)- 

2-R0  " 

x4(1/8  Ro8- 12/ 16  cos2  /Ro+80/12^  cosO'Ro3  +  . .  (3) 

Let  U  =90  -<* 

R=Ro-xsin  ©  +x2/2Rocos2€R  x8/2Ro  S<s&cos2§~x4/8Ro4 ( cos  20(l-5sin2))  (4) 
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With  x'  =  D 

x 

2 

and  u  =  kJD  sin  u  For  The  Linear  Term 

2 

Q  =  kcos3  up)2  For  The  Quadratic  Terra 
8R0 

kD3  sin  WCOS2^  For  The  Cubic  Term 

i6«;2 

and , 

o(  =  kn^  Cos3  0  ( I-5Sin3  u  )  For  the  Fourth  Order  Terra 
128R3 

1  j(ux-°x3-  x3  +  x^+...) 

E(u)  k  ^  F(x)e  dx 

-1 

equation  7  can  be  approximated  by 

l  w .. — u,.2_* 

A  J  \  li/v  A  W  A  ^ 

E(u)  k  ^  F(x)e  dx 

-1 

The  actual  far  field  Pattern  is  F.Q(u)  where, 
n1  jux 

E0(u)  =  k  Q  F(x)e  dx  =  g0(u) 

-1 

Since  Ro  is  chosen  so  that  8  are  less  than  1 
1  j  ux 

E(u)"k  \JF(x)e  ( l-j  hx3-j %  x3)dx 

-1 
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1  j  UX  n  1  jUX  1 

E(u)  ^  k  J  F(x)e  -  jpk^F(x)x2e  dx  -n^  F(x)x3 


jux 

Je  dx 


1  jux  1  jux 

g0(u)-j/?K  $  F(x)x2e  dx  -j'tf  x  j  F(x)3e  dx 

-1  -1 


Since, 


-i 


jux 

x^  F(x)e  dx 


1  dk  1  jux 

(j)^  duk  §  F(x)e  dx 

-1 


jp(2)  .  , 

■  80(“) - g0(u)-  g‘(u) 

where  g(n3  represents  the  derivative  to  the  power  the  corresponding  power 

pattern  is  P(u)  *  E(u)  E*(u) 

P(u)  | gQC ti)—  g£33(u)|2  +  .p2go(232  +  higher  order  terras 

-  g£(u)  +  'l  "(u)-  2>g0(u)g^^(ii)  +  j3*g0^)4  +  higher  order 

Using  the  assumption  that  the  ideal  far  field  Taylor  or  similar  far  field  Pattern 

can  be  approximated  by 

gQ(u)  =  a  sinu  bf (u)  cf ( u) 

u  u2  u3 

So  beyond  small  angle, 

gQ(u)  =s  a  slnu 


r<2) 


-sinu  -  2cosu  +2sinu 


(u)  a 


r- 


J 


; (  33 ( u) =  a 


l-cosu  +  3sinu  +6sinu  -6sinu' 


E(u)  ~  Fo(u3  -  P0(u) 

.,«*  u.,2  +  21  tC0SzJ  \  H £  + 

3ku2  V  u4  u2 

+sinucosua2y  -2Vu2^\ 
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Comparison  of  terms  indicates  the  2<  term  is  the  strongest.  Simplifying  further,  the 

u2 

error  envelope  is: 

^  -2^  a28P 

K(u) —  a2  and  , 


If  quadratic  phase  error  is  removed,  the  2~4  a2  becomes  the  dominant  term. 
kD~ 


f3- 


a3 


m. 


kD2 

kD2 

- KCOS2  ©  -1)  = 

_  . - 

00 

& 

o 

8R0 

Without  quadratic  phase  correction  the  dominant  error  term  is  proportional  to 
COS2  Q  so  the  distortion  is  greatest  at  ©=0.  With  the  corretion  this  error 
term  be  :oraes  proportional  to  sin2©  and  so  is  very  small  near  0=0.  As  might 
be  expected,  this  change  results  in  the  evaluation  of  the  ^  term  as  the  dominant 
error.  Examination  of  this  expression  gives  the  resulting  error  as  a  function 
of  the  array  size,  measurement  distance  and  wavelength.  This  error  is 


E(u)- 


and  is  the  absolute  resulting  error.  The  normalized  error  or  percentage  error 
is 

^  E(u) 

E(u)~ - -  2* 

go<u) 


or, 


kD3 


E(u)  = 


8R. 


4  M2 


Sin  0  cos2  0 


Sin©  cos2© 
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This  gives  the  residual  error.  If  a  is  also  the  acceptable  error 
allowed  by  this  measurement  technique  the  minimum  range  for  use  of  quadratic 
phase  only  is 

D3  1/2 

R  >  - 

.3>£ 

For  example  if  the  maximum  added  error  is  .ldb  (.023)  and  is  .16',  the 
R  30D3/2 


D 

R 

2D2/A 

8' 

1920' 

3072' 

lb' 

1920' 

3072' 

32' 

5430' 

12288' 

If  this  range  is  too  great  the  error  can  be  reduced  by  using  cubic  phase  correction 

repeated  several  times.  This  reduces  the  error  around  a  limited  angular  range.. 

The  normalized  error  is  then, 

D3 

E(u)  - (sin  ©cos2©  -Sin0  Cos2©  ) 

Ar2 

With  an  approximate  range  of  10°  and  the  third  order  correction  the  range  can  be 
reduced  by  a  factor  of  8. 


D 

R(cubic  correction) 

R(with  quadratic  only) 

8’ 

85' 

680' 

16' 

240' 

1920’ 

32' 

680' 

5430' 

Over  the  full  range  of  180°  the  phase  correction  must  be  applied  on  the 
order  of  9  times  since  the  correction  is  double  valued. 

CONCLUSION: 

The  use  of  both  quadratic  and  when  necessary  third  order  phase  correction 
permits  the  accurate  measurement  of  the  far  field  pattern  of  very  large  phased 
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arrays  at  ranges  much  less  than  the  normal  far  field  distance.  Thl3 


technique  applies  only  to  arrays  with  one  dimension  much  greater  than  the  other 
and  only  when  the  statistical  sidelobe  behavior  not  the  exact  (with  any  particular 
phase  shifter  setting)  is  required. 
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Abstract  -  A  methodology  was  developed  using  Holographic  Metrology  to  study  large 
aperture  parabolic  focusing  antennas  operating  in  the  near  field  at  1-12  GHz.  The  methodology 
included  the  development  of  an  algorithm  characterizing  the  RF  hologram.  RF  sources  and 
reflections  create  the  EMF  phase  and  amplitude  distribution  across  the  aperture,  which  are  used  to 
calculate  the  RF  hologram.  The  algorithm,  as  well,  calculates  the  radial  variation  associated  with 
intensity,  a  significant  property  when  analyzing  near  field  signals. 


A  computer  model  was  developed  to  implement  the  algorithm,  simulating  an  actual 
parabolic  focusing  antenna  with  an  adjustable  log-periodic  dipole  array  feed.  Extensive  field  tests 
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analysis  validated  the  use  of  the  algorithm,  proving  the  methodology  successful. 


I.  INTRODUCTION 

This  paper  describes  the  mathematical  model  used  to  determine  the  gain  of  very  high 
frequency,  large  aperture  antennas  used  for  electromagnetic  emission  testing.  The  model  was  used 
to  evaluate  the  gain  at  TEMPEST  distances  for  a  focusing  parabolic  reflector  over  a  frequency 
range  of  1.0  GHz  to  12.4  GHz. 

A  computer  model  was  established  to  replicate  the  mathematical  model.  The  characteristics 
of  a  specified  parabolic  antenna  with  a  log-periodic  antenna  at  its  focus  were  derived  and  integrated 
into  the  computer  model  as  a  test  of  accuracy  and  effectiveness. 

The  antenna  under  test  was  a  45-inch  diameter  Parabolic  Focusing  Antenna  with  an 
adjustable  log-periodic  dipole  array  feed.  Field  tests  were  performed  at  the  Rome  Air  Development 
Center,  Antenna  Test  Facility  in  Ipswich,  Massachusetts.  The  results  were  compared  to  those 
obtained  from  the  computer  model,  validating  the  algorithms  mathematical  model.  Comparative 
data  plots  will  be  presented  in  the  conclusion. 


II.  MODEL  DESCRIPTION 


Huygens  was  first  to  suggest  that  any  wavefront  could  be  regarded  as  an  array  of 
secondary  sources  and  he  postulated  that  the  field  of  the  secondary  source  should  be  entirely  in  the 
direction  of  the  advancing  wavefront.  In  addition,  the  Induction  Theorem  states  that  each 
secondary  source  is  a  combination  of  electric  and  magnetic  current  elements  whose  moments  are 
respectively  proportional  to  the  magnetic  and  electric  intensities  tangential  to  the  wavefront 
(Schelkunoff,  pp.  357-358).  What  this  means  is  that  once  we  know  the  electromagnetic  force 
phase  and  amplitude  distribution  across  an  aperture,  we  can  determine  by  the  radiation  equations 
what  the  three  dimensional  radiated  field  distribution  will  be.  When  this  is  accomplished  by  using 
filters  and  lasers  in  the  light  frequency  bands,  the  result  is  a  hologram.  Focusing  to  a  point  is  a 
relatively  simple  application  of  these  principles  where  the  object  is  to  obtain  a  three  dimensional 
EMF  pattern  that  converges  to  a  point  on  the  boresight  vector  line  of  an  antenna. 


Figure  1  -  Plane  aperture  of  diameter  "p’  with  a  phase  distribution  that  focuses  the 
energy  at  the  point  F. 
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III.  PARABOLA  EQUATIONS 


As  shown  in  Figure  1,  a  plane  circular  aperture  is  illuminated  by  a  point  radiator  at  F,  a 
point  on  the  boresight  WF  of  the  aperture.  The  electromagnetic  wave  arriving  at  W,  the  center  of 

the  aperture,  will  lead  in  phase  the  EM  wave  arriving  at  b,  the  point  on  the  aperture  a  distance  b-0 
from  the  center.  To  bring  into  focus  the  radiator  F  it  will  be  necessary  to  effectively  match  the 
phases  in  each  differential  portion  of  the  aperture.  In  effect,  it  will  be  focusing  on  one  element  of  a 
three  dimensional,  one  color  (frequency),  hologram.  The  phases  will  be  matched  by  means  of  a 
focusing  parabola.  (The  parabola  was  used,  even  though  it  was  not  the  optimum  aperture  for  the 
task,  because  it  was  available  to  validate  the  model). 


Figure  2  -  Section  through  the  center  of  a  focusing  parabola 


From  Figure  2  it  can  be  seen  that 

(1)  (F-X)/R  =  COS(0) 
and 

(2)  Y/R  =  SIN(9) 

But  for  a  parabola 
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(3)  F+W  =  R+W-X 


or 


(4)  R  =  F+X 

Substituting  Eq.4  into  Eq.  1  and  solving  for  X 


(5)  X=  F 


[i-cos(e)] 

[l+COS©)] 


Which  can  be  reduced  to 
(6a)  X  =  F[TAN2(0/2)] 


r-* :i i..  _ r*  £  r-  a _ a  r% 
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(6b)  R  =  F[SEC2(0/2)] 

And  then  using  Eq.2 
(6c)  Y  =  2F[TAN(0/2)] 

Thus  the  parabola  is  described  in  terms  of  F,  a  constant,  the  focal  distance,  and  the 
variable,  "0". 
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IV.  MOVING  THE  TRANSDUCER 


To  vary  the  phase  across  the  aperture,  the  transducer  is  moved  out  to  B,  the  ray  from  the 
transducer  at  B,  its  new  position,  now  makes  an  angle  "VF"  with  horizontal.  The  new  angle  ,"4/ ", 
can  be  obtained  in  terms  of  "  8"  by  the  following  expression  for  its  tangent. 

2  F  TAN  Q 

(7a)  TAN(V)  =  Y  /  (F+d-X )  _ 2 

F+d-FTAN2 

2 


So  that 

(7b)  V  =  TAN'  ’ 


2  F  TAN-®- 

A, 

F+d-F  TAN2-®— 
2 


Also  "P",  the  angle  the  ray  makes  with  the  ray  from  "0"  is  given  by 


(8)  p  =  (90°-'F)-(90°-  0)  ■  (0-  ¥) 


Figure  3  -  Phase  variation  across  the  aperture  when  then  transducer  is  moved  an  additional 
distance  d  out  from  F. 
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Using  Snell's  Equation  which  states  that  the  angle  of  incidence  equals  the  angle  of 
reflection,  the  reflected  wave,  Xb,  also  makes  an  angle  of  "P"  with  the  horizontal.  Thus 

(9)  TAN(p)  =  (Y-b)/(W-X) 

And  solving  for  b. 

(10)  b  =  F{  [TAN(e-'F)][TAN2(e/2)]+2TAN(e/2) )  -WTAN(G-'F) 

In  order  to  obtain  the  relative  phase  across  the  aperture,  it  is  necessary  to  find  the  relative 
time  delay.  In  free  space,  the  distance  delay  is  proportional  to  the  time  (or  phase)  delay.  As  a 
reference,  the  delay  will  be  related  to  the  normal  distance  delay  of  (F+D+W)  as  shown  in  Figure  4. 

The  relative  delay,  AL,  for  a  transducer  at  the  point  B  would  be 

(11)  AL  =  Li+L2-(F+d+W) 

Where 

COS(P)  WjK  -  W-Fl  AN2(Q/Z) 

L2  L2 


Or 

(12)  L2  =  W-F[TAN2(0/2)]/(COS(p)) 

Thus 

(13)  L2.  W-F1TAN2(0/2)1 

cos(e-T) 

Where  "'F"  is  shown  in  Eq.  7b. 


(14)  SIN('F)  =  Y/  Li 


So  that 


I 

(15)  Li  =  2FfTAN(0/2)l 
SIN('P) 


Solving  for  "AL" 

(16)  AL  «2ECIAN(fi/2H  +  W-FrTAN2te/2)1  -  (F+d+W) 
SIN('F)  COS(0  -'F) 


Converting  the  distance  to  radians  by,  multiplying  by  2nlX,  and  noting  the  phase 
variation  in  "AL"  as  "<!>",  we  obtain: 


(17)  2  n 

2  F  TAN(©/2)  +  W-FTAN  (0/2)  -  (F+d+W) 

X 

SIN(V)  COS(0  -NO 

V.  TRANSDUCER  EQUATIONS 


As  shown  in  Figure  4,  the  transducer  is  moved  out  to  the  point  B,  a  distance  d  from  the  fat 
field  focusing  point  F,  in  order  for  the  parabola  to  focus  at  the  required  point  C,  at  distance  D  from 
the  aperture  plane.  From  the  point  C,  the  signal  must  travel  a  distance  L  to  the  aperture  plane 

intercept  (W,b).  The  angle  the  line  L  makes  with  the  parabola  normal  is  noted  as  "4>"  so  that: 


(18)  L  =  D/(COS(<j>» 

However 

(19)  TAN(<())  =J2 

D 

Using  (10)  and  solving  for  "<J> " 

(20)  $  =  TAN-1  { [F/D]([TAN(©  -'F)]TAN2(9/2)  +  2TAN(0/2))  -  (W/D)TAN(0-vF)) ) 

Substituting  into  (18) 


(21)  L- 


COSITAN'1  [  (F/D)  ((TAN(  0  -  V  )]TAN2(  0  /  2)  +  2  TAN(  0/2) -(W/D*TAN(0  -  V)) 


Again  using  the  centerline  distance  of  travel  as  the  reference  distance,  the  distance  D  is 
s"btracted  from  L  to  get  the  relative  phase  variation  to  focus  at  the  point  C.  Converting  to  radians, 
and  calling  this  required  phase  variation  "F"  we  get: 


(22)  <*»  «  2kd 


COS(  TAN  *  [  (F/D)  ((TAN(  0  -  y)ITAI^ (  0  /  2)  +  2  TAN(  0  / 2)  -  (W/D)  TAN(9  -  f)  ] 


-  1 
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In  order  to  have  the  antenna  focus  on  the  point  C,  the  accumulative  EM  energy  at  the 
transducer  should  arrive,  with  the  minimum  phase  difference,  in  all  its  components  (Tliis  means 
that  the  energy  arriving  from  each  differential  area  on  the  aperture  should  add  to  the  energy  arriving 

from  every  other  differential  area  with  as  little  phase  difference  as  possible).  In  other  words,  O 

plus  <J>'  should  change  as  litde  as  possible  from  one  differential  area  on  the  aperture  to  the  next. 
When  this  condition  is  met  the  antenna  will  effectively  focus  at  the  required  point  C. 

VI.  TRANSDUCER  ANTENNA 

To  completely  define  the  received  vector,  the  relative  amplitude  of  each  differential  radiated 
wave  has  to  be  determined.  There  are  many  possible  transducer  antennas,  but  since  it  was 
desirable  to  field  test  the  methodology  a  log  periodic  dipole  used  with  the  antenna  to  be  tested  was 
used  in  the  calculations.  (Antenna  Research  Associates,  Inc.  antenna  LPD- 11 2/A).  Referring  to 

Figure  3,  "^"will  vary  with  Y  in  accordance  with 

(23)  TAN(T0  =  Y/(F+D-X) 

or 

(24)  ¥  =  TAN-1  [ Y/(F+D-X)] 

Thus  "  4y”  will  have  as  its  minimum 

(25)  ¥<>  =  0 

Its  maximum  value  will  be  its  value  at  the  edge  of  the  dish  (b  equal  to  p/2).  Calling  this 
value "  Yp", 

(26)  Tp  =  TAN-1  { p/[2(F+D-X)] } 

Thus  the  relative  amplitude  used  to  determine  the  relative  gain  of  the  incident  wave 

(effective  amplitude)  on  the  transducer  is  determined  from  the  transducer  pattern  of  width  2  vPp. 
Since  the  pattern  of  a  log  periodic  array  is  not  omnidirectional  but  instead  resembles  a  broadside 

dipole  pattern,  it  was  decided  to  use  a  cosine  taper  and  define  each  frequency's  taper  by  the 
available  manufacturer's  antenna  pattern. 
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Actually,  the  E  and  H  patterns  were  not  similar  nor  exactly  the  cosine  shape,  but  the 
differences  were  not  significant  enough  to  warrant  the  assumption  of  a  complex  taper.  The  major 
differences  were  above  10  Ghz,  but  even  there  the  differences  were  less  than  3  to  5  db. 

Defining  the  Cosine  taper  in  terms  of  the  radial  distance  b: 

(27)  assumed  Taper  -  COS(Jtb/Q) 


Where  Q  is  a  constant  for  each  taper  and  has  to  be  determined  from  the  log  periodic  antenna 
patterns.  At  b  equal  to  zero  the  cosine  is  equal  to  one.  At  the  edge  of  the  dish  were  b  is  equal  to 
r/2: 

(28)  COS(np/2Q)  =  Power  at  To  (From  the  pattern) 

Power  at  'Po  (From  die  pattern) 

Thus  the  constant  Q  for  a  particular  pattern  is  given  by: 


(29)  Q 


rc  p 


1 


2 


COS'1 


POWER  AT 
POWER  AT  V0 


Substituting  into  equation  27: 


(30) 


r  *1 

(—  * 

- 

TAPER  *  COS 

_2b  COS 

POWER  AT  Hfi 

P 

[POWER  AT  % 

VII.  WEIGHTING  FUNCTIONS 


The  aperture  distribution  of  assumed  reradiators  is  circularly  symmetrical  such  that  a 
differential  zone  plate  can  be  used  for  the  differential  reradiator.  This  yields  a  power  weighting 
function,  PWF,  where: 


(31)  PWF  =  2* 


COS  2  b“ 

f-  -1 

COS 

POWER  AT  Hfj 

— 

P 

POWER  AT 

As  the  transducer  is  moved  farther  away  from  the  dish,  the  power  also  decreases  with 
distance.  (Refer  to  Figure  2).  Let: 


(32)  5e  =  F  +  W 

A  r\A  • 

i  U1U> 

(33)  5y  =  F  +  W  +d 

Thus  the  power  density  correction  factor,  PDCF  is  given  by: 

(34)  PDCF  =  6e2/ V 

The  field  intensity  is  also  dependent  on  each  differential  reradiator  in  the  aperture  as  well  as 
the  distance  from  each  reradiator  to  the  transducer.  From  Figure  4: 

(35)  L  =  (D2  +  b2)i/2 

And  to  account  for  the  cosine  taper  that  is  being  used,  cos',vP"  is  incorporated  where: 

(36)  COS('P)  =  D/L  =  D/  ((D2  +  b2)i/2) 
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Taking  the  product  to  obtain  the  correction  factor: 


(37)  COS(40  =  D _  ■ _ 1 _ 

(D2+b2)l/2  (D2  +  b2)l/2 


Yielding: 


(38)  Pattern  Intensity  correction  factor  =  D/(D2+b2) 


VIII.  INTEGRATION  EQUATION 


The  relative  differential  intensity  assumed  to  be  reradiated  from  a  differential  zone  plate,  of 

width  db,  at  the  radial  distance  b  from  the  center  of  the  aperture  is  given  by  dEy  where  P0  is  a 
constant  determined  by  the  feed  antenna  and  the  distances  involved.  Its  value  is  unimportant  in  our 

case,  because  it  will  be  canceled  out  when  determining  the  relative  antenna  factor. 


1  D 

52e 

“cos 

2  b  P0  COS  ^  [POWER  AT  Nf) 

— 

- 

1  1  nktf 

_S2v 

P  [POWER  AT  % 

CD 

Integration  can  now  be  performed  for  the  complete  differential  vector  using  the  intensity 

from  equation  39  and  the  phase  angle,  from  equations  17  and  22.  Thus  the  field  intensity  at 
C,  Ec,  is  given  by 

ft-, 

(40)  Ec  a^dEW  fl>+<E' 

— * 

In  effect  the  phase  is  determined  by  the  relative  distances  the  signal  travels  from  the  point  of 
origin  to  the  transducer.  The  field  intensity  in  the  far  field,  Ef  is  obtained  by  assuming  C  in 
the  far  field.  Similarly,  the  field  intensity  at  one  meter,  Ea,  is  obtained  by  assuming  C  to  be  one 
meter. 
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The  change  in  antenna  factor  is  the  ratio  of  the  intensity  received  in  the  far  field  compared  to 
the  intensity  received  when  the  antenna  is  used  at  one  meter.  Thus,  the  change  in  antenna  factor, 
(CAF),  is  given  by: 

(41)  CAF  =  20  logio(Ef/En) 

Ef  is  the  relative  intensity  in  the  far  field  and  En  is  the  relative  intensity  at  one  meter. 

IX.  CONCLUSION 

A  final  conclusion  to  tills  study  is  that  the  mathematical  model  for  large  aperture  antennas, 
reconciled  with  the  experimental  work  conducted  on  the  antenna  under  test,  demonstrates  that  the 
methodology  is  valid  for  determining  the  intensity  for  a  given  radial  distance. 

Figure  5  through  figure  9  are  plots  comparing  the  field  test  data  to  the  mathematical  model 
output  for  a  similar  test  scenario.  The  coordinates  of  the  plots  are  gain  (db)  versus  the  point  of 
focus  (M).  The  results  show  remarkable  resemblance  given  the  fact  that  the  model  contained 
various  approximations. 

o 
-2 


-« 

•t 
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-12 

Q*ln(dB)  "14 
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-28 
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Figure  5  -  Comparison  of  Test  vs.  Model  at  2.6  GHz  (S-Band) 
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Figure  6  -  Comparison  of  Test  vs.  Model  at  4.0  GHz  (S-Band) 
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Figure  7  -  Comparison  of  Test  vs.  Model  at  4.0  GHz  (C-Band) 
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This  paper  demonstrated  that  Holographic  Metrology  can  be  used  to  analyze  the  near  field 
operator  of  large  aperture  antennas.  Computer  programs  (see  appendices  A  and  B)  for  algorithms 
based  upon  this  metrology  can  be  formulated  that  describe  the  three  dimensional  (  azimuth  and 
elevation )  distribution  of  sources  and  reflections  that  generate  a  specific  aperture  R.F.  Hologram. 
The  distribution  obtained  is  unique. 
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! Parabola  Program  "TASRJO"  -  CORRECTION  FACTOR 
! UPDATED  OCTOBER  27,  1987 

[THIS  PROGRAM  VARIES  THE  LOCATION  OF  THE  FEED 
! "SMD"  FOR  A  FIXED  POINT  OF  FOCUS  "D"  -  THE 
[OUTPUT  FILE  GENERATED  "ATE180"  IS  FORMATTED  AS 
1 POWER  INTENSITY  VERSUS  THE  POINT  OF  FOCUS  -  THIS 
[DATA  FILE  IS  THEN  INSTALLED  INTO  THE  PLOT  FILE 
1 "PLSMUZ" 

1 

[DEFINITION  OF  VARIABLES 
! 

■FREQ  -  FREQUENCY 

IRTOD  -  RADIANS  TO  DEGREES 

!DTOR  =  DEGREES  TO  RADIANS 

!TFI  =  TW°*PI 

!C  =  SPEED  OF  LIGHT 

!TPIL  =  TWO  -  PI  DIVIDED  BY  LAMBDA 

I SMD  =  FEED  LOCATION  WITH  RESPECT  TO  THE  FOF 

ID  ■=  POINT  OF  FOCUS 

!F  ■=  FOCAL  POINT 

[RHO  «  Y  -  INTERCEPT  OF  T)IE  DISH 

!W  =  X  -  INTERCEPT  OF  THE  DISH 

[CASE  *  FILE  CONTAINING  FEED  DISTRIBUTION 

! ATE180  =  OUTPUT  FILE  CONTAINING  INTENSITY 

!  VALUES  FOR  VARIOUS 

i  POINTS  OF  FOCUS  AT  A  GIVEN  FREQUENCY 

[FACTOR  =  THE  POWER  INTENSITY 
!POF  =  POINT  OF  FOCUS 

j 

CREATE  BDAT  "TV346", 200, 8 
ASSIGN  9Path2  TO  '"10/34  6" 

ASSIGN  8Path3  TO  "CASE3" 

C  -  300000000 
Tpi  =  2*PI 
Rtod  ■  180/PI 
Dtor  -  PI/180 
Freq  -  1.2E+10 
F  =  .50055 
W  -  6.61/39.37 
RHO  =  22.75/39.37 
Lambda  =  C/Freq 
Tpil  =  Tpi/Lambda 
DIM  Rhopwr  (60) 

Smd  =  .346 
FOR  X  «  1  to  60 

ENTER  @Path3,  X;  Rhopwr  (X) 

Mnvni  »» 

IILjJVI 

D  =  .6787 
FOR  T  =  1  TO  200 
Z  =  F +Smd 
Rel  =  0 
Imag  •=  0 
Delt  »  F+W 
Delp  »  F+W+Smd 
Sqdt  =■  Delt*Delt 
Sqdp  =  Delp*Delp 
Delaq  *  Sqdt /Sqdp 
Ref  =  Tpil* (F+Smd+W+D) 

N  »  0 
Taper  =  0 
Prime  =*  0 
Reffaz  =  0 
Capphi  =  0 
Tpwrth  =  0 
Relsq  =  0 
Imagsq  =  0 
Sumsq  =  0 
Magpsi  =  0 
Factor  =  0 


Db  *  .05 
r 

[ENTER  LOOP  TO  CALCULATE  THE  POWER  INTENSITY 
[OVER  THE  ANGLE  "THETA"  AT  A  GIVEN  POF 

l 

FOR  Theta  **  1.0  TO  60.0 
N  »  Nfl 

Hlfeta  «  ( ,5*Theta) *Dtor 
Tanhlf  **  TAN  (Hlfeta) 

Taraq  “  Tanhlf *2 
R  =  F /  (COS  (Hlfeta)  "2) 

X  -  F*Tanaq 
Y  -  2*r*Tanhlf 
Pal  ■=  ATN  ( ( Y )  /  (Z-X)  ) 

Pal  =  Pal/Dtor 
Beta  =  Theta-Pal 

Beta  =»  Beta*Dtor 
Psi  =  Psi*Dtor 

B  =  F*  (TAN  (Beta)  *Tan.sq+2  ‘Tanhlf)  -W*TAN  (Beta) 
Wft=*  W-X 

Capphi  -  Tpil*  ((Y)/SIN(Pai)  +  (Wft)  /COS  (Beta)  ) 
Baq  =  B*B 
Daq  ■  D*D 
Bplsa  =  Baq+Dsq 
Rdcomp  *  D/Bpl-id 
Sqdb  =  SQR(Bplsd) 

Pprime  =  Tpil* (Sqdb) 

Reffaz  =  (Pprime+Capphi) -Ref 
Pwr  =  Rhopwr (Theta) 

Taper  =  Tpi*B*Db*COS ( (2*B/Rho) *Acoa (Pwr) ) 
Pv?rtht  —  SQP.  (Tspsi**Ds  J.sq)  *  Rdcomp 
Tpwrth  =  Tpwrth+Pwrtht 
Intpsi  -  SQR (Taper*Delsq) *Rdcomp 

Rel  =  Rel+Intpai*COS (Reffaz) 

Imag  *  Imag+Intpai*SIN (Reffaz) 

NEXT  THETA 

; 

i  EXIT  "THETA  LOOF"  AND  DETERMINE  "FACTOR" 

i 

Relsq  ■*  Rel*Rel 
Imagaq  =  Imag* Imag 
Sumsq  =  Relaq* Imagaq 
Magpsi  *  SQR (Sumsq) 

Factor  =  Magpsi /Tpwrth 
OUTPUT  @Pnth2,T; Factor 
D  =  D+.005 
NEXT  T 
END 
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AV  VEND  XX  B 


1THIS  PROGRAM  PLOTS  THE  OUTPUT  FILE  GENERATED 
!  FROM  "TASRJO”  -  THE  DATA  IS  PLOTTED  AS 
I  POWER  INTENSITY  VESUS  POINT  OF  FOUCUS 
ASSIGN  @Path  TO  "ONEOOO" 

INPUT  "ENTER  #  OF  DATA  POINTS",  Sata 

INPUT  "ENTER  X  -  AXIS  LENGTH  IN  METERS’’,  Scat 

DIM  Dumy  (200) 

FOP.  X  -  J.  TO  Sata 
ENTER  @Path, X; Dumy (X) 

NEXT  X 

PLOTTER  US  705, "HPGL” 

PEN  5 

VIEWPORT  10,110,10,95 

FRAME 

AXES  4,10 

CSIZE  4, .27 

G  »  H 

Test  ■  . 1 

MOVE  0,0 

FOR  X  -  1  TO  5 

MOVE  1, G 

LABEL  TEST 

TEST  *  TEST+.2 

G  »  G+20 

NEXT  X 

Testl  *  .4545 
Testot  «•  Testl 
G  -  9 

FOR  Y  *  1  TO  11 

MOVE  G,  2 

LABEL  Testot 

Testot  «  Testot+Testl 

G  =  G+11.5 

NEXT  Y 

PEN  5 

GRID  12,10 
PEN  5 

MOVE  70,80 
CSIZE  4,  .3 
LABEL  "INTENSITY" 

MOVE  70,75 

LABEL  "AMP  VS.  METERS" 

MOVE  70,70 

LABEL  "FREQUENCY  «  1.0  GHZ" 

MOVE  70,65 

*  mm  neAT r'/’V'nr-  t  r*  AnTTMV'TPn  PAO  TUP  XT'  T>  H> 
linDLjlj  C  vill  1  L/i.  L  i.  VI  iilii.ULil'  t  SXi  %  i  ‘  ^  'VI' 

FIELD" 

MOVE  0,0 

FOR  T*1  TO  Sata 

Dumx  =  T* (133/Scat) 

Dumx  ■  Dumx*  (Scat/Sata) 

Dumy(T)  ■  Dumy  (T)  *100 
DRAW  Dumx, DUMY  (T) 

NEXT  T 
PEN  0 
END 
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ABSTRACT 


Reflector  antennas  on  earth  orbiting  spacecrafts  generally  cannot  be  desciibed  analytically.  The 
reflector  surface  is  subjected  to  large  temperature  fluctuation  and  gradients,  and  is  thus  warped  from  its 
true  geometrical  shape.  Aside  from  distortion  by  thermal  stresses,  reflector  surfaces  are  often 
purposely  shaped  to  minimize  phase  aberrations  and  scanning  losses.  To  analyze  distorted  reflector 
antennas  defined  bv  discrete  surface  points,  a  numerical  technique  must  be  applied  to  compute  an 
interpolator  surface  passing  through  a  grid  of  discrete  points.  In  this  paper,  the  distorted  reflector 
surface  points  are  approximated  by  two  analytical  components;  an  undislortcd  surface  component  and  a 
surface  error  component.  The  undislorled  surface  component  is  a  best  fit  paraboloid  polynomial  for  the 
given  set  of  points  and  the  surface  error  component  is  a  Fourier  series  expansion  of  the  deviation  of  the 
actual  surface  points,  from  the  best  fit  paraboloid.  By  applying  the  numerical  technique  to  approximate 
the  surface  normals  of  the  distorted  reflector  surface,  the  induced  surface  current  can  be  obtained  using 
physical  optics  technique.  These  surface  currents  are  integrated  to  find  the  far  field  radiation  pattern. 


ANALYTICAL  APPROXIMATION  OF  A  DISTORTED  REFLECTOR 
SURFACE  DEFINED  BY  A  DISCRETE  SET  OF  POINTS 


Roberto  J.  Acosta  and  Afroz  A.  Zaman 


National  Aeronautics  and  Space  Administration 
Lewis  Research  Center 
Cleveland,  Ohio  44135 


L  Litrad  action 

Reflector  antennas  in  a  space  environment  are  subjected  to  continuous  variation  in  temperature 
distribution,  and  arc  thus  distorted  from  its  true  geometrical  shape  (typically  parabolic,  hyperbolic, 
elliptical,  etc.).  The  distorted  reflector  surface  has  in  general  a  very  complicated  shape  and  hence  can 
not  be  represented  with  an  analytical  expression.  The  analysis  of  a  distorted  reflector  antenna  defined 
by  a  set  of  discrete  surface  points  require  the  use  of  numerical  techniques.  Many  numerical  techniques 
for  analyzing  the  performance  of  reflector  antennas  defined  by  a  set  of  points  have  been  extensively 
reported  in  the  open  literature  (refs.  1-4).  The  most  popular  of  these  techniques  represent  the  reflector 
surface  either  globally  or  locally  by  using  polynomial  splines  (refs  5,6).  The  polynomial  splines  and 
other  techniques  require  the  surface  points  to  be  orderly  labeled.  This  is  an  undesired  characteristic 
because  it  can  result  in  a  nonunique  approximation  for  the  desired  reflector  surface.  The  order  of  the 
spline  polynomial  necessary  to  best  fit  the  reflector  surface  points  is  in  general  unknown  a  priori,  and 
therefore  it  involves  a  trial  and  error  procedure  for  checking  the  accuracy  of  the  interpolation. 

The  distorted  reflector  surface  points  can  be  best  approximated  by  two  analytical  components,  an 
undislorlcd  surface  component  end  a  surface  error  component.  The  undistorted  surface  component  is 
a  best  fit  paraboloid  polynomial  for  the  given  set  of  points  and  the  surface  error  component  is  the 
deviation  of  the  actual  surface  points  from  the  best  fit  paraboloid.  This  residual  error  component  is 
then  described  with  a  sinusoidal  Fourier  series  expansion.  This  approximation  technique  is  insensitive 
to  the  labeling  of  the  reflector  surface  points  and  can  describe  the  surface  errors  very  accurately. 
Similar  to  a  lime  signal  the  spatial  spectrum  of  the  surface  error  component  is  unique  to  the  reflector 
under  consideration.  Therefore  spatial  spectra  can  be  utilized  as  a  performance  index  for  comparing 
distortion  profiles  in  reflector  antennas. 
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The  best  analytical  representation  of  the  distorted  reflector  antenna  surface  that  uniquely  identifies 
the  surface  errors  can  be  obtained  as  follows:  the  reflector  surface  points  are  separated  into  two 
components,  a  best  fit  paraboloid  component  and  a  sinusoidal  Fourier  series  expansion  of  the  residual 
(ref.  7).  Figure  1  illustrates  a  conceptual  layout  of  the  problem  under  consideration. 

In  analyzing  large  reflector  antenna  performance  it  is  necessary  to  accurately  characterize  the 
reflector  surface  points.  Any  deviation  from  its  ideal  geometry  causes  the  antenna  performance  to 
degrade.  The  surface  error  component  provides  an  independent  performance  index  against  which 
distorted  reflector  antennas  can  be  compared.  The  surface  error  component  by  definition  is  a 
sinusoidal  Fourier  series  expansion  of  the  difference  between  the  actual  reflector  surface  points  and  the 
best  fit  paraboloid  geometry.  Information  such  as  root  mean  square  value  (RMS),  peak  surface  error 
and  two-dimensional  distortion  profile  can  be  calculated  from  the  surface  error  component.  Fourier 
coefficients  in  the  series  expansion  represent  the  spatial  spectrum  that  uniquely  identifies  the  distorted 
reflector  under  study.  The  best  fit  paraboloid  surface  represents  the  reflector  antenna  surface  in  an 
average  sense.  When  the  surface  error  is  zero  the  best  fit  paraboloid  surface  reduces  to  the  ideal  or  the 
design  surface  geometry. 

3.  tjjjmehcal  flcsults.andlJisciissimi 

Figure  2  represents  a  block  diagram  of  a  computer  simulation  of  the  above  stated  problem.  The 
reflector  surface  points  are  usually  obtained  from  a  holographic,  photogramctric  or  any  other  surface 
detection  technique.  The  computer  algorithm  was  tested  with  known  distortion  profiles  and  TRASYS- 
SYNDA-NASTRAN  (Ref.l)  simulated  thermal  distortions  superimposed  into  a  reflector  anlenna 
geometry. 

The  distortion  profiles  considered  for  the  simulation  arc  described  in  Table  1  end  the  reflector 
geometry  illustrated  in  figure  3.  Case  A  in  Table  1  describes  a  small  (X/20  or  less)  distortion  profile  and 
case  B  illustrates  a  large  (in  the  order  of  several  wavelengths)  distortion  profile.  The  frequency  used  for 
the  simulation  was  10  GHz.  A  set  of  equally  spaced  data  points  (100,100)  were  needed  to  generate  the 
distorted  reflector  surface.  The  estimated  surface  errors  and  best  fit  paraboloid  surface  geometry  are 
described  in  Table  2. 

Figures  4(a)  and  (b)  show  the  estimated  surface  profile  corresponding  to  case  A  and  case  B  and  their 
respective  far  field  radiation  patterns.  The  results  presented  in  table  2  are  in  good  agreement  with  the 
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input  distortion  by  Fourier  coefficients  in  Table  2.  It  was  found  that  the  Fourier  matrix  had  higher 
order  coefficients  with  non  zero  values  but  they  were  one  order  of  magnitude  less  than  the  lowest 
amplitude  of  the  Fourier  coefficient  presented  in  Table  2. 

The  technique  was  compared  with  results  obtained  by  using  a  spline  polynomial  fit  for  approximating 
the  distorted  reflector  surface  points  (ref .  1)  In  this  case  the  distorted  reflector  surface  points  were 
obtained  by  simulating  thermal  deformation  with  a  TRASYS-SYNDA-NASTRAN  computer 
programs.  In  brief,  TRASYS  and  SINDA  are  used  to  charcterize  the  in-orbit  thermal  environment; 
NASTRAN  calculates  the  thermally  induced  mechanical  distortions.  Figure  5  shows  the  temperature 
distribution  on  the  reflector  antenna  surface  for  the  case  under  consideration.  The  neflector  geometry 
input  to  the  thermal  programs  is  presented  in  fig.  6.  The  frequency  considered  was  28.75  GHz  and 
(100,100)  surface  points  were  used  for  the  analysis.  The  far  field  radiation  pattern  corresponding  k;  the 
thermal  simulation  case  is  presented  in  figure  7.  The  continuous  line  pattern  corresponds  to  the 
polynomial  spline  algorithm  and  dotted  line  pattern  corresponds  to  the  best  fit  paraboloid  and  Fourier 
scries  expansion  approximation.  There  arc  no  major  differences  between  the  beam  direction  and 
sidclobe  levels,  indicating  a  good  agreement  between  the  two  techniques.  The  best  fit  paraboloid  and 
Fourier  series  algorithm  was  very  slow:  3  hr.  of  c.p.u.  time  in  an  IBM  370  computer.  The  spline 
polynomial  algorithm  lakes  about  1/2  hr.  on  a  CRAY  XMP  computer.  The  long  compulation  can  be 
justified  as  a  trade-off  to  obtain  vital  information  about  the  distorted  surface  characteristics.  These  arc 
the  amplitude  spectra  of  the  surface  error,  the  distorted  surface  profile,  RMS  value  and  the  largest 
deviation  cn  the  reflector  antenna  surface.  These  are  not  directly  available  from  any  of  the  other 
existing  techniques. 

4.  Concluding  Rswarks 

One  advantage  of  the  developed  technique  is  that  it  can  be  easily  implemented  to  any  existing 
reflector  antenna  secondary  pattern  computational  method.  It  can  easily  be  extended  to  nonparabolic 
reflectors  surfaces  (spherical,  planar,  hyperbolic,  elliptical,  etc.)  by  modifying  the  least  square 
polynomial  approximation.  In  applications  involving  the  fabrication  and  design  of  precision  reflector 
antennas  the  technique  can  be  used  as  a  computer  aided  tool.  Information  such  us  the  average  focal 
length,  a  root  mean  square  of  the  surface  error ,  surface  error  profile  and  the  amplitude  spectra  of  the 
reflector  antenna  under  consideration  are  "easily  accessible  from  the  algorithm.  One  draw  back  of  the 
algorithm  is  that  it  is  slow.  This  can  be  improved  bv  using  a  faster  computer  such  as  a  CRAY  and 
optimizing  the  computer  codes. 
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1  Distorted  reflector  surface  points  separated  into  an  undislorted  surface  component  and  a  surface 
error  component. 
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2  Computer  implementation  for  obtaining  an  analytical  representation  for  a  reflector  surface  defined 
by  t.  set  of  points. 
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3  Offset  parabolic  reflector  geometry 
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TABLE  1.  -  DISTORTION  PROFILE  PARAMETERS 
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TABLE  2  -  ESTIMATED  DISTORTION  PROFILE  PARAMETERS 
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Cett  A  3.290109  ft 
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METERS 


CONTOUR  CONTOUR 
LABEL  VALUE 


A 

B 

C 

D 

E 

F 

6 

H 

I 

J 

K 

L 

n 

N 

o 


0.008000 
0.004000 
0.002000 
0.001000 
0.000500 
0.000250 
0.000125 
0.000000 
-0.000125 
-0.000250 
-0.000500 
-0.001000 
-0.002000 
-0.004000 
-0  008000 


DEG 


<*>  Cat*  A  (suH  PUtortlon  profile). 


4  E-plane  radiation  pattern  and  their  estimated  distortion  profile. 
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ABSTRACT 

This  paper  will  describe  the  Adaptive  Spatial 
Signal  Processing  Test  Facility  and  some  of  the 
experimental  results  obtained  using  the  facility.  The 
test  facility  consists  of  an  anechoic  chamber,  antenna 
pattern  measurement  equipment,  and  the  Flexible 
Adaptive  Spatial  Signal  Processor  (FASSP). 

Methods  of  implementing  antenna  arrays  and 
spatially  distributed  RF  si gnal /i nterf erence  scenarios 
using  the  anechoic  chamber  are  described.  The  use  of 
antenna  pattern  measurement  equipment  to  aid  in  the 
configuration  of  antenna  arrays  and  positioning  of 
si gnal /i nterference  sources  is  also  presented. 

The  capabilities  of  the  FASSP  are  described.  The 
different  hardware  modes  of  the  FASSP  and  the  methods 
of  adaptive  beamforming  and  nulling  are  described. 
Data  from  different  arrays,  scenarios,  adaptive 


algorithms  and  FASSP  hardware  configurations  is 
presented , 

The  test  bed  has  recently  been  enhanced  by  the 
addition  of  a  systolic  array  processor.  The  function 
and  capabilities  of  this  adaptive  processor  wi  be 
descri bed . 

An  effort  to  upgrade  the  test  facility  by 
employing  optical  fiber  transmission  lines  is  also 
discussed  in  the  paper. 

1.  INTRODUCTION 

In  an  effort  to  employ  the  interference 
suppression  capabilities  of  adaptive  antenna 
technology  in  military  communications  systems,  many 
DoD  labs  have  funded  a  large  number  of  research  and 
exploratory  development  efforts  toward  this  end.  For 
many  applications  of  this  technology,  the  specific 
design  and  development  of  an  antenna  array,  signal 
processing  hardware  system,  and  algorithm  to  meet  a 
specific  user's  needs  is  required.  Although  the 
interference  suppression  capability  of  such  systems  is 
impressive,  many  applications  have  not  been 
implemented  because  they  were  not  cost  effective. 

Over  the  years,  RADC  has  attempted  to  reduce  the 
development  costs  of  adaptive  antenna  systems  by 


various  means.  In  1978,  RADC  purchased  and  Installed 
an  anechoic  chamber  and  an  antenna  analyzer  system. 
Significant  cost  savings  on  a  number  of  adaptive 
antenna  efforts  was  realized  using  this  test  facility. 
For  the  first  time  at  RADC,  antenna  arrays  could  be 
designed  and  tested  in  a  cost  effective  manner  using 
the  test  facility. 

The  Flexible  Adaptive  Spatial  Signal  Processor 
( FASSP )  was  installed  and  integrated  with  the  anechoic 
chamber  and  antenna  analyzer  In  1985,  greatly 
enhancing  the  test  facility.  The  flexibility  of  the 
FASSP  hardware  and  the  ability  to  program  algorithms 
made  the  complete  design  and  testing  of  adaptive 

antenna  systems  relativity  inexpensive. 

The  Systolic  Processor  Array  for  Radar  and 
Communications  (SPARC)  was  delivered  and  Integrated 
into  the  adaptive  test  facility  in  May  of  1988.  The 

impressive  computational  power  of  this  digital 
processor  performs  a  specific  algorithm  about  400 

times  faster  than  previously  developed  all-digital 
adaptive  processors. 

RADC  is  presently  looking  at  the  ways  In  which 

recent  developments  in  photonics  can  enhance  the 
adaptive  test  facility.  A  few  efforts  are  planned  to 
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explore  optical  iinks,  optical  beamforming  and  optical 
signal  processors. 

2.  ANECHOIC  CHAMBER  AND  ANTENNA  PATTERN  ANALYZER 

The  anechoic  chamber  was  installed  at  RADC  in  the 
RF  Lab  in  1978.  The  chamber  is  rectangular  in  shape 
and  can  be  used  over  a  wide  range  of  frequencies  from 
300  MHz  to  10  GHz  and  above.  The  receive  end  of  the 
chamber  has  a  pedestal  that  can  be  rotated  in  azimuth 
or  elevation  either  manually  or  under  computer 
control.  Antennas  are  mounted  on  the  pedestal  in  the 
center  of  a  6  ft.  diameter  spherical  quiet  zone.  The 
transmit  end  of  the  chamber  is  adjacent  to  an 
equipment  room  where  interference  sources  and 
communications  systems  can  be  installed  to  implement 
many  different  test  configurations. 

The  Scientific  Atlanta.  measurement  system 
provides  a  means  of  determining  antenna/antenna-array 
patterns.  The  antenna  analyzer  system  can  be  operated 
manually  or  automatically.  For  adaptive  antenna  array 
tests  it  is  necessary  to  locate 
interference/communications  sources  at  precise 
positions  with  respect  to  the  array.  A  computer 
program  can  be  used  to  predict  the  expected  pattern  of 
a  given  array  based  on  the  number  of  elements  and 


geometry  of  the  array.  A  laser  can  be  mounted  on  the 
pedestal  and  used  to  precisely  locate  a  broadside  CW 


source . 

Then 

the  array 

can  be 

installed  on 

the 

pedestal 

and  an 

actual 

pattern  can 

be  obtained  using 

the  CW 

source  and  the 

antenna  analyzer  system. 

The 

actual 

pattern 

and  the 

predicted  pattern  can 

be 

compared 

.  Once 

the 

two 

patterns 

agree,  the  actual 

pattern 

can  be 

used 

to 

locate 

other  sources 

for 

adaptl ve 

array  testing 

( see 

figures 

1  and  2). 

1 


FIGURE  1  -  ACTUAL  ANTENNA  PATTERN 

(Four  element  linear  horn  array) 
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FIGURE  2  -  COMPUTER  GENERATED  ANTENNA  PATTERN 
(Four  element  linear  horn  array) 


3.  FLEXIBLE  ADAPTIVE  SPATIAL  SIGNAL  PROCESSOR 

Syracuse  Research  Corporation  (SRC)  installed  and 
tested  the  FASSP  in  the  RADC  adaptive  test  bed  in 
April  of  1984.  Essentially,  the  FASSP  is  a 
hardware/software  reconf  i  gurabl e/programmabl  e  analog 
and  digital  adaptive  beamf ormi ng/i nterf erence-nul 1 i ng 
processor.  The  FASSP  has  12  SHF  receivers  which  down- 
convert  to  an  IF  of  540  MHz.  An  IF  interconnect  panel 
allows  12  IF  signals  to  be  connected  to  any  of  the  12 
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analog  and/or  digital  weighting  modules.  The  baseband 
receivers  down-convert  the  26  real  (13  complex)  IF 

channels  to  baseband  and,  using  eight  bit  sign 
magnitude  A/D  conversion,  produce  26  real  digital 

sampled  data  streams.  All  the  adaptive  data 

processing  is  done  in  an  AP120B  floating  point 
processor  and  the  digital  weights  are  converted  to 

analog,  using  D/A  conversion,  for  application  to  the 
analog  weighting  modules  or  directly  applied  to  the 
digital  weighting  modules.  The  whole  process  is  under 
the  control  of  an  PUPil/44  DEC  computer  (see  figure 
3). 


FASSP  FUNCTIONAL  DIAGRAM 


FIGURE  3 


FASSP  SPECIFICATIONS 


1  Phv"C®l/Mfch*niC*l  Laboratory  Environment.  Commert.i,il 
Cabinets.  Semi -portable 

2.  Operating  System  Software -Based.  Use-  interactive 

3.  Adaptive  Algorithms  Closed- Loon  Gradient.  Conjugate  Gradient. 

Guided  Random  Search.  Oi>Jn  Loon  Sample  Matrix  Inversion. 

4.  Number  Channels  1?  Adaptively -Wnqhted  Channels,  Derived 

From  6  RF  Input  Channels.  Via  Manually  Patched  Delay  Lines 

5  Operating  Frennency  Band  7  0  to  8  5  GH/ 

6.  Input  Dynamic  Range  -90  dBm  to  -70  dBm  Total 

40  dB  Instantaneous 

7.  Adaptive  Nuliinq  Performance  -  Mode  1 

4C  dB  Over  40  To  liO  MHj  Bandwidth  iCIosed  •  Loop) 

20  dB  Over  40  To  50  MHr  Bandwidth  (Open  Loop) 

8  Adaptive  Nulling  Performance  -  Mode  2 

40  dB  Over  5  MH/  Bandwidth  (10  MM/  Gn.iD 

9  Convergence  Speed  Suitable  For  Operation  From  A  Scanning 

Antenna  PedrsMl.  At  ?  To  10  RPM 


FIGURE  4 


FASSP  ALGORITHM  PERFORMANCE 
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FIGURE  5 
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For  a  list  of  various  performance  specifications  and 
available  algorithms  see  the  appropriate  charts  (see 
figure  4  and  5).  The  FASSP  is  capable  of  producing 
charts  and  graphs  of  antenna  patterns,  beamformer 
residue  output  power  and  Eigen  vector/value  analysis 
of  the  array  correlation  matrix  (see  figure  6). 


Cf  H  MUCH  HJ  MM/ 


FIGURE  6 


4.  SYSTOLIC  PROCESSOR  ARRAY  FOR  RADAR  AND  COMM 

The  Hazeltine  Corp.  delivered  to  RADC  the 
Systolic  Processor  Array  for  Radar  and  Communication 
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In  the  first  quarter  of  1988.  Integration  into  the 
adaptive  test  facility  and  acceptance  testing  was 
completed  In  the  second  quarter  of  1988  (see  figure 
7).  SPARC  Is  controlled,  via  an  IEEE  -  488  bus, 

through  Its  Interface  by  the  same  computer  as  the 
FASSP. 


FIGURE  7  -  SPARC 


The  processor  is  a  linear  systolic  array  which 
solves  an  M  x  12  system  of  equations  by  performing  QR 
-factorization  using  Givens  rotations.*  The  processor 
operates  directly  on  the  data  which  flows,  in 
parallel,  from  the  12  complex  channels  of  the  FASSP 
through  the  SPARC  interface  to  the  SPARC  processor. 
The  SPARC  solves  the  system  of  equations  and  provides 
a  solution  every  50  microseconds.  The  processor  does 
all  computation  in  3j  bit  IEEE  standard  arithmetic  but 
supplies  the  solutions/weights  in  8  or  12  bit  integer 
format. 

Many  of  the  SPARC  parameters  are  programmabl e . 
Some  of  the  parameters  that  can  be  selected  are  data 
source,  batch  or  iterative  mode,  pilot  vector,  active 
elements  and  scaling  for  the  integer  weight  output. 

Using  data  from  the  12  element  array  implemented 
in  the  anechcic  chamber,  the  SPARC  in  concert  with  the 
FASSP  can  remove  11  spatially  distributed  interference 
sources  from  the  residue  output  of  the  FASSP  digital 
beamf ormer . 

5.  OPTICAL  TECHNOLOGY  APPLICATIONS 

One  of  the  reasons  for  upgrading  the  adaptive 
test  facility  by  employing  fiber  optical  transmission 
lines  is  to  preserve  the  performance  of  the  anechoic 


chamber.  Installation  of  the  fiber  optical  cable  will 
reduce  the  overall  amount  and  size  of  existing  cable 
and  allow  the  removal  of  equipment  now  employed  within 
the  anechoic  chamber  to  perform  phased  array  antenna 
measurements.  The  FASSP  receiver  assembly  wi  also 
be  removed  from  the  chamber. 

There  are  some  impressive  advantages  to  fiber 
optical  systems.  Fiber  optics  offers  very  low, 
constant  attenuation  at  all  frequencies,  bandwidths 
that  are  extremely  wide,  impervious  to  EMI  and  RFI, 
and  the  reduced  size  and  weight  of  components  to 
perform  the  same  tasks. 

A  four  element  array  will  be  constructed  in  the 
normal  fashion  to  provide  a  means  of  comparison  to  a 
four  element  optically  linked  array  (see  figure  8). 
The  components  necessary  for  the  four  element  link 
have  been  purchased  and  are  being  evaluated.  The 
components  will  be  assembled  into  links  and  time  and 
phased  matched  as  required  by  the  beamforming/nulling 
adaptive  process.  The  optically  linked  array  will  be 
evaluated  in  terms  of  noise  figure,  dynamic  range, 
linearity,  signal  to  noise  and  cancellation  ratios  and 
the  antenna  patterns  of  the  arrays.  These  performance 
measures  can  be  compared  to  the  baseline  standard 


configured  four  element  array. 


FIGURE  8  -  SIX  ELEMENT  RECEIVE  ARRAY 

Design  and  development  of  an  optical  adaptive 
processor  to  be  installed  into  the  FASSP  system  is 
currently  underway.  The  goal  of  this  effort  is  to 
exploit  the  capability  of  optical  technology  in 
performing  beamforming  and  nulling.  The  optical 
adaptive  processor  will  be  a  digita  system.  It  will 
be  capable  of  performing  digital  beamforming  and 
nulling  with  greater  dynamic  range  than  previously 


developed  analog  optical  systems.  It  also  will  be 
capable  of  performing  the  adaptive  weighting  and 
processing  requirement  much  faster,  using  less  power 
than  its  electronic  counterparts  with  the  same  degree 
of  accuracy.  Design  specifications  for  the  optical 
adaptive  processor  to  be  installed  in  the  FASSP 
include  a  processing  bandwidth  of  10  MHz,  4  degrees  of 
freedom,  and  null  depths  of  40  dB  with  a  goal  for  the 
weight  update  rate  of  2.5  microseconds  (see  figure  9). 


FIGURE!  9  -  OPTICAL  LOOKUP’  TABLE 
✓ 


As  a  comparison  of  the  performance  of  an  optical 
processor  against  an  electronic  one,  it  has  been 
estimated  that  if  the  same  architecture  were  used  for 
both  types,  the  optical  processor  could  achieve 
efficiencies  of  9.15  x  10^  operations/sec.  W ,  whereas 
the  state  of  the  art  in  GaAs  would  be  limited  to  8.9  x 

o  2 

10  operations/sec.  W. 


FIGURE  10  -  OPTICAL  MULTIPLEXER  DEMULTIPLEXER 


Wavelength  Division  Multiplexing  (WDM)  is  a 
technique  where  each  signal  is  asssigned  its  own  color 
(wavelength).  The  WOM  system  is  relatively  easy  to 
implement  with  commercially  available  components.  The 


system  includes  (see  figure  10)  a  laser  diode 

transmitter  for  each  signal,  a  multiplexer  and 

demultiplexer  (these  are  typically  a  grating  which 
diffract  each  signal  at  a  different  angle  as  a 

function  of  its  wavelength),  optical  fiber,  and  a 

photodetector  receiver. 

There  are  numerous  system  applications  which 
would  benefit  by  the  WDM's  high  speed,  compact  and 
reliable  system,  and  non-radiating  and  EMI  impervious 
fiber  transmission  line.  An  immediate  application  at 
RADC  is  the  anechoic  chamber  where  the  WDM  system 
would  eliminate  the  potentially  RF  radiating,  12 
microwave  cables,  as  well  as  its  bulky  preamplifiers 
and  would  result  in  a  more  accurate  and  efficient 
system . 

WDM  cwttgmc  nf  1500  signals  are  theoreti cal  1 y 
possible  and  work  is  progressing  rapidly  toward  design 
and  realization  of  these  systems.  When  these  systems 
become  available  it  will  become  possible  to  control 
antenna  arrays  with  hundreds  of  elements  with  a  single 
WDM  system. 

6.  CONCLUSIONS 

Over  the  years  the  adaptive  test  facility  has 
been  upgraded  many  times  and  will  continue  to  be 


upgraded  in  the  future.  The  facility  has  been  a  cost 
effective  tool  for  adaptive  antenr.a  system  design. 
Development  efforts  for  beamf ormi ng/i nterf erence - 
suppression  for  satellite  communications,,  radar  and 
antenna  array  designs  have  made  use  of  this  facility 
to  reduce  design  and  development  costs. 
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A  RAPID,  ACCURATE  TEST  METHOD  FOR  MEASURING  ELECTRICAL 

PHASE  LENGTH 
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ABSTRACT 

This  paper  first  presents  an  overview  of  the  types  of  errors  present  in  phase 
length  measurements.  Next,  a  minimum-error  test  method  is  formulated.  The 
first  test  results  using  this  test  method  are  included,  demonstrating  accuracies 
on  the  order  of  ±3  electrical  degrees  at  18  GHz.  Measurement  time  using 
currently  available  ANAs  is  accomplished  in  less  than  1/2  minute  when  run  in 
the  reflective,  time-domain  mode. 

Key  words:  phase  length  measurements,  minimum-error  test  method 


1.0  INTRODUCTION 

For  optimum  performance  of  advanced  multibeam  phased-array  systems, 
the  design  criteria  of  electrical  phase  must  be  satisfied  since  any 
deviations  from  these  criteria  result  In  a  reduction  of  system  radiated 
power  (Figure  1),  beam-brpadening>  beam-steering  errors,  and  higher 
sidelobes  (Figure  2).  Because  ail  of  these  factors  are  considered 


o  20  40  GO  80  100  120  140 

a  In  Degrees 

Figure  1.  Transmitter  with  Eight  TWYs  Degradation  with  Phase  Error 
(Gaussian  Distribution) 


Phase  Error  (Degrees) 

Figure  2.  Small  Array  (N  Equals  10)  First  Sldeiobe  Level  versus 
Phase  Error 

This  paper  presents  an  overview  discussion  of  the  types  of  errors 
present  in  phase-length  measurements;  considers  the  test  options,  and 
then  formulates  a  fast,  minimum-error  test  method. 


2.0  REDUCING  ERROR  SOURCES 

Figure  3  illustrates  the  system  flow  and  errors  to  be  considered.  (Figure 
3  is  based  on  Hewlett-Packard  8510  General  Information,  pp  1-89.) 
Random  errors  are  typified  by  connector-mating  uncertainties 
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and  flexure  (bending)  effects  when  cables  are  tested  .  Machine  errors 
are  those  resulting  from  internal  errors.  They  can  be  combined  in  RSS 
fashion  to  specify  analyzer’s  proficiency.  Under  controlled  conditions, 
without  regard  to  mating  uncertainties,  this  RSS  error  represents  the  best 
tolerance  to  which  a  phase  measurement  can  be  made. 


c 


E,p  =  Sin  - 1  ( V  r/S  i  ,)  +  1fa(  phase)  +2Sh 
where 

V,  =  Sr  +  vVTx,VY?7i7 

Sr  =  Systematic  Error  =(D  +-S,,)  +  T,Sn  +(MS  +  S, , }S ,  ,2  4  (M,  i  S,2> 

S21S12  +  AS,, 

Wr  =  Random  low-level  noise  =  3N, 

Xr  =  Random  high-level  noise  =  3NhSn 

Y,  =  Random  Port-1  repeatability  =  R(,  i  2Rt)Sl1  ,  R,,Sn? 

2r  =  Random  Port-2  repeatability  =  R,2S2iS,2 

T[d  =  Phase  stability  with  temperature  =  |0  1  r  001  ♦  IGh^]  C 

S„  =  Port  1  cable  transmission  stability  =  (0  05)  x  t^H/ 

Figure  3.  System  Flow  and  RSS  Errors 

Systematic  errors  ore  typified  as  bias  effects.  For  example,  for  a  phase 
measurement  relying  on  comparison  to  a  reference  channel  cable,  the 
difference  in  temperatures  between  the  test  and  reference  channels 
produces  a  bias  error.  In  addition,  a  bias  error  will  exist  in  converting 


electrical  length  to  physical  length  if  the  velocity  of  light  used  in  the 
computation  is  not  corrected  for  local  ambient  conditions. 

The  correction  for  calculating  the  speed  of  light  relative  to  local  ambient 
conditions  requires  knowledge  of  temperature,  humidity,  and  air 
pressure.  Given  these  values,  and  using  an  NBS-furnished  equation,  a 
correction  can  be  effected. 

Connector-mating  uncertainties  can  be  reduced  by  using  precision 
connectors  and  adapters.  Flexure  uncertainties  in  cabie  tests  can  be 
avoided  using  a  cable-holding  jig. 

Using  a  frequency  synthesizer  locked  to  an  ultra-stable  source  minimizes 
the  frequency  error  so  that  the  remaining  error  is  associated  to  the  errors 
in  determining  time  delay.  It  has  been  found  that  the  best  results  are 
achieved  by  linearizing  the  network  analyzer  over  the  widest  bandwidth 
(i.e.,  2  to  18  GHz  is  typical). 

Following  linearization  and  transformation,  the  machine  is  commanded 
to  measure  the  distance  difference  between  reflection  from  a  short  at  the 
output  and  then  from  the  same  short  positioned  at  the  end  of  the  cable- 
under-test. 


Errors  in  the  resulting  measurement  include  the  intrinsic  network 
analyzer  error  and  the  uncertainty  errors  associated  with  connecting  to 
the  cable-under-test. 


3.0  TEST  METHOD 

Looking  at  the  two  choices  for  measuring  electrical  length,  it  appears  that 
less  error  is  Incurred  in  a  reflection-type  measurement  (see  Table  1). 
Figure  4  shows  a  simplified  schematic  of  the  test  configuration  used  to 
automatically  determine  the  electrical  length  of  the  cable.  Either  network 
analyzer  can  be  programmed  Kind  used  to  find  the  time  delay  to  a  short  at 
the  end  of  the  cable-under-test.  Internal  to  the  network  analyzer,  the 
measurement  is  performed  in  the  frequency  domain  and  the  value  is  then 
converted  to  time  delay.  A  final  computation  is  made  to  determine 
electrical  length, 

Electrical  Length  s  360  x  F  x  T  (1 ) 

where: 

F  =  Frequency  in  GHz 

T  =  One-way  cable  delay  in 
nanoseconds 
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TABLE  1.  MEASUREMENT  CHOICES 


Transmission  Reflection 

Requires  the  fewest  X 

number  of  adapters 


Permits  measurement  with-  X 

out  distorting  the  cable 


Transits  through  the  cable  X 

N  =  2,  averaging  effect 


(a)  Test  Schematic  (b)  Illustrative  Measurement 


Figure  4.  Measurement  Concept 


Rather  than  attempting  to  quantify  each  error  and  computing  an  RSS 
error,  the  following  method  was  used: 

1.  First,  using  beadless  air  lines,  the  mechanical  length  is  accurately 
found. 

2.  Next  the  electrical  length  is  computed  from  these  dimensions, 
taking  into  account  local  ambient  temperature,  humidity  and  pressure 
which  slightly  affect  the  speed-of-light  to  be  used  in  computing  length. 

3.  The  air  line  is  then  measured  using  the  network  analyzer. 

4.  Results  are  then  compared. 

5.  The  connection  and  measurement  (steps  2  through  4)  are  ihen 
repeated  several  times  to  determine  a  quasi-probable  error  (many  more 
measurements  over  a  time  span  are  required  to  obtain  a  statistical 
probable  error). 

Following  the  above  calibration,  the  electrical  length  for  the  cable  is 
measured,  first  from  J1  to  J2  and  then  from  J2  to  J1.  Nex\,  an  average 
electrical  length  is  computed.  Following  the  test,  the  calibration 
intactness  is  re-verified  with  the  air  line. 


Two  sets  of  results  are  given. 


V 


4.1  AIR  LINE  MEASUREMENTS 

Following  the  above  procedure,  air  lines  were  tested  at  both  Hewlett- 
Packard  and  Wiltron.  First  the  mechanical  lengths  were  determined. 
Next  the  lengths  were  determined  using  the  respective  ANAs.  Values 
found  by  each  ANA  measurement  were  then  corrected  for  local  ambient 
conditions.  Table  2  shows  the  results. 


TABLE  2.  AIR  LINE  RESULTS 


Vendor 

e*: _ * 

r  ioi 

Article 

(Beadless 

Air  Lines) 

mcasuicu 

Mechanical 

Length 

(MM) 

MedSUrdu 

Length 

By  ANA 
(MM) 

Corrected 

ANA 

Measurement 

(MM) 

A 

7.5  CM 

74.930  ±0.013 

75.064 

75.044 

B 

10.0  CM 

100.012  ±0.013 

100  040 

100.008 

■  Vendor  A  results  =  0.11-'.  milli-meters  or  about  2.5  electrical  degrees  at  18  GHz 


■  Vendor  B  results  =  0.028  milli-meters  or  about  0.6  electrical  degrees  at  18  GHz 


■  Repeat  test  for  vendor  A  showed  1.5'’  at  18  GHz 
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The  difference  for  Vendor  A  showed  0.114  mm,  which  is  about  2.5 
electrical  degrees  at  18  GHz.  The  difference  for  Vendor  B  showed  0.028 
mm,  which  is  about  0.6  electrical  degree  at  18  GHz.  In  addition,  similar 
measurements  were  performed  by  a  cable  manufacturer  using  Vendor 
A's  equipment.  The  resulting  error  was  1.5  electrical  degrees  at  18  GHz. 

An  interesting  fact  is  that  for  stainless  steel  airlines  of  a  nominal  length  of 
7.5  cm,  a  temperature  change  of  ±3°C  changes  the  length  about  ±0.i 
electrical  degree  at  18  GHz. 


4.2  CABLE  MEA3VPSMSNTS 

Following  the  air  line  tests,  the  electrical  lengths  for  four  0.141  mm  semi¬ 
rigid  cables  were  determined  by  each  manufacturer.  Table  3  shows  the 
results. 

From  this  tabulation  it  may  be  observed  that  the  electrical  lengths  agree 
within  6  electrical  degrees,  or  approximately  0.1  percent.  The  final 
operation  is  to  find  the  mechanical  length  for  the  above  cables. 


ft 


TABLE  3.  CABLE  RESULTS 


Vendor 

Test  Article 
(0.141  Semi- 
Rigid  Straight 
Cables) 

Measured 
Time  Delay 

In  Pico 
Seconds 

ANA  Calculated 
Electrical  Length 

In  Degrees 

At  18  GHz 

A 

1 

1665  053 

5395 

A 

2 

1667.098 

5403 

A 

0 

1667.368 

5402 

A 

i  /■> 

*c.c-y.  ton 

IVUv  Vv 

r.oov 

1  B 

1 

1 606  64  4 

5409 

B 

2 

1667  644 

5403 

0 

5 

1868.644 

5406 

B 

6 

1607.644 

5403 

Observe  mat  the  electrical  lengths  agree  within  6''  or 
approximately  0.1  percem 


The  following  expression  relates  the  /arables, 

Electrical  Length  =  360°  x  L  \  x  Ve 
where: 

L  \  Physical  Length 
x  =  Wave  Length  =  C/F 
e  -  Dielectric  Constant 
C  =  Velocity  of  Light 
F  =  Frequency 
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In  using  equation  (2)  both  x  and  L  must  be  in  consistent  units,  and 
accurate  knowledge  of  the  dielectric  constant  (as  it  exists  in  the  cable)  is 
required.  This  can  be  found  by  testing  a  long,  known  length  of  cable  and 
computing  the  phase  velocity,  which  is  1/Ve.  In  addition,  the  speed  of 
light  should  be  corrected  for  air  temperature,  humidity,  and  pressure. 


For  example,  in  a  vacuum,  the  speed  of  light  is  given  as: 

C  =  2.997928  ±0.000004  x  1010  CM/Second 
(Physics  and  Chemistry  Handbook) 

Using  an  NBS  computation,  the  effective  dielectric  constant  of  air  is  found 
to  be  1.000649  for  760  rr.m  pressure,  50  percent  humidity,  and  20°C. 

Thus  the  speed  of  light  used  in  equation  (2)  should  be  multiplied  by, 

-1-  s  0.9996757 

VeAIR 

Finally,  the  frequency  must  be  accurately  known,  and  frequency  stability 
maintained  throughout  tests. 


5.0  CONCLUDING  REMARKS 


1.  A  reliable  and  fast  test  method  for  accurately  determining 
electrical  length  has  been  described. 

2.  By  rigorously  following  the  procedure,  the  electrical  lengths  can  be 
found  to  within  ±3  electrical  degrees,  provided  connector  uncertainties 
do  not  exceed  ±2  electrical  degrees. 

3.  In  tests  with  both  analyzers,  it  was  found  that  the  phase  length 
measurement  time  was  under  1/2  minute. 
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ABSTRACT 


Because  of  the  long  wavelength  of  microwave  radiation,  the  imaging 
of  targets  requires  very  large  antennas.  This  paper  describes  three 
methods  of  imaging  aircraft  from  the  ground  with  microwave  radar 
with  quality  suitable  for  aircraft  target  recognition.  The  imaging 
methods  are  based  upon  a  self-calibration  procedure  called  adaptive 
beamforming  (ABF)  that  compensates  for  the  severe  geometric  distor¬ 
tion  inherent  in  any  imaging  system  that  is  large  enough  to  achieve 
the  high  angular  resolution  necessary  for  two-dimensional  target 

nudging*  itit;  OlgUCXJL  pi  wcooill^  luiliu  lo  cui*^  /v  baiiu 

cm)  wavelength  experiments  demonstrate  its  success  on  commercial 
aircraft  flying  into  Philadelphia  International  Airport. 


1 .  INTRODUCTION 


This  paper  describes  three  methods  of  imaging  aircraft  from  the 
ground  with  microwave  radar  with  quality  suitable  for  aircraft 
target  recogi  it  ion.  Illustrations  are  shown  in  the  upper  portions 
of  Figure  1-3*  Plan-view  drawings  of  the  airplanes  are  included  to 
facilitate  comparison.  The  imaging  methods  are  based  upon  a  self- 
calibration  procedure  called  adaptive  beamforming  (ABF)  that 
compensates  for  the  severe  geometric  distortion  inherent  in  any 
imaging  system  that  is  large  enough  to  achieve  the  high  angular 
resolution  necessary  for  two-dimensional  target  imaging. 

The  lower  portions  of  the  figures  illustrate  the  failure  of 
conventional  electronic  scanning  to  form  useful  images;  in  all  three 
cases  the  distortions  ir:-  the  microwave  apertures  grossly  exceed  the 
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tolerances  stipulated  by  diffraction  theory*  The  same  data  sets  were 
used  to  form  the  upper  and  the  lower  images  in  each  experiment.  ABF 
was  introduced  only  when  the  upper  images  were  being  formed.  It  is 
evident  from  the  ease  of  target  recognition  from  these  images  that 
the  procedure  phase -corrected  the  radar  echos  and  thereby  compen¬ 
sated  for  the  distortions  in  the  apertures.  The  result  (upper)  is 
high  resolution,  diffraction-limited  imagery  from  which  length  and 
wing  span  can  be  measured,  the  sweep  of  the  wings  observed,  the 
question  of  whether  or  not  there  are  engines  in  the  wings  answered, 
and  so  forth.  In  short,  target  classification  becomes  possible. 

2.  ADAPTIVE  BEAMFORMING 

ABF  is  a  data  adaptive  feedback  control  procedure  that  adjusts 
the  complex  weight  vector  in  a  large  phased  array  so  as  to  compen¬ 
sate  for  geometric  and  electrical  distortion  in  the  system.  The 
control  information  is  derived  from  measurements  of  the  complex 
radiation  field  from  the  targets  being  imaged. 

2.1  Dominant  Scatterer  Algorithm 

The  basic  procedure  assumes  that  there  exists  somewhere  in  the 
field  of  view  of  the  imaging  system  a  point-like  scatterer  cr  source 
having  large  radar  cross-section  or  source  strength.  A  corner 
reflector  proves  to  be  an  excellent  target  for  the  beamforming 
procedure.  The  theory  governing  the  requirements  on  its  characteris- 
tics  is  published. ’ c .  There  it  is  shown  that  the  physical  size  of 
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the  source  can  be  no  larger  than  AR/2L,  where  A  is  the  wavelength..  R 
is  the  raaar-to-target  distance  and  L  is  the  size  of  the  imaging 
aperture;  and  that  its  echo  strength  must  exceed  the  total  backseat  - 
ter  from  all  the  other  scatterers  in  its  range  bin  by  at  least  4  dB. 
Such  a  target  or  source  radiates  a  nearly  spherical  wavefront,  which 
induces  a  simple  phase  variation  across  the  array.  Provided  that  the 
array  is  linear  or  planar,  the  phase  variation  is  linear  if  the 
target  or  the  source  is  in  the  far  field,  or  approximately  quadratic 
when  the  source  is  in  the  near  field.  Deviations  from  such  simple 
behavior  indicate  geometric  distortion  in  the  array,  electrical 
mistunmgs,  or  wavefront  distortion  due  to  turbulence  in  the 
propagation  medium. 

2 

Figure  4  models  the  problem  .  A  transmitter  in  the  upper  right 
launches  a  short  RF  pulse  toward  a  point  3Catterer.  The  reradiation 
illuminates  the  receiving  phased  array  with  a  spherical  wavefront. 
The  array  is  shown  distorted  .in  two  dimensions:  the  antenna  elements 
are  displaced  from  a  datum  line  denoting  the  location  of  the  array, 
and  the  spacings  between  the  elements  are  not  equal.  Thus  tne  phases 
of  the  signals  received  by  these  elements  are  not  those-  determined 
by  the  distances  from  tne  point  target  to  the  intersections  of  a 
periodically  sampled  straight,  line  (the  array)  and  a  circle  (cross 
section  of  the  spherically  expanding  wave).  Instead,  the  measured 
phases  are  altered  by  the  displacements  of  the  antenna  elements  from 
their  design  positions.  The  phase-error  tolerance  is  smaller  than 
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one-tenth  of  a  wavelength  if  the  gain  of  the  array  is  to  be  main¬ 
tained  to  within  about  one  dB,  and  smaller  still  if  reasonable 

7 

control  is  to  be  exercised  over  the  side  radiation  pattern  . 
The  tolerable  distortion,  therefore,  is  very  small.  Similarly  small 
are  the  allowed  mistunings  in  the  receiver  channels  (not  shown) 
behind  each  antenna  element  or  the  phase  distortion  in  the  wavefront 
induced  by  variations  of  the  refractive  index  of  the  propagation 
medium. 

Adaptive  beamforming  compensates  for  this  problem  in  the  fol¬ 
lowing  manner:  The  amplitudes  of  the  echoes  from  the  point  target, 
as  measured  in  the  receivers,  are  little  perturbed  by  the  geometric 
errors  in  the  array  or  the  electrical  errors  in  the  receivers.  Phase 
information,  on  the  other  hand,  is  substantially  destroyed  when  the 
displacements  of  the  antennas  exceed  a  small  fraction  of  the 
wavelength.  The  signal  processor  uses  the  fact  that  the  variance  of 
the  amplitude  of  the  echoes  from  the  point  target,  as  measured  in 
the  receivers,  is  small:  it  searches  in  range  for  that  range  bin, 
called  the  reference  range,  in  which  the  echoes  have  the  smallest 
amplitude  variance.  Corrections  are  easily  made  to  the  array  by 
phase  shifting  all  the  signals  received  from  the  point  scatterer  at 
the  reference  range  to  a  common  phase.  This  process  is  called  phase 
conjugation  ’  and  is  schematically  indicated  by  the  feedback  net¬ 
work  in  the  lower  right  of  the  figure.  A  superior  correction  is  a 
complex  'weight  at  each  location  that  is  the  complex  conjugate  of  the 
received  signal.  Inis  correction  results  in  a  radiation  pattern  that 


is  approximately  the  matched  spatial- filter  response  of  the  system 
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to  the  target  or  source  ’  .  This  means  that  the  magnitude  of  the 
radiation  pattern  closely  follows  the  modulus  of  the  reflection 
profile,  in  the  transverse  or  azimuth  direction,  of  the  point-like 
scatterer.  In  this  way,  a  very  narrow  beam  is  formed  from  a  large, 
distorted  antenna  array  provided  that  the  adaptive  beamforming 
target  driving  the  process  reradiates  like  a  point  source,  A  corner 
reflector  has  been  found  to  emulate  such  a  target  very  well. 

Table  1  describes  the  algorithm.  Following  compensation  for 
known  gain  variations  in  the  antenna  elements  (step  2),  the  system 
searches  through  range  for  that  range  tin  in  which  the  normalised 
variance  of  the  amplitudes  of  the  echoes  across  the  array  is  the 
smallest  (step  3).  It  designates  that  range  bin  as  the  reference 
bin  and  calls  its  range  R  .  It  then  conjugates  the  phases  at  kQ  to 
form  the  weight  vector,  phase  shifts  ail  eences  by  the  weight  vector 
(steps  A  and  5),  focuses  the  array  at  all  ranges  (step  6),  and  scans 
the  array  in  azimuth  to  form  the  two-dimensional  image  (steps  7  ana 
8).  This  algorithm  is  sometimes  called  the  minimum  variance  algo¬ 
rithm  (MVA)  arid  sometimes  the  dominant  scot"  e’v-r  algorithm  (DSA). 

2.2  Reference  Target  of  Opportunity 

In  imaging  aircraft  from  the  ground,  it  is  net  possible  to 

implant  an  artificial  target  such  as  a  corner  reflector  in  the 
general,  area  of  tne  target.  Nor  can  it  be  assumed  that  point-like 
targets  of  opooroun.  .y  will  be  in  the  neighborhood  of  the  airplane. 
Hence  the  procedure  d  >per<j3  upon  tne  existence  of  one  or  more 
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TABLE  1.  STEPS  IN  RADIO  CAMERA  IMAGING  (FROM  [27], [23]) 


Step 

1  Mcasuic  and  store  complex  en¬ 
velopes  of  echo  samples 

2  Correct  amplitudes  by  dividing 
by  element  pattern  estimate  /„ 

3  Find  R0  such  that  A0n  «  A, 
all  n 

4  Phase  rotate  at  R0  by  phase 
conjugate  in  relation  to  refer¬ 
ence  dement, exp  j(if>un  ~  t on) 

5  Phase  rotate  at  all  range  ele¬ 
ments 

6  Focus  at  each  range  R , 

7  Phase  shift  linearly  with  angle 

8  Sum  at  each  range  element 


range  bin— k1  ^ — element  number 


Ae'*0’ 

Ae‘*m 
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suitable  reflectors  on  the  aircraft  itself. 

It  is  evident  that  the  smaller  the  range  cell,  the  larger  the 
probability  that  a  range  bin  containing  a  single  dominant  reflector 
exists.  100  experiments  on  approximately  30  aircraft  were  conducted 
to  determine  the  statistics.  Data  we  re  taken  at  X-band  (1  =  3  cm) 
at  the  Valley  Forge  Research  Center  of  the  University  of 
Pennsylvania.  Hie  aircraft,  flying  at  distances  from  3  to  30  km 
from  the  laboratory,  were  primarily  commercial  airplanes  flying  into 
the  Philadelphia  International  Airport;  some  were  military  aircraft 
associated  with  the  Naval  Air  Station  in  Willow  Grove,  PA.  When  the 
range  cell  was  1  m,  a  suitable  beamforming  source  was  available  in 
approximately  80%  of  the  experiments.  When  the  range  cell  was 
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increased  to  3  m,  the  success  rate  dropped  to  approximately  25%  .  A 
1  m  range  ceil  was  employed  in  the  imagery  shown  in  Figures  1-3- 


2.3  Efficacy  of  Adaptive  Beamforming 

Figures  1 '3  show  the  crucial  importance  of  adaptive  beamforming 
in  large  imaging  systems.  Their  lower  portions  show  imagery  from  the 
same  data  sets  without  the  self-calibration  procedure.  The  image  of 
each  such  unphase-corrected  echo  sequence  is  a  hopeless  jumble  of 
pixels. 

3.  SPATIAL  RADIO  CAMERA 

g 

The  term  "radio  camera"  was  introduced  to  capture  the  idea  of 
a  picture-taking  instrument  at  radio  or  radar  wavelengths.  It  is 
characterized  by  an  aperture  large  enough  to  provide  the  desired 
azimuthal  or  transverse  resolving  power  but,  because  of  its  large 
size,  is  intrinsically  distorted  beyond  the  tolerance  permitted  by 
diffraction  theory.  Compensation  is  provided  by  the  adaptive  beam- 
forming  procedure. 

3.1  32-Receiver  Experiment 

A  bistatic  radar  configuration  was  employed  in  the  experiment 
that  produced  the  aircraft  image  of  Figure  1 .  The  X-band  transmitter 
radiated  a  1  kw  peak -power  pulse  of  7  ns  duration.  The  range  cell, 
therefore,  was  1  m,  and  the  radiated  and  received  bandwidth  was  150 
MHz.  The  antenna  was  a  1  .2  ,i  paraboloid  mounted  on  a  mechanically 
steerable  pedestal.  The  receiving  array  consisted  of  32  X-band 
receivers  deployed  on  a  laboratory  roof  in  an  irregular  line.  The 
receiver  spacings  were  approximately  1  m.  The  size  of  array,  there¬ 
fore,  was  approximately  1000  times  the  wavelength.  Diffraction 


theory  teaches  that  the  beamwidth  cf  a  diffraction-limited  antenna 

is  (in  radians)  about  the  reciprocal  of  its  size  as  measured  in 
3  1 0 

wavelengths'  ’  .  After  ABF,  v/hioh  converted  the  system  from  a 

highly  irregularly  placed  set  of  antennas  to  a  diffraction-limited 
phased  array,  the  radiation  pattern  that  was  formed  exhibited  a 
beamwidth  of  10  rad  or  1  mr,  as  expected.  At  the  nominal  target 
distance  of  3  km,  the  transverse  dimension  of  the  resolution  cell 
was  3  m.  The  receiver  cabinets  were  propped  up  on  cinder  blocks  to 
an  elevation  of  approximately  30°.  Each  cabinet  had  a  horn  antenna 
19  x  1 M  cm  and  a  horizontal  beamwidth  of  12°.  The  pointing  angles 
and  pattern  orientations  uii  iereu  somewhat  from  receiver  to 
receiver. 

Each  antenna  in  the  receiving  array  received  a  modulated  RF 
wave  consisting  of  the  echoes  of  targets  and  clutter  from  the  trans¬ 
mitted  pulse.  The  receivers  were  fully  coherent  in  the  sense  that 
all  oscillators  were  frequency-and  phase-locked  to  a  120  MHz  master 
oscillator  in  the  system.  Its  waveform,  is  frequency-multiplied 
several  times  to  provide  the  transmitter  RF  and  the  required  local 
oscillator  waveforms.  The  received  RF  waveforms  were  detected  in 
each  receiver,  coherently  quadrature  demodulated,  sampled  at  a  200 
megasample  per  second  rate  and  each  quadrature  sample  was  converted 
to  an  8-bit  digital  word. 

3.2  Signal  Procct'-ing 


The  sampled  and  stored  data  are  the  complex  field  amplitudes 
received  by  the  antennas  and  measured  in  the  receivers.  The  function 
of  the  sign  processor  is  to  transform  these  data  into  an  image. 
The  large  size  of  the  aperture  precludes  the  possibility  that  it  is 
undistorted.  As  a  consequence,  the  argument  of  each  measured  complex 
field  amplitude  is  in  error  by  the  altered  phase  delay  from  scat- 
terer  to  antenna  due  to  the  misplacement  of  each  antenna  from  its 
design  position.  The  adaptive  beamforming  algorithm  described  ear¬ 
lier  corrects  this  error.  Without  such  correction  the  measured 
phases  are  garbled  and  the  resulting  image  is  a  jumble  of  pixels. 
The  lower  picture  in  Figure  1  i3  an  example. 

H.  TEMPORAL  RADIAL  CAMERA  (ISAR  +  A3F) 

1  i  -|  ?  1  g  1  9 

Synthetic  aperture  radar  (SAR)  ’  was  a  monumental 
invention  of  the  early  1950's.  Deviations  from  nonaccelerating 
flight  due,  for  example,  to  wind  gusts  distort  the  synthetic 
aperture.  Measurements  can  be  made  within  the  airplane  with  pen¬ 
dulums  or  other  inertial  devices  to  note  the  instantaneous 
acceleration  of  the  platform,  from  which  the  instantaneous  vector 
deflection  from  straight  line  flight  can  be  calculated.  In  turn, 
these  measurements  can  lead  to  phase  corrections  made  in,.the  local 
oscillator  of  the  receiver  or  in  the  data  set  to  correct  for  the 
aperture  distortions.  In  this  way,  SAR  systems  are  able  to  compen¬ 
sate  for  deviations  from  straight  line,  constant  speed  flight. 


Not  so  with  inverse  synthetic  aperture  radar  (ISAR),  however. 
Although  the  concept  was  a  product  of  the  3ame  decade,  the  problem 
of  phase  compensation  is  intrinsically  more  complicated  and  has 
delayed  ISAR's  application.  In  contrast  to  SAR  the  aircraft  is  the 
target  and  the  radar  is  on  the  ground  or  on  some  arbitrary  platform. 
The  target  motion  is  used  to  synthesize  an  extended  aperture  which 
can  then  image  the  airplane.  As  a  consequence,  the  radar  has  no 
means  for  obtaining  the  mechanical  measurements  of  instantaneous 
aircraft  position  and  correcting  the  radar  data  accordingly.  This  is 
exactly  what  adaptive  beamforming  can  do,  however.  Thus  radio  camera 
processing  is  directly  applicable  to  ISAR.  Time  replaces  the  spatial 
dimension.  That  is,  the  radar  waveform  received  as  a  result  of  the 
nth  transmission  (of  N)  is  treated  in  the  signal  processor  as  if  it 
were  the  nth  antenna  element  in  an  array  of  N  antennas. 

Figure  2  shows  an  airplane  imaged  in  this  fashion.  The  aircraft 
is  a  Boeing  727  flying  at  3.2  km  from  VFRC.  The  same  antenna  and 
transmitter  described  in  the  earlier  section  were  used  for  this  ISAR 
experiment.  Only  a  single  receiver  was  employed,  however,  and  it  was 
driven  by  the  same  1.2  m  dish  used  for  the  transmitter.  The  radar 
was  operated  in  a  conventional  monostatic  fashion.  The  pulse  length, 
as  before,  was  1  m.  The  lower  portion  of  the  figure  shows  the  image 
from  the  same  data  sequence  without  the  phase  corrections  introduced 
by  adaptive  beamforming.  As  before,  the  image  is  useless  when  phase 
compensation  for  the  geometric  distortion  is  not  applied. 


5.  SPATIAL- TEMPORAL  RADIO  CAMERA  (COMBINED  ISAR  AND  RADIO  CAMERA) 

Because  the  same  signal  processing  is  applied  to  both  the 
spatial  and  the  temporal  systems,  it  is  not  unreasonable  to  assume 
that  they  may  be  combined  into  a  common  spatial-temporal  imaging 
system. 

Figure  5  illustrates  the  procedure  with  two  receivers.  A  trans¬ 
mitter  is  colocated  with  one  receiver.  The  target  is  an  airplane. 
Synthetic  apertures  are  formed  by  each  receiver.  The  combined  data 
from  the  two  antennas  form  the  total  data  set.  Adaptive  beamforming 
compensates  for  flight  path  deviations  in  each  of  the  synthetic 
apertures,  and  also  compensates  for  the  unknown  jump  in  phase  be¬ 
tween  the  two  data  sequences.  The  upper  image  of  Figure  3  was  formed 
in  this  manner.  The  two  receiving  antennas  were  25  m  apart,  the 
aircraft  was  flying  at  approximately  250  knots  transverse  to  the 
line-of-sight,  the  interpulse  period  was  2.5  ms  and  the  number  of 
pulses  integrated  was  85.  The  radar  wavelength  was  3  cm. 

7.  SUMMARY 

Microwave  raaar  has  now  matured  to  the  point  where  aircraft 
free  of  a  clutter  background  can  be  imaged  in  two  dimensions  with 
sufficient  resolution  to  recognize  the  target.  While  the  experi¬ 
ments  were  performed  from  tne  ground,  there  is  nothing  in  trie  nature 
of  the  techniques  described  in  this  paper  to  preclude  their  ap¬ 
plicability  to  airborne  or  space  platforms.  High  resolution  in  the 
range  dimension  requires  a  large  bandwidth.  For  example,  15C  MHz  is 
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needed  to  achieve  one  meter  in  range.  Such  bandwidth  can  be 
achieved  with  existing  design  procedures  and  components. 

The  same  is  not  true  in  the  orthogonal  dimension,  however. 
High  angular  or  cross-range  resolution  requires  an  aperture,  so  huge 
as  to  be  intrinsically  distorted,  often  making  it  impossible  to  form 
a  highly  directive  main  lobe. 

Electrical  distortion  (primarily  phase  errors)  will  accrue 
through  mistunings,  impedance  mismatches  and  component  drift  with 
age  and  temperature.  A  self-calibrating  technique  called  adaptive 
beamforming  (ABF)  successfully  overcomes  such  distortion  and  cor¬ 
rectly  reconstructs  the  main  beam  properties.  Such  an  instrument  is 
called  a  radio  camera.  The  same  procedure  compensates  for  flight- 
path  distortion  in  airborne  targets  that  are  in  the  clear,  thereby 
permitting  ISAR  imaging  to  work  satisfactorily.  Because  the  signal 
processor  is  the  same  in  the  two  systems,  the  combination  of  ISAR 
and  ABF  is  called  a  temporal  radio  camera.  Further,  because  the 
same  adaptive  beamforming  procedure  solves  both  problems,  they  may 
be  combined  to  form  highly  thinned,  huge  imaging  arrays.  Figures  1  - 
3  demonstrate  the  quality  currently  being  achieved  by  these  three 
imaging  configurations. 
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FIGURE  1.  MICROWAVE  HUGE  OF  BOEING  727  OBTAINED 
WITH  32m  X-BAND  RADIO  CAMERA.  AIRPLANE 
IS  FLYING  AT  DISTANCE  OF  3.2  Km  FROM 
VALLEY  FORCE  RESEARCH  CENTER. 
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FIGURE  2.  MICROWAVE  IMAGE  OF  DOE  1  NO  72  7  OBTAINED  AT  X-BAND 
WITH  MONOSTATIC  RADIO  CAMERA.  AIRPLANE  IS  FEY  INC 
AT  DISTANCE  OF  2.7  Kim  FROM  YAEI.EY  FORGE  RESEARCH  CENTER 
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FIGURE  3.  MICROWAVE  IMAGE  OF  LOCKHEED  L-1011  OBTAINED  WITH  x-BAND 
RADIO  CAMERA  CONSISTING  OF  TWO-ELEMENT  INTERFEROMETER 
PLUS  INVERSE  SYNTHETIC  APERTURE  PROCESSING.  THE 
AIRPLANE  IS  FLYING  AT  A  DISTANCE  OF  2.5  Km  FROM 
VALLEY  FORGE  RESEARCH  CENTER. 
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FIGURE  4.  THE  DISTORTED  ARRAY  SAMPLES  THE  WAVEFRONT  AT  IRREGULAR  INTERVALS. 

MEASURED  PHASES  ARE  IN  ERROR  BECAUSE  OF  DISPLACEMENTS  OF  ANTENNA 
ELEMENTS  FROM  THEIR  DESIGN  POSITIONS.  ADAPTIVE  BEAMFORMING  COM¬ 
PENSATES  FOR  THESE  ERRORS. 
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SPATIAL-TEMPORAL  RADIO  CAMERA.  TWO  OR  MORE  RECEIVERS  FORM 
SYNTHETIC  APERTURES  SIMULTANEOUSLY.  ADAPTIVE  BEAMFORMING 
PERMITS  THE  DATA  TO  BE  COMBINED  TO  REMOVE  GRATING  LOBES 
AND  PROVIDE  FINER  RESOLUTION.  IN  ADDITION,  THE  INDIVIDUAL 
IMAGES  ARE  COMBINED  FOR  DIVERSITY  ENHANCEMENT. 
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RAVC  plant  and  executes  research,  de.ve.lo pme.nl,  test 
and  selected  acquisition  programs  In  support  of 
Command,  Control,  Communications  and  Intelligence 
( C 3 1)  activities.  Technical  and  engineering 
support  within  areas  oft  competence  ts  provided  to 
ESV  Program  Offices  (POs)  and  other  ESV  elements 
to  perform  elective  acquisition  of  C^I  systems. 

The  areas  o technical  competence  Include 
communications,  command  and  control,  battle 
management,  Information  processing ,  surveillance 
sensors,  Intelligence  data  collection  and  handling, 
solid  state  sciences,  electromagnetics ,  and 
propagation,  and  electronic,  maintainability , 
and  compatibility. 


